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PREFACE 


This book embodies a course given by the writer for a number of 
years in the Mathematical Laboratory of the Massachusetts Institute 
of Technology. It is designed as an aid in the solution of a large num- 
ber of problems which the engineer, as well as the student of engineering, 
meets in his work. 

In the opening chapter, the construction of scales naturally leads to 
a discussion of the principles upon which the construction of various 
slide rules is based. The second chapter develops the principles of a 
network of scales, showing their application to the use of various kinds 
of codrdinate paper and to the charting of equations in three variables. 

Engineers have recognized for a long time the value of graphical 
charts in lessening the labor of computation. Among the charts devised 
none are so rapidly constructed nor so easily read as the charts of the 
alignment or nomographic type —a type which has been most fully 
developed by Professor M. d’Ocagne of Paris. Chapters III, IV, and V 
aim to give a systematic development of the construction of alignment 
charts; the methods are fully illustrated by charts for a large number 
of well-known engineering formulas. It is the writer’s hope that the 
simple mathematical treatment employed in these chapters will serve to 
make the engineering profession more widely acquainted with this time 
and labor saving device. 

Many formulas in the engineering sciences are empirical, and the 
value of many scientific and technical investigations is enhanced by the 
discovery of the laws connecting the results. Chapter VI is concerned 
with the fitting of equations to empirical data. Chapter VII considers 
the case where the data are periodic, as in alternating currents and volt- 
ages, sound waves, etc., and gives numerical, graphical, and mechanical 
methods for determining the constants in the equation. 

When empirical formulas cannot be fitted to the experimental data, 
these data may still be efficiently handled for purposes of further 
computation, — interpolation, differentiation, and integration, — by the 
numerical, graphical, and mechanical methods developed in the last 
two chapters. 

Numerous illustrative examples are worked throughout the text, 
and a large number of exercises for the student is given at the end of 
each chapter. The additional charts at the back of the book will serve 
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as an aid in the construction of alignment charts. Bibliographical 
references will be found in the footnotes. 

The writer wishes to express his indebtedness for valuable data to 
the members of the engineering departments of the Massachusetts 
Institute of Technology, and to various mathematical and engineering 
publications. He owes the idea of a Mathematical Laboratory to 
Professor E. T. Whittaker of the University of Edinburgh. He is 
especially indebted to Capt. H. M. Brayton, U.S. A., a former student, 
for his valuable suggestions and for his untiring efforts in designing a 
large number of the alignment charts. Above all he is most grateful to 
his wife for her assistance in the revision of the manuscript and the 
reading of the proof, and for her constant encouragement which has 
greatly lightened the labor of writing the book. 


JOSEPH LIPKA. 


CAMBRIDGE, MAss., 
Oct. 13, 1918. 
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CHAPTER VI. 
EMPIRICAL FORMULAS — NON-PERIODIC CURVES. 


68. Experimental data. — In scientific or technical investigations we | 
are often concerned with the observation or measurement of two quanti- 
ties, such as the distance and the time for a freely falling body, the volume 
of carbon dioxide dissolving in water and the temperature of the water, 
the load and the elongation of a certain wire, the voltage and the current 
of a magnetite arc, etc. The results of a series of measurements of the 
same two quantities under similar conditions are usually presented in the 
form of a table. Thus the following table gives the results of observa- 
tions on the pressure p of saturated steam in pounds per sq. in. and the 
volume v in cu. ft. per pound: 

p= 10 20 30 40 50 60 
v= 337.80 19.72 13.48 10.29 8.34 6.62 


We represent these results graphically by plotting on codrdinate paper 
the points whose codrdinates are the corresponding values of the measured 
quantities and by drawing a smooth curve through or very near these 
points. Fig. 68 gives a graphical representation of the above table, 
where the values of # are laid off as abscissas and the values of v as ordi- 
nates and a smooth curve is drawn so as to pass through or very near the 
pletted points. 

I20 
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The fact that a smooth curve can be drawn so as to pass very near the 

plotted points leads us to suspect that some relation may exist between 
the measured quantities, which may be represented mathematically by 
the equation of the curve. Since the original measurements, the plotting 
of the points, and the drawing of the curve all involve approximations, 
the equation will represent the true relation between the quantities only 
approximately. Such an equation or formula is known as an empirical 
ormula, to distinguish it from the equation or formula which expresses 
a physical, chemical, or biological law. A large number of the formulas 
in the engineering sciences are empirical formulas. Such empirical for- 
mulas may then be used for the purpose of interpolation, 7.e., for comput- 
ing the value of one of the quantities when the value of the other is given 
within the range of values used in determining the formula. 

It is at once evident that any number of curves can'be drawn so as to 
pass very near the plotted points, and therefore that any number of 
equations might approximate the data equally well. The nature of the 
experiment may give us a hint as to the form of the equation which will 
best represent the data. Otherwise the problem is more indeterminate. 
If the points appear to lie on or near a straight line, we may assume an 
equation of the first degree, y = a + bx, in the variables. But if the 
points deviate systematically from a straight line, the choice of an equa- 
tion is more difficult. Often the form of the curve will suggest the type 
of equation, parabolic, exponential, trigonometric, etc., but in all cases, 
we should choose an equation of as simple a form as possible. Before 
proceeding any further with this choice we may test the correctness of 
the form of the equation by ‘‘rectifying’’ the curve, 7.e., by writing the 
assumed equation in the form 


fo) =a ON), or (2) y= ae be, 


where y’ = f(y) and x’ = F(x), and plotting the points with x’ and y’ 
as coordinates; if the points of this plot appear to lie on or very near a 
straight line, then this line can be represented by equation (2) and hence 
the original curve by equation (1). We shall use the method of rectifica- 
tion quite freely in the work which follows. 

Having chosen a simple form for the approximate equation we now 
proceed to determine the approximate values of the constants or co- 
efficients appearing in the equation. The method of approximation 
employed in determining these constants depends upon the desired degree 
of accuracy. We may employ one of three methods: the method of 
selected points, the method of averages, or the method of Least Squares. Of 
these, the first is the simplest and the approximation is close enough 
for a large number of problems arising in technical work; the second re- 
quires a little more computation but usually gives closer approximations; 
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while the third gives the best approximate values of the constants but the 
work of determining these values is quite laborious. All three methods 
will be illustrated in some of the problems which follow. 

After the constants have been determined the formula should be 
tested by performing several additional experiments where the variables 
lie within the range of the previous data, and comparing these results 
with those given by the empirical formula. 

We shall now work two illustrative examples to indicate the general 
method of procedure. 


(I) THE STRAIGHT LINE. 


69. The straight line, y = bx.— The following table gives the results 
of a series of experiments on the determination of the elongation £ in 
inches of annealed high carbon steel wire of diameter 0.0693 in. and gage 
length 30 in. due to the load W in pounds. 


W E EW w E} aes Eg Ar an Aut 
Olli) ° ONO ° ° ° ° ° 
50 | 0.0130 0.650 2,500 | 0.0130 | ©.013I | 0.0131 O}-— 1) —=s 

I0O | 0.0251 2.510 10,000 | 0.0260 | 0.0261 | 0.0262) — 9 | —10 | —11 

I50 | 0.0387 5.805 22,500 ||| 010300) |) 2 O).0302, 190.0393] "— 35s 5 a) 

200 | 0.0520] 10.400 | 40,000 | 0.0520 | 0.0522 | 0.0524 o|/+2!1—4 

225 | 0.0589 TBD 50,625 | 0.0585 0.0587 | 0.0589] + 4 | + 2 fe) 

250 | 0.0659 16.475 62,500 | 0.0650 | 0.0653 | 0.0655 + 9| +6] +4 

260 | 0.0689 17.914 67,600 | 0.0676 0.0679 | 0.0681] +13 | +10] + 8 

2 1235 | 0.3225 | 67.007 | 255,725 38 36 34 
2+8= 4.8 4.5 4.3 
ZA? = 356 270 25 


The plot. — The data are plotted on a sheet of codrdinate paper about 
10 inches square and ruled in twentieths of an inch or in millimeters. If 
we wish to express the elongation as a function of the load, we plot the 
load on the horizontal axis or as abscissas, if the load as a function of the 
elongation we plot the latter as abscissas. In Fig. 69 we have plotted 
the values of W as abscissas and the values of E as ordinates. The scales 
with which these values are plotted are generally chosen so that the 
length of the axis represents the total range of the corresponding vari- 
able, and so that the line or curve is about equally inclined to the two 
axes. There is no advantage in choosing the scale units on the two axes 
equal. Care should be taken not to choose the units either too small or 
too large; for in the former case the precision of the data will not be 
utilized, and in the latter case the deviations from a representative line 
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or curve are likely to be magnified. The drawing of a good plot is evi- 
dently a matter of judgment. It is best to mark the plotted points as 
the intersection of two short straight lines, one horizontal and one 
vertical. 


p 
Ss 
Dd 
LoS] 


Se aie 
Besse 


100 wy 200 250 


FIG. ih 


The representative curve and its equation. — We now draw a smooth 
curve passing very near to the points of the plot, so that the deviations 
of the points from the curve are very small, some positive and some 
negative. In Fig. 69, the points seem to fall approximately on a straight 
line. This should be tested by moving a stretched thread or by sliding 
a sheet of celluloid with a fine line scratched on its under side among the 
points and noting that the points do not deviate systematically from this 
thread or line. Having decided that a straight line will approximate 
the plot, we assume that an equation of the first degree, E = a + dW, will 
approximately represent the relation between the measured quantities. 
In this example we may evidently assume that EH = dW since a zero load 
gives a zero elongation. 
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The determination of the constant. — We shall now determine the con- 
stant bin the equation E=bW. This may be done in several ways. The 
three methods which are generally employed are as follows: 

I. Method of selected points. — Place the sheet of celluloid on the 
codrdinate paper so that the scratched line passes through the point 
W = 0, E = 0, and then rotate the sheet until a good average position 
among the plotted points is obtained, 7.e., until the largest possible num- 
ber of points lie either on the line or alternately on opposite sides of the 
line, in such a manner that the points below the line deviate from it by 
approximately the same amount as the points above it. Then note the 
values of W and £& corresponding to one other point on this line, prefer- 
ably near the farther end of the line. Thus we read W = 250, E = 
0.0650. Substituting these values in the equation E = bW, we have 
0.0650 = 250), and hence b = 0.000260, and finally E = 0.000260 W. 
Since the choice of the ‘“‘best’’ line is a matter of judgment, its position, 
and hence the value of the constant, will vary with different workers and 
often with the same worker at different times. 

II. Method of averages. — The vertical distances of the plotted points 


from the representative line are called the residuals; these are the differ- 
ences between the observed values of E and the values of E calculated 
from the formula, or E — E., where E. = DW; some of these residuals 
are positive and others are negative. If we assume that the “‘best’’ line 
is that which makes the algebraic sum of the residuals equal to zero, we 
have 
2 (E-—6W).=o or LE—bIW=o0, 

BE 440-3225 
SW 1235 
and we may call this an average value of 6. By this method it is no 
longer necessary to shift the line among the points so as to get an average 
position. 

III. Method of Least Squares. — In the theory of Least Squares * it is 
shown that the best line or the best value of the constant is that which 
makes the sum of the squares of the differences of the observed and cal- 
culated values a minimum, 7.e., 


x (EF — bW)? = minimum. 


Hence the derivative of this expression with respect to must equal zero, 
or 


hence : i 


= 0.000261, 


d 
5,2 (E-bW)?=0, or BW(E-bW) = 
LEW 
or LWE — ba W?2 = O; and |b = >We 


* See Bartlett’s “‘ The Method of Least Squares,” or any other book on this theory. 
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We form two columns, one giving the values of EW and the other the 
values of W?, and adding these columns, we find 


b = 67.007/255,725 = 0.000262. 


We may now compare the results obtained by each of the three 
methods. For this purpose we complete the table by computing the 
values of E from the formulas 


I. £ = 0.000260 W; II. E = 0.000261 W;° III. E = 0.000262 W. 


These are marked E, EY, E2™, in the table. To discover how closely 
the computed values agree with the observed values we form the residuals 


At =f Be Ed, Au = K-— ES, Au = K— EX, 


Disregarding the signs of these residuals, we add them and divide by 
their number, 8, and find the average residual to be 0.00048, 0.00045, 
0.00043, respectively. We also find the sum of the squares of the resid- 
uals to be 356, 270, 254, respectively. We may therefore draw the fol- 
lowing conclusions: all three methods give good results; the method of 
Least Squares gives the best value of the constant but requires the most 
calculation; the method of averages gives, in general, the next best value 
of the constant and requires but little calculation; the graphical method 
of selected points requires the least calculation but depends upon the 
accuracy of the plot and the fitting of the representative line. 

70. The straight line, y = a + bx. — For measuring the temperature 
coefficient of a copper rod of diameter 0.3667 in. and length 30.55 in., the 
following measurements were made. Here, C is the temperature Centi- 
grade and r is the resistance of the rod in microhms. 


& r Cc rc Te rit r ,it At Au Aut 
19.1 76.30 364.81 | 1,457.33] 76.190] 76.19] 76.26 |+o0.11\/-++0.11] +0.04 
25.0 77.80 625 .00-| 1,945.00|-77.91| 77.92) 77.96 |—o.11|—0.12| —o.16 
20ur 79.75 906.01 | 2,400.48] 79.39] 79.41] 79.43 {+0.36|/+0.34| +0.32 
36.0 80.80 | 1296.00 | 2,908.80) 81.11] 81.14) 81.13 |—0.31/—0.34| —0.33 
40.0 82.35 | 1600.00 | 3,294.00] 82.27] 82.31; 82.28 |+0.08]/+0.04| +0.07 
45.1 83.90 | 2034.01 | 3,783.89] 83.75} 83.80} 83.76 |+o0.15|/+0.10} +0.14 
50.0 85.10 | 2500.00 | 4,255.00] 85.18) 85.24) 85.16 |—o.08/—o.14| —o.06 
Z 245.3 | 566.00 | 9325.83 | 20,044.50 120} ere D9 Ef gts 
U7 = 0.171 | 0.1701)" 0.160 

’ DA*2= 2852 | 2869 2646 


The plot (Fig. 70) appears to approximate a straight line, so that we 
shall assume the relation r= a-+0C. We shall determine the con- 
stants, a and 0, by the three methods. 

I. Method of selected points. — Use a sheet of celluloid to determine 
the approximate position of the best straight line, and note two points 
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on this line; thus, C = 20, 7 = 76.45, and C = 48, r = 84.60. Substi- 
~~ tuting these values in the equation 7 = a + OC, we get 
76.45 =a+20b and 84.60 =a-+ 48), 
from which we determine 
"2 a= 70.63 and b = 0291, 
so that our relation becomes 
a r = 70.63 + 0.291 C. 


~\ 


MS 10 75 20 25 #30 85 40 50 
, (0) 
Fic. 70. 


\— IIL. Method of averages. — Since we have to determine two constants, 
\._ we divide the data into two equal or nearly equal groups, and place the 
sum of the residuals in each group equal to zero, 7.e., 


i x (r —a— 6C) = 0 or ar = na + bxC, 


_ where 7 is the number of observations in the group. Thus, dividing the 
“above data into two groups, the first containing four and the second three 
sets of data, and adding, we get 


314.65 =4@a@+ 11026 and 251.35 = 3a + 135.1 8, 
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from which we determine 
a= 70.59 and 6b = 0.203, 
so that our relation becomes A. 
r = 70.59 + 0.293 C. 
III. Method of Least Squares. — The best values of the constants are 
those for which the sum of the squares of the residuals is a minimum, 7.e., 


2 (r — a — dC)? = minimum; hence the partial derivatives of this 
expression with respect to a and 6 must be zero; thus, 


£3 (r-a-sC) =0, SE (r—a- 00} =0, 


or Xl[2(r—a-—dbC)(-1)])=0, Yl2(r—a-—dC)(—-C)j =0, 
or Lr = an + bXxC, 
BrC = alc + bC, 


where 7 is the number of observations. We solve these last two equations 
for a and b. (Note that these equations may be formed as follows: 
substitute the observed values of r and C in the assumed relation r = a 
+ bC; add the m equations thus formed to get the first of the above 
equations; multiply each of the m equations by the corresponding value 
of C and add the resulting equations to get the second of the above 
equations.) 

We now compute the values of rC, C?, UC, XrC, and XC?, and substi- 
tute these in the equations for determining a and b. We thus get 


566.00 = 7a + 245.3), 
20,044.50 = 245.3. @ + 9325.83 ), 
from which we determine 
a= 70.76 and b = 0.288, 
so that our relation becomes 
r = 70.76 + 0.288 C. 


Comparison of results. — We note that the various results agree very 
well with the original data and with each other. We compute the resid- 
uals and find that the average residual is smallest by the third method 
and is approximately the same by the first two methods. The computa- 
tion necessary in applying the method of Least Squares is very tedious. 
The method of selected points requires the fitting of the best straight line, 
and this becomes quite difficult when the number of plotted points is 
large. We shall therefore use the method of averages in most of the 
illustrative examples which follow. 
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(II) FORMULAS INVOLVING TWO CONSTANTS. 


71. Simple parabolic and hyperbolic curves, y = ax’. — As stated in- 
Art. 68, when the plotted points deviate systematically from a straight 
line, a smooth curve is drawn so as to pass very near the points; the shape 
of the curve or a knowledge of the nature of the experiment. may give us 
a hint as to the form of the equation which will best represent the data. 

Simple curves which approximate a large number of empirical data 
are the parabolic and hyperbolic curves. The equation of such a curve 
is y = ax, parabolic for 6 positive and hyperbolic for 6 negative. In 
Fig. 71a, we have drawn some of these curves for a = 2 and b = —2, 
—I, —0.5, 0.25, 0.5, 1.5,2. Note that the parabolic curves all pass 


Fic. 71a. 


through the points (0, 0) and (1, a) and that as one of the variables 
increases the other increases also. The hyperbolic curves all pass through 
the point (1, @) and have the codrdinate axes as asymptotes, and as one 
of the variables increases the other decreases. 

There is a very simple method of verifying whether a set of data can 
be approximated by an equation of the form y = ax*. Taking loga- 
rithms of both members of this equation, we get log y = loga + b log x, 
and if x’ = log x, y’ = log y, this becomes y’ = log a + 6x’, an equation 
of the first degree in x’ and y’; therefore the plot of («’, y’) or of (log x, 
log y) must approximate a straight line. Hence, 
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If a set of data can be approximately represented by an equation of the 


se mean an esha 


1G y = ax’, then the plot of (log x, log y) approximates a straight line. 

Instead of plotting (log x, log y) on ordinary coordinate paper, we 
may plot (x, y) directly on logarithmic codrdinate paper (see Art. 13). 
We determine the constants a and 0b from the equation of the straight 
line by one of the methods described in Art. 70. 

Example. The following table gives the number of grams S of anhy- 
drous ammonium chloride which dissolved in 100 grams of water makes 
a saturated solution of 6° absolute temperature. 


log 6) 
2.43 2.44 245 246 2.47 2.48 2.49 pete 252 253 254 255 256 257 


\A 
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6 iS log 0 log S St - Sot At At 

273 29.4 2.4362 1.4684 29.7 20.7 —0.3 —0.3 
283 33-3 2.4518 1.5224 Bana Bene ONL +o.1 
288 B5e2 2.4594 1.5465 35.0 BGolt =-0).2 +0.1 
293 B72 2.4669 1.5705 37.0 37.0 --O.2 +o.2 
313 45.8 2.4955 1.6609 Wixesy ape Oa qr Oa +0.5 
333 55-2 2.5224 1.7419 54.9 54.9 +0.3 +0.3 
353 65.6 2.5478 1.8169 65.7 65.8 —o.1 —0.2 
373 Whee 2.5717 1.8882 tipo) 78.0 —o.6 —0.7 
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The points (6, S) are plotted in Fig. 71b. The curve appears to be 
parabolic, i.e., of the general form illustrated in Fig. 71a. We therefore 
plot (log @, log S)-and note that this approximates a straight line, so that 
we may assume 

S=aée or log S = loga+ blog ?é. 
We shall first determine the constants by the method of selected 
points. We note two points on the line whose codrdinates are 
log 0 = 2.445, log S= 1.50 and logé@ = 2.555, log S = 1.84, 
hence we have 
1.50 = loga + 2.445 J, 
1.84 = log a + 2.555 0d. 
b = 3.09, loga = —6.0550 = 3.9450 — 10, @ = 0.000,000,881. 
log S = —6.0550 + 3.09 log #, or S = 0.000,000,881 6%. 


We shall now determine the constants by the method of averages. 
We divided the data into two groups of four sets, and adding, we have 


6.1078 = 4loga+ 9.8143 dD, 

7.1079 = 4 loga + 10.1374 Db. 
*. b = 3.09, log a = —6.0546 = 3.9454 — 10, @ = 0.000000882. 
“. log S = —6.0546 + 3.09 log@ or S = 0.000000882 43%, 


We complete the table by computing S, the residuals, and the average 
residual. The agreement between the observed and computed values of 
5 is quite close. 

Example. The following table gives the pressure p in pounds per 
sq. in. of saturated steam corresponding to the volume v in cu. ft. per 
pound. (From Perry’s Elementary Practical Mathematics.) 


v Dp log v log pD Ds A 
53-92. 6.86 1.7318 0.8363 6.85 aIpOLOL 
26.36 I4.70 I.4210 I .1673 14.69 +o.o1 
14.00 28 .83 1.1461 1.4599 28 .85 —©).02 

6.992 60.40 0.8446 1.7810 60.49 —0.09 
4.280 IOI .9 0.6314 2.0082 102.1 —On 
2.748 163.3 0.4390 2.2130 163.7 —0.4 
1.853 250.3 0.2679 2.3984 249.2 SPs 


The points (v, p) are plotted in Fig. 71c. The curve appears to be 
hyperbolic on comparison with Fig. 71a. Hence we plot (log v, log p) 
and note that this approximates a straight line, so that we may assume 


p=av, or log p = loga+ dlogy. 


We shall use the method of averages to determine the constants a and b. 
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Dividing the data into two groups, the first four and the last three sets, 
and adding, we have 
5:2445 = 4 log a + 5.1435 8, 
6.6196 = 3 log a + 1.3383 bd. 
mn 0N= — 10662, logy a =!2.0822, @ = 481: 1. 
log p = 2.6822 — 1.0662 logyv, or pv! = 481.1. 


We now compute f and A and note the close agreement between the 
observed and calculated values. 


(log v) 
0.2 0.4 0.6 0.8 1.0 1.2 1.4 16 1.8 

260 
240 
220 es 2.4 
200 i 2 
180 L 2.0 
160 k& 1.8 

Sle y 
140 os 1.63 
(P) e) — 
120 1.4 
100 1,2 

rm 

80 : 1.0 
60 eS 0.8 
40 
20 7 22) 


0 
GO. 5240 1820 725 Cy 85 40 45 50 55 
v 
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72. Simple exponential curves, y = ae’”.— Other simple curves that 
approximate a large number of experimental results are the exponential 
or logarithmic curves. The equation of such a curve may be written in 
the form y = ae’*, where e is the base of natural logarithms; the form 
y = ab* is sometimes used. In Fig. 72a, we have drawn some of these 
curves for a = 1 and 6} =—2, —I, —0.5, 0.5, 1, 2. Note that these 
curves all pass through the point (0, a) and have the x-axis for asymptote. 
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There is a very simple method of verifying whether a set of data can 


be approximated by an equation of the form y = ae??. 


Taking logarithms 
of both members of this equa- 
tion we get logy = loga + 
(b log e) x, and if y’ = logy, 
this equation becomes y’ = 
log a + (blog e) x, an equation 
of the first degree in x and y’; 
therefore the plot of (x, y’) or 
of (x, log y) must approximate 
a straight line. Hence, 

If a set of data can be ap- 
proximately represented by an 
equation of the form y = ae, 
then the plot of (x, log y) ap- 
proximates a straight line. 

Instead of plotting (x, logy) 
on ordinary coérdinate paper, 
we may plot (x, y) directly 
on semilogarithmic codrdinate 
paper (see Art.14). Thecon- 
stants a and 6 are determined 
from the equation of the 


straight line by one of the methods described in Art. 70. 

Example. Chemical experiments by Harcourt and Esson gave the 
results of the following table, where A is the amount of a substance re- 
‘maining in a reacting system after an interval of time ¢. 


t A log t log A ae A 

2 94.8 0.3010 1.9768 94.9 —-O.0 
5 87.9 0.6990 I.9440 87.7 =-0.2 
8 81.3 ©.9031 I.QIOI 81.0 “0.3 
TEL 74.9 I .O414 1.8745 74.8 OnE 
14 68.7 I.1461 1.8370 69.1 —=0.4 
17 64.0 I .2304 1.8062 63.8 +o.2 
27 49.3 1.4314 1.6928 49.0 0.3 
31 44.0 I .4914 1.6435 44.1 Ont 
35 39.1 1.5441 1.5922 39.6 —0.5 
44 31.6 1.6435 1.4997 31.2 +0.4 


The points (¢, A) are plotted in Fig. 720. 


ZA +10 = 0.26 
This curve appears to be 


exponential, so that we plot (¢, log A) and (log t, A); it is seen that the 


plot of (¢, log A) approximates a straight line. 


an equation of the form 


A = ae 


We may therefore assume 


or log A = loga + (b loge) t. 
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(log t) 
OR OE US ONS BO A CS Bs 
—- 
95 
90 2.00 
85 1,95 
\ 
80 ie 1.90 
75 = 1.85 
I 
70 %, | | 1.80 
+' 
lel K 
65 | cS 1.75 
{A) | |_| _|(log A) 
60 1.70 
1 (Be a 
55 pals Zz 11.06 
-50 ; 1.60 
T 
45 7 ei Paths 1.55 
al 
40 7.50 
Be ; iN 
35 1.45 
ae 
30 
0 5 10 15 20 ai, A) 30 35 40 45 
Fic. 720. 
We shall use the method of averages to determine the constants. Divid- 


ing the data into 2 groups and adding, we get 


9.5424 = 5 loga + 40 (0 loge), 
8.2344 = 5 loga + 154 (0 loge). 


lth 


log A = 2.0005 — 0.0115 f, 


— 0.0265, 


blog e = —0.01I5, 


log @ = 2.0005. 
@ = 100.1, since log ¢ = 0.4343. 
or A = 100.1 ¢7?.%§?, 


We now compute the values of A and the residuals, and note the close 
agreement between the observed and the calculated values of A. 
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Example. The following table gives the results of measuring the 
electrical conductivity C of glass at temperature 6° Fahrenheit. 


Bienen 
PEELE 


80’ 100 120 140 160 180 200-7220 240° 
| (6) 


FIG. 72c. 

0 (G; log 6 log C (on A 

58 ° 1.7634 —o ©.0019 

86 ©.004 1.9345 7 .6021—10 © .0039 +0 .0001 
148 0.018 DEL 703 8.2553—10 0.0185 —0.0005 
166 0.029 2.2201 8 .4624—10 0.0292 —0.0002 
188 0.051 2.2742 8.7076—10 0.0510 ° 

202 0.073 2.3054 8.8633 —I0 0.0728 +0.0002 
210 0.090 2.3222 8.9542—10 0.0891 —90.0010 
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In Fig. 72c, the potnts (6, C) and (6, log C) are plotted; the latter 
plot approximates a straight line. We may therefore assume the equa- 
tion 

C = ae, or log C = loga+ (6b loge) 0. 


We use the method of averages to determine the constants. Omitting 
the first set and dividing the remaining data into two groups of three 
sets, we get 
24.3198 — 30 = 3 loga + 4oo (0 loge), 
26.5251 — 30 = 3 loga + 600 (0 loge). 
bloge = 0.0110, loga = 6.6399 — Io. 
b = 0.0253, @ = 0.000436. 
log C = 6.6399 — 10 + 0.01108; or C = 0.00436 e€°-02588, 


We now compute the values of C and the residuals and note the re- 
markably close agreement between the observed and computed values 
of C. 

73. Parabolic or hyperbolic curve, y = a + bx” (where nis known). — 
In using this equation, it is assumed that from theoretical considerations 
we suspect the value of n. It is evident that 

If a set of data can be approximately represented by an equation of the 
form y = a+ bx", where n is known, then the plot of (x", y) approximates 
a straight line. 

Example. A small condensing triple expansion steam engine tested 
under seven steady loads, each lasting three hours, gave the following 
results; J is the indicated horse-power, w is the number of pounds of 
steam used per hour per indicated horse-power. (From Perry’s Ele- 
mentary Practical Mathematics.) 


I w wl DW, A 
36.8 1255 460.0 12.6 Oni 
ies 12.9 406.4 Te fe +o0.1 
26.3 ite) at 344.5 T3350 -+o.1 
21.0 13.3 279.3 wge4 —O.1 
15.8 Ae 222.8 I4.0 —-o.r 
12.6 14.5 182.7 14.6 —o.I 
8.4 TOES 136.9 16.1 +o.2 


YA +7 =0.11 


Fig. 73a gives the plot of (J, w). This is not a straight line. But if 
we plot (J, wl), 7.e., the total weight of steam used per hour instead of 
the weight per indicated horse-power, we find that this plot approximates 
a straight line. Hence, we may assume the linear relation wl = a + OI. 
This relation may also be written w = 6 +<a/T, so that the plot of (1/J, w) 
also approximates a straight line. We use the method of averages to 
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letermine the constants. Dividing the data into two groups, the first 
three and last four sets, and adding, we have 
1210.9 = 3a + 94.6), 
821.7 = 4a+ 57.8). 
b= 11.6, a= 37.8. 


8 
wl = 37.8 + 11.6 If. Or w= 116 +35 x 


We now compute the values of w and the residuals. 


FIG. 73a. 


Example. For a parachute or flat plate falling in air we have the 
following observations; v is the velocity in ft. per sec. and is the pressure 
in pounds per sq. in. 


v D v De A 
7.87 0.2 61.94 0.187 —OnOrS 
II.50 0.4 ROS JOS ©.401 +o .001 
16.40 0.8 268 .96 0.815 —0.015 
22.60 1.6 510.76 1.548 +0.052 
32.80 Sa 1075.84 3.260 —0o.060 
en ene ee ee ee ree Bete ee ee 


ZA + 5 = 0.028 


, x 
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In Fig. 730, we have plotted (v, p). It is surmised that for low veloci- 

ties, the pressure and the square of the velocity are linearly related, i.e., 

p=a-+ bv. We verify this by plotting (v, p) and noting that this 

approximates a straight line. We use the method of averages to deter- 
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mine the constants. Dividing the data into two groups, the first three 
and the last two sets, and adding, we have 


1.4=3a+ 463.15), 
4.8 = 2a + 1586.60 b. 


b = 0.00303 and @ = —O.OOIII. 
p = —O.O0III + 0.00303 v. 


We may with good approximation take a = 0, so that p = 0.00303 v”, i.e., 
the pressure varies directly as the ae of BS velocity. 


74. Hyperbolic curve, y = ae Frege es = a-+ bx. —This equation 
represents the ordinary hyperbola with Sana x = —a/b and y = 
1/6; as illustrated in Fig. 74a for values of a = 0.2, b = 0.2; a = O.I, 
b = 0.2; a = —0.1, b =0.2; a = —0.2, b = 0:2. Quite a large number 


of experimental results may be represented by an equation of this type. 
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; : I a 
The equation may also be written in the forms 9 ==") sent so that the 


plots (x, *) and @ : approximate straight lines. Hence, 
y . 
If a set af data can be approximately represented by an equation of the 
I : 
form y = oe Ah a + bx then the plot of (x, 4 or of (2, 4 approxt- 


mates a straight jae 


FIG. 74a y= 


x 
a+ bx 
Example. From a magnetization or normal induction curve for 
iron we find the following data; H is the number of Gilberts per cm., a 


measure of the field intensity, and B is the number of kilolines per sq. cm., 
a measure of the flux density. 


H B H/B Be A 
2.5 3-5 ©.714 7-97 

3.0 5.0 ©.600 8.78 

Bow 7.5 0.413 8.91 

3.8 I0.0 0.380 9.8 +o. 
TAO T25 ©.560 12.4 +o.1 
9.5 ToS 0.703 13.6 —OFd 
ge 14.0 © .808 14.0 ° 
WH eS 15.0 cai Toee —o.1 
31.5 16.0 1.97 TOR —On2 
45.0 16.5 BLD TORY One 
64.0 T7O) 3.76 L7LO fo) 
95-2 17-5 5-43 17-3 +0.2 


ae eS ee ee 


ZA 0) — Oi 


In Fig. 746, (H, B) is plotted. The curve appears to be of the type 
illustrated in Fig. 74a. Furthermore, an important quantity in the 
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theory of magnetization is the reluctance, H/B, and if we plot (H, H/B), 
we note that this plot approximates a straight line for values of H > 3.1. 
(We may similarly introduce the permeability, B/H, and note that the 
plot of (B/H, B) approximates a straight line.) Hence, we assume a 


relation of the form = a+bH. Using the method of averages, 


B 
| a 

; L 5.5 
rai sisal Sal anlieat 5.0 

18 | 4.5 
16 114.9 

>) 
14 oe i 3.5 
ri <a ae 
(B) !? | ; oh 3.0(#) 
| 8 
10 2.5 
poe | 
8 | a T 2.0 
| Bee 
6 1 | 1.5 
| he] 
4 T 7.0 
a 
2 0.5 
09-40 20.30 4050 60.70 80 90 100° 
(7) 
Fic. 740. 


omitting the first three values of H, and dividing the remaining data 
into two groups containing five and four sets respectively, we get the 
equations 
3.621 = 5a + 49.1 B, 
13.88 = 4a + 235.5 0. 
= 0.0560, @ = 0.174. | 
H 
0.174 + 0.0560 H 
We now compute B and the residuals and note the close agreement 
between the observed and computed values. 


oe a= 0.174 + 0.0560H or B= 
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(III) FORMULAS INVOLVING THREE CONSTANTS. 


75. The parabolic or hyperbolic curve, y = ax’ + c.—It is often im- 
possible to fit a simple equation involving only two constants to a set of 
data. In such cases we may modify our simple equations by the addition 
of a term involving a third constant. Thus the equation y = ax’ may be 
modified into y = ax’+ c. If 6 is positive, the latter equation repre- 
sents a parabolic curve with 
intercept c on OY; if bd is 
negative, the equation rep- 
resents a hyperbolic curve 
with asymptote y=c. In 
Fig. 75a, we have sketched the 
curves y=2 x95, y=2 °*+ 2, 
Ay as 2 On 2 ee ee 
illustrate the relation of the 
simple types to the modified 
types. 

In Art. 71 it was shown 
that if we suspect a relation 
of the form y = ax’, we can 

0 Oo. Snake 10 verify this by observing 

©) whether the plot of (log x, 

og ee log y) approximates a straight 

EIS 152. line. Now the form y = ax® 

+c may be written log (y — c) = loga + db log x, so that the plot of 

(log x, log (y — c)) would approximate a straight line. To make this 

test we shall evidently first have to determine a value of c. We might 

attempt to read the value of c from the original plot of (x, y). In the 

parabolic case we should have to read the intercept of the curve on OY, 

but this may necessitate the extension of the curve beyond the points 

plotted from the given data, a procedure which is not safe in most cases. 

In the hyperbolic case, we should have to estimate the position of the 
asymptote, but this is generally a difficult matter. 

The following procedure will lead to the determination of an approxi- 
mate value of c for the equation y = ax + c. Choose two points (x, y:) 
and (x, ye) on the curve sketched to represent the data. Choose a third 
point (xs, y3) on this curve such that x3 = Vx,x2, and measure the value 
of y3. Then, since the three points are on the curve, their codrdinates 
must satisfy the equation of the curve, so that 


MW = ax +c, Yo = axe + 6, 3 = axs? +c. 
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Now, since x3 = V x1%., 


therefore x? = V xx, and axe = VaxP- aX2, 
or DE Se Vy — 6) Goi Os 
= 2 
and therefore ee ee eee 
Vi ys = 23 


It is evident that the determination of c is partly graphical, for it 
depends upon the reading of the codrdinates of three points on the curve 
sketched to represent the data. The curve should be drawn as a smooth 
line lying evenly among the points, 7.e., so that the largest number of the 
plotted points lie on the curve or are distributed alternately on opposite 
sides and very near it. 

Having determined a value for c, we plot (log x, log (y — c)). If this 
plot approximates a straight line, the constants a and 6 in the equation 
log (y — c) = log a+ b log x may then be determined in the ordinary 
way. 

Example. Ina magnetite arc, at constant arc length, the voltage V 
consumed by the arc is observed for values of the current 7. (From 
Steinmetz, Engineering Mathematics.) 


1 V V — 30.4 log (V — 30.4) logi Ve A 
Ons 160 129.6 220 9.6990 — I0 158.8 = eee? 
I 120 89.6 Tob 28 ©.0000 -- 10 120.8 Ono 
2 94 63.6 1.8035 0.3010 — Io 94.0 ° 
4 75 44.6 I .6493 0.00213 — TO Tis —o.1 
8 62 31.6 I.4997 On003; aro 61.9 -+o.1 

12 56 25.6 1.4082 1.07.92)— 10) 56.0 fo) 


We plot (z, V) and note that the curve appears hyperbolic with an 
asymptote V = c, and hence we assume an equation of the form V= az? + ¢. 
To verify this we must first determine a value for c. Choose two points 
on the experimental curve; in Fig. 75), we read 7 = 0.5, Vi = 160 and 
ig = 12, V2 = 56. Choose a third point such that i = Vive = V6 = 
2.45, and measure V3 = 88. Then 


pees ee (100) (56) (88) 1276 
VitV2—2V3; 160+56—2(88) 40 30-4. 


Now compute the values of V — 30.4 and log (V — 30.4) and plot 
(log 7, log (V — 30.4)). This last plot approximates a straight line so 
that the choice of the equation V = az’ + c is verified. 

To determine the constants in the equation 


log (V — 30.4) = loga + b log z, 
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we use the method of averages, dividing the data into two groups of three 
sets each, and find 


5.8684 = 3 loga, 
4.5572 = 3 loga + 2.5844 b. 
b = —0.507, loga = 1.9561, @ = 90.4. 
log (Vv. — 30.4) = 1.9561 — 0.507 logz, or V = 30.4 + 90.47°°™, 


Finally, we compute the values of V and the residuals. 


log 7) 
ene 9.7 9.8 9.9 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 Oh 


76. The exponential curve, y = ae” + c. — The simple exponential 
equation y = ae* may have to be modified into y = ae* + c in order to. 
fit a given set of data. In the latter curve, the asymptote is y =. 
In Fig. 76a, we have sketched the curves y=2 17, y=2 @l7+1, 
y= 2e0l2 y= Qe%le4 1, 

In Art. 72 it was shown that if we suspect a relation of the form 
y = ae", we can verify this by observing whether the plot of (x, log y) 
approximates a straight line. Now y = ae*=+c may be written 
log (y — c) = loga + (0 log e) x, so that the plot of (x, log (y — o)) 
would approximate a straight line. Evidently we shall first have to 
determine a value for c. We proceed to do this in a manner similar to 
that employed in Art. 75. Choose two points (x, 41) and (x2, ye) on 

’ 
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the curve sketched to represent the data, and then a third point (x3, y3) on 
this curve such that x3 = $ (x1 -+ x2) and measure the value of 3. Since 
the three points are on the curve.. 

Bigeye, Vor sder ic) oS ye = ae FC, 


or log a = (bloge) x1, log a es (b log e) x2, log wees = (blog e) x3. 


Fic. 76a. y = ae*+c 


Now, since x3 = 4 (x1 + x), 
therefore (b log e) x3 = 4 [(6 log e) x1 -+ (6 log e) x], 
of ee eee Ne Seiace ee (| a ye 
and log eva [toe Are + log ae | log ee Lee 
— 2 
Hence y3 —c = V(yi—c) (2 —c), and c= BREADS iat NEE 


Vie Ye — 23 


If the data are given so that the values of x are equidistant, 7.e., so 
that they form an arithmetic progression, we may verify the choice of 
the equation y = ae’* + c and determine the constants a, 6, and c in the 
following manner. Let the constant difference in the values of x equal h. 
If we replace x by x + h, we get y’ = ae") +c, and therefore, for the 
difference in the values of y, 


Ay = y' — y = ae) — ae = ae” (e — 1), 
and log Ay = log a (e — 1) + (b loge) x. 
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This last equation is of the first degree in x and log Ay so that the plot of 
(x, log Ay) is a straight line. To apply this to our data, we form a 
column of successive differences, Ay, of the values of y, and a column of 
the logarithms of these differences, log Ay, and plot (x, log Ay); if the 
equation y = ae’ +c approximates the data, then this last plot will 
approximate a straight line. We may then determine }loge and 
log a (& — 1) and hence a and 6 in the ordinary way, and finally find an 
average value of c from Ly = ade + nc, where n is the number of data. 

Example. In studying the skin effect in a No. 0000 solid copper 
conductor of diameter 1.168 cm., Kennelly, Laws, and Pierce found the 
following experimental results; F is the frequency in cycles per second, 
L is the total abhenrys observed. 


F 1b; L — 51,860 log (L — 51,860) Jb, A 
60 53,912 2052 Bes iz2 53,952 —40 
306 53,767 1907 3.2804 53,668 +99 
888 53,143 1283 3.1082 53,140 + 3 
1600 52,669 809 2.9079 52,699 —30 
2040 52,499 639 2.8055 52,506 a 
3065 SOG 355 2.5502 Keke + 3 
3950 2,082 222 2.3464 52,068 +14 
5000 51,965 105 2.0212 51,972 —7 


In Fig. 760, the points (F, L) are plotted; the curve appears to be 
exponential with an asymptote L = c. We shall try to fit the equation 
L=ae" +c. First determine an approximate value for c by choosing 
two points on the experimental curve, Fi; = 875, Li = 53,140, and F) = 
5000, Lz = 51,980, and a third point F; = 4 (F, + Fy) = 2938, L3 = 

LL, — L;? 
Ii + L,—21; 
and log (L — 51,860), and plot (F, log (L — 51,860)); this plot approxi- 
mates a straight line, thus verifying the choice of equation. We deter- 
mine the constants in the equation log (L — 51,860) = log a + (b log e) F 
by the method of averages. Dividing the data into two groups of four 
sets each and adding, we have 


52,250.- Then ¢ = = 51,860. Now compute (L — 51,860) 


12.6087 = 4 loga + 2854 b loge, 
9.7233 = 4 loga + 14,055 b loge. 


blog e = —0.0002576, log a = 3.3360, 
and b = —0.0005931, @ = 2168. 
log (L — 51,860) = 3.3360 — 0.0002576 F, 
or L = 51,860 + 2168 e~0-0005931 F 
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We now compute L and the residuals, and note the close agreement 
between the observed and computed values except for the first two values 
of F. If we omit these two values in computing a and 3, these constants 
have slightly different values, but the agreement between the observed 
and computed values of L is about the same. 


Fic. 760. 


77. The parabola, y = a+ bx + cx?. — The equation of the straight 
line y = a + bx may be modified by the addition of a term of the second 
degree to the form y= a+ bx + cx?. This is the equation of the ordi- 
nary parabola. We may verify whether this equation fits a set of experi- 
mental data by one of the following methods. 
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(1) Choose any: point (x,,y,) on the experimental curve; then 
Ye = a + bxy + cx”, and 
y — yn = O(x — xx) +(x? — x7), oF ae = = (b+ cx.) + cx. 


This last equation is of the first degree in x and a= so that the plot of 
med) 


(x, : = - ‘) will approximate a straight line. 
Se VK 


(2) If the values of x are equidistant, 7.e., if they form an arithmetic 
progression, with common difference h, then if we replace x by x + / in 
the equation, we get y’ =a+b(xa+h)+c(x+h)? and Ay=y'—y 
= (bh + ch?) + 2chx. This last equation is of the first degree in x and 
Ay, so that the plot of (x, Ay) will approximate a straight line. 

Hence, if a set of data may be approximately represented by the equation 
y=a+ bx + cx’, then (1) the plot of (x, » = 2), where (Xz, Ye) are the 
coérdinates of any point on the experimental curve, will approximate a 
straight line, or (2) the plot of (x, Ay), where the Ay’s are the differences in 
y formed for equidistant values of x, will approximate a straight line. 

The following examples will illustrate the method of determining the 
constants. 

Example. In the following table, 6 is the melting point in degrees 
Centigrade of an alloy of lead and zinc containing x per cent of lead. 
(From Saxelby’s Practical Mathematics.) 


x 6 x — 36.9 6 — 181 Rae 6. A 
87.5 292 50.6 III 2.20 295 —3 
84.0 283 Apion 102 ay 05) 285 —2 
77.8 270 40.9 89 2.18 268 +2 
63.7 235 26.8 54 2.01 234 cfed 
46.7 197 9.8 16 1.63 199 —2 
36.9 181 ° ° 182 —I 


In Fig. 77a, we have plotted (x, @). We shall try to fit an equation 
of the form @ = a + bx + cx? to the data. To verify this choice, observe 
that the curve passes through the point x, = 36.9, 6, = 181, and plot the 
points (« eet | 

(f—1300 
plotting the ordinates for the straight line a scale unit ten times as large 
as that used for the ordinates of the experimental curve has been used: 
any further increase in the scale unit would simply magnify the devia- 


) ; this last plot approximates a straight line. (In 
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6 — 181 
<= 36:9 
the method of averages to determine the constants. Dividing the data 
into two groups of three and two sets respectively and adding, we get 


tions.) We may now assume the relation =a’ + b’x, and use 


6.55 = 3a’ + 249.3 0’, 

3.64 = 2a’ + 110.40’. 

6’ =.0.0130, @® = 1.10. 
g EL a aaa Ox, or 0= 141.4 + 0.620x + 0.0130 x? 
ST Were rae 01204; 0 = PAT.4 .620 x O140M- 


We now compute @ and the residuals. 


‘300 


y _ 7 436.9 ft TY 2.0 


ae LS) 
che A apo 1 co 


70 89 90 


30 40 50 60 
(x) 
Fic. 77a. 
Example. The following table gives the results of the measure- 


ments of train resistances; V is the velocity in miles per hour, R is the 
resistance in pounds per ton. (From Armstrong’s Electric Traction.) 
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V R AR v2 Ie A 
20 ais 3.6 400 5.70 OAS 
40 Ot 5.8 1,600 9 .08 +0 .02 
60 14.9 7.9 3,000 14.82 +0.08 
80 DD fe) 10.5 6,400 22.86 —o.06 
100 B3n3 127 10,000 33-22 +0.08 
120 46.0 14,400 45.90 +o.10 
Zz 420 I31I.6 36,400 


In Fig. 770, the plot of (V, R) appears to be a parabola, R = a + bV 
+ cV?. Since the values of V are equidistant, we shall verify our choice 
of equation by a plot of (V, AR); this last plot approximates a straight 
line. We may therefore assume AR = (bh + ch?) + 2chV, where h = 20. 


50 


45 af = 


40 


35 ae 


30 15 


(R) 25 10 


0 20 40 60 700 120 


60 
(V) 
Fic. 770. 


We determine the constants in this last equation by the method of 
averages, using the five sets of values of V and AR. Dividing these data 
into two groups of three and two sets respectively and adding, we get 


17.3 = 3 (bk + ch?) + 120 (2 ch), 

23.2 = 2 (bh + ch?) + 180 (2 ch). 

fon? Ch =O) Oech = 08, 
c¢ = 0.0029, 6 = —0.004. 

R =a — 0.004 V + 0.0029 V?. 
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We determine a by substituting the six sets of values of V and R, and 
summing, thus 
XR = 6a — 0.004 UV + 0.0029 LV, 
or 131.6. = 6a — 0.004 (420) + 0.0029 (36,400), 
and therefore @ ="4.62. 
Hence, finally, R = 4.62 — 0.004 V + 0.0029 V?. 


We now compute the values of R and the residuals; the agreement 
between the observed and porate values of R is very close. 


78. The hyperbola, y = + c.— This equation is a modifica- 


abs 


tion of the equation y = discussed in Art. 74. In the latter 


ea 


equation, x = 0 gives y = 0, while in the former, x = o gives y = c. We 


may verify whether the equation y = +c fits a set of experi- 


a 
a+ bx 
menta! data as Ss Choose any point (xz, y,) on the experimental 


+c, and 


curve; then y, = 


ere 


x OAL Xs) (De ere b 
ee ba) a on)" or Pe a a 


This last equation is of the first degree in x and 


of (x, 


Hence, if a set of data may be pcre represented by the equation 


i peta =, the plot of (x, = 


of a point on the experimental an will approximate a straight line. 
Example. The following table gives the results of experiments on 
the friction between a straw-fiber driver and an iron driven wheel under 
a pressure of 400 pounds; y is the coefficient of friction and x is the slip, 
per cent. (From Goss, Trans. Am. Soc. Mech. Eng., for 1907, p. 1099.) 


= aS so that the plot 
— Vk 


= =) will approximate a straight line. 
ae Ke, 


=), where (xXx, Ye) are the codrdinates 
k 


% y x — 0.65 y — 0.129 oa = Yo ye 
0.65 0.129 fo) ° 0.129 0.129 
0.87 0.217 0.22 0.088 2.50 0.253 0.228 
0.88 0.228 0.23 0.099 Boge) 0.256 0.232 
0.90 0.234 Qw2s 0.105 2.38 0.264 0.238 
0.93 0.275 0.28 ©.146 I.Q2 0.274 0.248 
1.16 0.318 O.51 0.189 270 0.326 0.304 
1.80 0.400 Tess 0.271 4.25 0.394 0.388 
By) ©.410 1.47 0.281 Ges 0.410 ©.411 
3.00 0.435 2335 © .306 7.68 | 0.435 0.451 
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In Fig. 78 we have plotted the points (x, y); the experimental curve 
: x 
appears to be an hyperbola with an equation of the form y = eae + ¢. 


To verify this we note the point x = 0.65, y = 0.129 on the curve, and 


F — 0.6 
plot the points (x, See eo! 


) This last plot approximates a straight 
y — 0.129 


— 0.65 


Uae 
line. We may therefore assume the relation =a-+ bx, and 


Vv — 0.129 


0.4 
0.3 
(y) 
0.2 
0.1 
0 
0 0.5 1 1.5 mare 2.5 3 
(x) 
Fic. 78. 


we shall determine the constants by the method of averages. As the 
first three points do not lie very near this straight line, we shall use only 
the last five sets of data, and dividing these into two groups of three and 
two sets respectively and adding, we get 
8.87 = 3a + 3.89), 
12.91 = 2a+ 5.12). 
b= 2.77, a= —0.64. 
Ki 0.05 = pee kes 
y— 0129 _ Raat Fh 2.77 x — 0.64 
If we had used all eight points in determining the constants, we should 
have obtained 


+ 0.129. 


9.12 = 4a + 3.58 8, 
19.86 = 4a + 8.08 bd. 
b = 2.39, a= 0.14; 
et es TOS 
~ 2.39 x + 0.14 


x — 0.65 


Vie ODO. 0.14 + 2.39% or y 


+ 0.120. . 


Arr. 79 - THE LOGARITHMIC CURVE, log y =a-+ bx + cx® I51 


We have computed both y and y’ and note that the agreement with 
the observed values is probably as close as could be expected. 

79. The logarithmic or exponential curve, log y = a + bx + cx? 
or y = ae’”*+°”, — These equations are modifications of the logarithmic 
form log y = a + bx and the exponential form y = ae”. The equation 
y = ae=+e may be written log y = log a + (6 log e) »+ (c log e) x’, and 
so is equivalent to the form log y = a + bx,+ cx®. ‘This last equation 
is similar in form to the equation y = a + bx + cx? discussed in Art. 77, 
and the equation may be verified and the constants determined in a 
similar way. 

Hence, if a set of data may be approximately represented by the equation 
logy =a+ bx + cx’, then (1) the plot of (x, ey Eh) where 
(xz, Vx) are the codrdinates of a point on the experimental curve, will approxt- 
mate a straight line, or (2) the plot of (x, A log y), where the A log y are the 
differences in log y formed for equidistant values of x, will approximate a 
straight line. 

Example. The following table gives the results of Winkelmann’s 
experiments on the rate of cooling of a body in air; @ is the excess of tem- 
perature of the body over the temperature of its surroundings, ¢ seconds 
from the beginning of the experiment. 


t 0 log é log 6 — log 118.97 1088 — Tog 138.97 6; A 

° 118.97 2.07544 ° 118.97 ° 
Ue at 116.97 2.06808 —0.00736 —o0.000608 116.99 —0.02 

251.8 II4.97 2.06059 —0.01485 —0.000576 114.97 ° 
Ain7, II2.97 2.05296 —0.02248 —0.000539 112.90 +0.07 
SORT II0.Q7 2.04520 —0.03024 —0.000507 110.90 | +0.07 
82.0 108 .97 2.03731 —0.03813 —0.000465 108 .go +0.07 
109.0 106.97 2.02926 —0.04618 —0.000424 107.15 — yan its) 


In Fig. 79 we have plotted the points (¢, @). According to Newton’s 
law of cooling, 6 = ae” or log 9 = a + Ot, and so we have also plotted the 
points (¢, log 6); this last plot has a slight curvature. We shall therefore 
assume the law in the form log @=a-+bt+c#. To verify this, we 
note the point t = 0, 0, = 118.97 on the experimental curve, and plot 
cen “8 eet); this plot approximates a straight 

log 0 — log 118.97 _ 


line, so that we may assume ; =b+ct. We use the 


method of averages to determine the constants. Dividing the data into 
two groups of three sets each and adding, we get 

—0.001723 = 3b+ 79.66, 

—0.001396 = 30 + 250.7 ¢. 


the points (« 
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C = 0.0000019II, b = —0.000625. 
log 6 — log 118.97 _ 
t 
or log 8 = 2.07544 — 0.000625 ¢ + 0.0000019I11 #. 
We now compute @ and the residuals and note the close agreement 
between the observed and calculated values. 


— 0.000625 + O.00000IQ9II t 


a 


log 118 


t 


( log 6— 


S 
Ss 
S 
Q) 


0 10 20 30 40 50 és 70 80 90 100 110 


FIG. 79. 


(IV) EQUATIONS INVOLVING FOUR OR MORE CONSTANTS. 


80. The additional terms ce” and cx*.— It is sometimes found that 
a simple equation will represent a part of our data very well and another 
part not at all, z.e., the residuals yy — y, are very small for one part of 
our data and quite large for another part. Geometrically, this is 
equivalent to saying that the plot of the simple equation coincides 
approximately only with a part of the experimental curve. In such 
cases a modification of the simple equation by the addition of one or 
more terms will often cause the curves to fit approximately throughout. 
Such terms usually have the form ce? or cd*, and added to our simple 
equations give the forms 


y=a-+ bx + ce*, y=a+t bx + cx4, 
Aye de Ge, y = ax cx", 

= Os dz = x d 

y ac hy sine at y Tie etc. 


We shall give a few examples to illustrate some of these cases. 
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81. The equation y = a+ bx + ce”.—If a part of the experimental 
curve approximates a straight line, we may fit an equation of the form 
y =a-+ bx to this part of the curve. The deviation of this straight 
line from the remainder of 
the experimental curve (Fig. 
81a) will be measured by the 
residuals 7 = yo — y. = y — 
(a + bx). We now plot (x, 7) 
and study the nature of this 
plot. We may be able to 
represent this plot by means 
of the simple exponential r = 
ce?*, where the values of the 
constants c and d are such 
that the value of + is negli- 
gible for that part of the plot 
to which the straight line has 
been fitted. The entire ex- 
perimental curve can thus be 
represented by ce’ = y — (a ieee 
+ bx) ory =a-+ bx + ce®. 

The equation y = a+ bx + ce# may fit an experimental curve 
although no part of the curve is approximately a straight line; this 
means that the values of the term ce“ are not negligible for any values 
of x. If the values of x are equidistant, we may verify that this equa- 
tion is the correct one to assume by the following method. Let the 
constant difference in the values of x be h. If we replace x by x +h, 
we get 


y =atb(ut+h) + cere, 
and, therefore, for the difference in the values of y, 
Ay = y' — y = bh + ce* (e# — 1). 
If Ay and Ay’ are two successive values of Ay, then 
Ay’ = bh + cet@+ (e% — 1), 
and the difference in the values of Ay is 


A?y = Ay’ — Ay = ce* (e** — 1)?, 


Hence, 
log A?y = log c (e* — 1)? + (dloge) x. 


The last equation is of the first degree in x and log A’y so that the plot 
of (x, log A?y) will approximate a straight line. From this straight 
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line we may determine the constants log c (e# — 1)? and dloge and 
therefore c and d in the usual way. We now write the equation in the 


y=atdbet+ co 


form y — ce” =a -+ bx, and 
from the straight line plot of 
(x, y — ce’), we determine the 
constants a and b. 

In Fig. 816 we have plotted 
the equations 


y=0.5 +x, 

Viz 0:5 wk — O01, 
y = 0.5 + x — 0.001 e, 
yV=0.5 +274 0.016, 
y= 0.5 + 4+ -0.00L-¢. 


Example. The following 
data are the results of experi- 
ments made with a gasometer 
by means of which the amount 
of air which passes into a re- 
ceiving tank can be measured; — 
x is the vacuum in the tank in 
inches of mercury, y is the 
number of cu. ft. of air per 


} minute passing into the tank. 


(Experiments made by W. D. 
Canan at the Mass. Inst. of 


Fic. 810. Tech.) 

Zz y y' r=y'—y log r Te Ye 4 

8 1.17 1.49 0.32 9.5051 —10 | 0.322 I.17 Giri 
se) 37 1.55 0.18 9.2553 —1I0 | ©.179 Te Si ° 

12 I.50 1.61 O.II 9.0414 — IO ©.099 Tee Sa —O.O 
14 62 Oy 0.05 8.6990 — 10 | 0.055 1.61 +0.01 
16 1.71 Te76 0.02 0.031 Ti7O —-OLOr 
18 1.80 E79 —0O.01 0.017 TES viy +0 .03 
20 1.85 1.85 ° © .009 1.84 =1-O).OL 
22 I.QI I.QI ° ©.005 I.go =OLOk 
24 I.96 1.97 ©.O1 0.003 1.97 —0.0r 
26 2.02 2.03 0.01 © .002 2.03 —0:0L 
28 2.10 2.09 —0 OL ©.0O1 2.09 -++o.o1 


In Fig. 81c we note that the plot of (x, y) approximates a straight 
line for values of x > 14, and we shall fit an equation of the form 
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y’ =a + bx to this part of the data. Using the method of averages 
and dividing the data into two groups of four and three sets, we have 
7-27 =4a+ 76), 
6.08 = 3a+ 78), 
b= 0.03, @= 1.25 
and yy’ = 1.25 + 0.03%. 


eel [af] |S SO eT SNE 
Deane Sees 


ao 


& 


EU (ee TF ee | 


18 
(@) 
FiG.., 81c¢. 

Now compute the values of y’ and the residuals r = y’ — y (by 
taking ry = y’ — y instead of r = y — y’, the residuals are positive and 
easier to handle in the subsequent calculations). Plot («, r) for values 
of x < 14 and study the nature of this plot; this seems to be a simple 
exponential, 7 = ce**; verify this by plotting (x, log r) and note that 
this plot approximates a straight line. Using the method of averages 
determine the constants in the equation log r = log c + (d loge) x; thus 

8.7604 — 10 = 2loge + 18dloge, 
7.7404 — 10 = 2logc + 26d loge. 
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d log e = 9.8725 — 10 = —0.1275, loge = 0.5277. 
d= — 0.204, 16 = 13:37: 
log Y = 0.5277 — 0.1275 %, and r= 2.37, e70.2942 


The final equation is 
y = 1.25 + 0.03% — 3.37 € 0". 


Now compute y and the residuals, and note the close agreement be- 
tween the observed and calculated values. 

82. The equation y = ae’” + ce*”.—A part of the experimental 
curve may be represented by a simple exponential y = ae”, 1.¢., a part 
of the plot of (x, log y) approximates a straight line. We then study 
the deviations, r = yo — ye = y — ae’, of this exponential curve from 
the rest of the experimental curve. The plot of (x, 7) may be repre- 
sentable by another exponential, r = ce”, where the values of 7 are 
negligible for that part of the experimental curve to which y = ae’? has 
been fitted. The entire curve can then be represented by the equa- 
tion y = ae + ce, 

The equation y = ae** + ce may fit an experimental curve although 
no part of the curve can be approximated by the simple exponential 
y = ae, If the values of x are equidistant, we may verify that this 
equation is the correct one to assume by the following method. Let 
the constant difference in the values of x be #. Consider three succes- 
sive values x, x + h, x + 2h and their corresponding values y, y’, y’’. 
We evidently have 

y =a ae" er, 
yy! = aed (th) 4 ced(eth) = aber! + cetzeth, 
yy!’ = aeb@t2h) A ced (et2h) = gedrg2rh 4 codzg2dh, 
Now eliminate e’* and e** from these three equations by multiplying 
the first equation by e@t%’, the second by —(e + e#), and adding the 
results to the third equation. We get 
yl" — (eh + eft) y! 4+ ebtdhy = 0, 
yl : y’ 
or OL me (eb ob em) ie e(otd) h 
Bf y. 
This is an equation of the first degree in y’/y and y’’/y so that the plot 
of (y’/y, y”’/y) will approximate a straight line. From this straight line 
determine the constants e’ + e” and e@+®*, and hence b and d as usual. 
We now write the original equation ye = ae@-%*+-¢, This is a linear 
equation in e@-* and ye so that the plot of (e@-9+, ye-@) would 
approximate a straight line. From this straight line determine the 
values of the constants a and c. 
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In Fig. 82a, we have plotted the equations y = e~*, y = e-* +05 75%, 
y= Ca 0.5 Cae y= ew a Gre Vee ee? — e722, 


2.0 
1.8 
1.6 
1.4 
2 
i = Seis “5 
\ & 
Q\ (ee 
(y) 1.0 —\ 
ce can 
0.8} SSMS, 
0.6 I~ o> 
"ON 
0.4 | ar 
0.2 @1 =P 
a) = 
0 
0 0.5 1.0 1.5 2.0 2.5 
(x) 
= ae ae cel® 
Pie. 82a, 


Example. The following are the measurements made on a curve 
recorded by an oscillograph representing a change of current 7 due to a 
change in the conditions of an electric circuit ¢. (From Steinmetz, 
Engineering Mathematics.) 


t a log z OM r=v-i log r Ty OF A 
° 2 LOe Os 222 4.94 2.84 0.4533 2.85 2.09 +o0.01 
One 2.48 | 0.3945 4.44 1.96 0.2023 1.96 2.48 ° 
0.2 2.66 | 0.4249 3.99 1 2 0.1239 1.34 2.65 +0.01 
0.4 2.58 | 0.4116 Bee 0.64 g .8062 — I0] 0.63 2.59 —0.01 
0.8 | 2.00 | 0.3010 2.10 0.10 Q .0000 — IO} 0.14 1.96 +0 .04 
2 i307 O11335 Tey ©.O1 0.03 1.34 Fin. O2 
ELON |, O4907|.0).9542)—10|' 10.89 —O1OL 0.01 0.88 +0.02 
QAOme O25 & \\e9n 7034) — 10 0.58 ° ° 0.58 ° 
aks || Ciszek Wojciak)» Koleeysl ° ° 0.34 ° 
BE! 80).20' | 19-2010 — TO|| 0.20 ° ° 0.20 ° 

| 


In Fig. 826 we note that the right-hand part of the plot of (¢, 7) appears 
to be exponential. We verify the choice of 2’ = ae” by plotting (¢, log 2) 
and noting that this plot approximates a straight line for values of 
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t> 0.8. We therefore assume log i’ = log a + (0 log e) ¢, and using the 
method of averages for the values of t > 0.8, we have 
9.8511 — 10 = 3loga+ 4.8 d loge, 
8.8325 — 1o = 2loga+ 5.5 db loge. 
blog e = 9.5356 — 10 = —0.4644, loga = 0.6934, 
b= -—1.07, &@ = 4.94, 
and log i’ = 0.6934 — 0.464414, or 1 = 4.94 eh! 


seeseces cat 
ALI op eS 
iil 


log 

ARBRE Ree oss 

ea es tata 
2.0 : 3.0 


. Fic. 820. 


Now find the values of 7’ and the residuals r = i’ — 7; these residuals 
are practically negligible for values of t > 0.8. We plot (¢, r) and try 
to fit an equation to this curve. This again appears to be exponen- 
tial and we verify this by plotting (¢, log 7); the plot approximates a 
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straight line, except for ¢ = 0.8. We therefore assume r = ce” or log r = 
log c + (dloge) t. Using the method of averages for ¢ < 0.8, we have 
0.7456 = 2logc + o.1 dloge, 
9.9301 — 10 = 2logc + 0.6dloge. 
dloge = —1.6310, logc = 0.4544. 
“id= —3.76,. € = 2.85; 

and logr = 0.4544 — 1.63101, or 7 = 2.85 e% 7, 
The final equation is 


I 


co 4.94 e 1.07 t ewe 2.85 E7876 E | 

We now compute 7 and the residuals and note the very close agree- 
ment between the observed and computed values of 7. 

83. The polynomial y = a + bx + cx? + dx’ + - - - .—The equa- 
tion y = a+ bx + cx® may be modified by the addition of another 
term into y = a + bx + cx? + dx*. If the values of x are equidistant, 
we may verify the correctness of the assumption of the last equation 
by the following method. Let the constant difference in the values of 
x beh. Then the successive differences in the values of y are 

Ay = (bh + ch? + dh’) + (2ch + 3. dh?) x +3 dhx?, 

A’y = (2ch? + 6dh?) + 6dh’x, ~ 

My = 6 dh’. 
Hence the plot of (x, A’y) will approximate a straight line, and the values 
of A*y are approximately constant. From the equation of the straight 
line we may determine the constants c and d, and writing the original 
equation in the form (y — cx? — dx’) = a + bx, the plot of (x, y —cx?— dx?) 
will approximate a straight line, from which the constants @ and b 
may be determined. Another method of determining the constants 
a, 6, c, d in the equation y = a + bx 4+ cx? 4+ dx} consists in selecting 
four points on the experimental curve, substituting their codrdinates 
in the equation, and solving the four linear equations thus obtained for 
the values of the four quantities a, b, c, and d. 

In a similar manner the polynomial y =a-+ bx -+cx? ++ +--+ + kx" 
may be determined so that the corresponding curve passes through 
nm +1 points of the experimental curve; it is simply necessary to sub- 
stitute the codrdinates of these m + I points in the equation and to solve 
the nm + 1 linear equations for the values of the m + 1 quantities, a, B, 
c,...,k. Ifthe values of x are equidistant, we can show that the plot 
of (x, A*-1y) is a straight line and that A”y is constant, where A™-ly 
and A”y are the (x — 1)st and nth order of differences in the values of y. 
Thus, if a sufficient number of terms are taken in the equation of the 
polynomial, this polynomial may be made to represent any set of data 
exactly; but it is not wise to force a fit in this way, since the deter- 
mination of a large number of constants is very laborious, and in many 
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cases a much simpler equation involving fewer constants may give much 


more accurate results in subsequent calculations. 
We shall work a single example to illustrate the method of determin- 


ing the constants. 
Example. We wish to fit a polynomial equation to the following 


data: 


ase 
pei 
ee 
ee 
Pe 
Selo 
pA 
Paina 
ee 
ae 
ale 


0 1 0.2 30.3 0546-5" 0:6 057, 0:8 10:9 7.0" 
kB) 
Fic. 83 
z y Ay A’y A’y Ye a 
° ° 0.212 ©.039 0.019 ° ° 
O.1 ©.212 0.251 0.058 0.014 0.210 +0 .002 
0.2 © .463 0.309 ©.072 ©.019 0.463 ° 
0.3 On772 0.381 ©.0QgI 0.019 0.770 * +0 .002 
0.4 Teese 0.472 ©.110 0.018 PeL5s +0.00o1 
0.5 1.625 0.582 0.128 0.021 I .625 ° 
0.6 20207 0.710 0.149 0.014 2.209 —0.002 
0.7 2.917 0.859 0.163 °.018 2.920 —0.003 
0.8 3.776 I .022 0.181 3.776 ° 
0.9 4.798 I .203 4.797 +0 .001 
1.0 6.001 5.998 +9 .003 
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In Fig. 83 we have plotted («, y). We form the successive differences 
and note that the third differences are approximately constant, and 
that the plot of (x, A’y) approximates a straight line (Fig. 83). We 
may therefore assume an equation of the form y = a + bx + cx? + dx’, 
or y = bx + cx’ + dx’, since the curve evidently passes through the 
origin of codrdinates. To determine the constants b, c, and d, select 
three points on the experimental curve; three such points are (0.2, 0.463), 
(0.5, 1.625), and (0.8, 3.776). Substituting these codrdinates in the 
equation, we get 

0.463 = 0.2b + 0.04.¢ + 0.008 d, 

1.625 = 0.56+0.25c¢ + 0.125 d, 

3.776 = 0.86 + 0.64¢ + 0.512 d. 
Solving these equations for b, c, and d, we have 

b= 1.989, ¢ = 1:037, d = 2:972 
and hence the equation is 

y = 1.989 x + 1.037 x? + 2.972 x*. 
We now compute the values of y and the residuals. 

84. Two or more equations. — It is sometimes impossible to repre- 
sent a set of data by a simple equation involving few constants or even 
by a complex equation involving many constants. In such cases it is 
often convenient to represent a part of the data by one equation and 
another part of the data by another equation. The entire set of data 
will then be represented by two equations, each equation being valid 
for a restricted range of the variables. Thus, Regnault represented the 
relation between the vapor pressure and the temperature of water by 
three equations, one for the range from —32°F. to 0° F., another for 
the range from 0° F. to 100° F., and a third for the range from 100° F, 
to 230° F. Later, Rankine, Marks, and others represented the rela- 
tion by a single equation. The following example will illustrate the 
representation of a set of data by two simple equations. 

Example. The following data are the results of experiments on the 
collapsing pressure, p in pounds per sq. in. of Bessemer steel lap-welded 
tubes, where d is the outside diameter of the tube in inches and ¢ is the 
thickness of the wall in inches. (Experiments reported by R. T. Stew- 
art in the Trans. Am. Soc. of Mech. Eng., Vol. X XVII, p. 730.) 


4 t 

7 P log 5 log P Wee A 
0.0165 225 8.2175 — 10 2) UY: 230 —5 
0.0194 383 8.2878 — Io 2.5832 381 +2 
0.0216 524 8.3345 — 10 2.7193 533 —9 
©.0214 536 8.3304 — 10 2.7292 517 +19 
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% 0.01 0.02 0.03 0.04 ce hoe 0.07 0.08 0.09 0.10 


/a 
Fic. 84. 
t t 
. P P, A Q P P A 
0.0228 592 570 ae 0.0370 1779 1821 = Ale 
©.0250 670 764 — 84 0.0374 1860 1856 + 4 
0.0253 870 790 ap ee 0.0375 1879 1865 + 14 
0.0277 928 1002 —) 74 © .0400 2147 2085 = 62 
0.0298 964 1187 —= 222 © .0403 2224 2112 ir? 
0.0299 1184 1196 — 12 0.0436 2280 2403 —123 
© .0309 T2250 1284 == ee} ©.0442 2441 2455 — 14 
0.0316 1319 1346 Ba 0.0477 2962 2764 +198 
0.0309 I419Q 1284 SSeS 0.0527 3170 3204 — 34 
0.0343 1680 1583 +197 0.0628 4095 4194 — 99 
0.0349 1762 1636 +126 0.0815 5560 5741 —181 
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It should be noted that a set of corresponding values of t/d and P 
are not the results of a single experiment but the averages of groups 
containing from two to twenty experiments. 

Following the work of Prof. Stewart, we have plotted (¢/d, P), Fig. 84, 
and note that the experimental curve approximates a straight line for 
all values of ¢/d except the first four, 7.e., for values of t/d > 0.023. 


We may therefore assume P = a + »(;). If we use the method of 


selected points to determine the constants a and b we may choose the 
points t/d = 0.065, P = 4250, and t/d = 0.030, P = 1215 as lying on 
the straight line; we then have 

4250 = a+ 0.065), 

1215 =a+0.0300. 


b = 86,714, = —1386 
and P = 86,714 (5) — 1386. 


This result agrees with that given by Prof. Stewart. If we use the 
method of averages to determine the constants a and 6 we divide the 
last 22 sets of data into two groups of I1 each, and get 


12,639 = 11 ad + 0.3231 0, 
30,397 = I1a + 0.5247 b. 
b = 88,085, a= —1438, 


L 
and P = 88,055 (5) — 1438. 


In our table we have given the values of P computed from this last 
formula. The values of P computed from the first formula agree very 
closely with these. It is seen that the percentage deviations are in 
general quite small though large in a few cases, varying from 0.2 per 
cent to Io per cent, which is to be expected from the nature of the 
experiments. 

We now attempt to fit an equation to the first four sets of data. The 


, aN 
addition of a modifying term of the form (5) or ce . to the above 


formula is not successful here. We shall therefore follow Prof. Stew- 
art’s work and attempt to fit an equation of the parabolic form, P = 
b : t : 
a(5) . We verify this choice by plotting (Ioe 7 log P) and observing 
that this plot approximates a straight line. (The fewness of the ex- 


periments for values of t/d < 0.023 is a handicap here.) Assuming 


log P = loga+ blog(5); 
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and using the method of averages, we find 
4.9354 = 2 loga + (6.5053 — 10) 8, 
5.4485 = 2 loga + (6.6649 — 10) b. 
= 3.11, @ = 80,580,000 
t 3.11 
and P = 80,580,000 (;) : 
We compute the values of P from this formula. 


The entire set of data have thus been represented by means of two 
simple equations, each valid for a restricted range of the variables.* 


EXERCISES. 


[Note. The exercises which follow are divided into two sets. The type of equation 
that will approximately represent the empirical data is suggested for each example in 
the first set. For the examples in the second set, the choice of a suitable equation is left 
to the student.] 

1. Temperature coefficient; 7 is the resistance of a coil of wire in ohms, @ is the tem- 
perature of the coil in degrees Centigrade. [y = a + bx] 

10.939 DL.22 | 11.799 | 12.242 12.668 
29.49 42.70 60.01 75.51 QI.05 

2. Galvanometer deflection; D is the deflection in mm., J is the current in micro- 
amperes. [Y = a + bx] 

IB) psy 92.0 | 118 | 140 165 

‘I | 0.0493 0.154 OslO 7a). 0.234 0.274 


3. Volt-ampere characteristic of 118 volt tungsten lamp; e is the terminal voltage, z 
is the current. [y = ax®] 


r 


) 


10.421 
10.50 


48.2- 
0.0821 


72.7 
0.123 


199 
0.328 


| 


é 2 4 8 16 25 32 50 64 100, | 125 
1 | 0.0245 | 0.0370 | 0.0570 | 0.0855 | 0.1125 | 0.1295 | 0.1715 | 0.2000 | 0.2605 | 0.2965 
e 150 180 200 218 
4 | 0.3295 | 0.3635 | 0.3865 | 0.4070 
e Pressure-volume of saturated steam; v is the volume in cu. ft. of 1 pound of 
stedm, p is the pressure in pounds per sq. in. [y = ax®] 
v | 26.43 22.40 19.08 16.32 14.04 | 12.12 10.51 | 9.147 7.995 
p | 14.70 17.53 20.80 24.54 28.83 33570 39.25 45-49 52.52 


5. Chemical concentration experiment; x is the concentration of hydrogen ions, y is 
the concentration of undissociated hydrochloric acid. [y = ax®] 
1.68. 0.784 0.047 
D382 0.216 0.00315 


x 
Sy, 


1.22 
0.676 


0.426 
0.074 
6. Vibration of a long pendulum; A is the amplitude in inches, ¢ is the time since it 
was set swinging. [y = ae>a] 
MT meres oe 
4.97 2.47 1.22 0.61 0.30 0.14 


0.092 
0.0085 


0.0096 | 0.0049 
0.00036 | 0.00014 


0.00098 
0.000018 


oO 
10 


A 


* Prof. Peddle in ‘‘ The Construction of Graphical Charts”’ has fitted the equation 
t/d = 0.00274 VP + 0.0000000011 P? to Prof. Stewart’s data. 
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ewton’s law of cooling; @ is the excess of the temperature of the body over the 
temperature of its surroundings, ¢ is the time in seconds since the beginning of the experi- 
ment. [y = ae>~] 
t | 0) | 3-45 


ads 10.85 | 19.30 
19.9 18.9 


16.9 | 14.9 


28.80 | 40.10 | 53-75 | 70.95 


0 12.9 10.9 8.9 6.9 


8. Barometric pressure; p is the pressure in inches of mercury, # is the height in ft. 
above sea level. [y = ae°*] 


A o | 886 | 2753 | 4763 6942 10,593 
30 29 27 25 23 20 


P 
g. Electric arc of length 4 mm.; V is the potential difference in volts, 2 is the current 


in amperes. [» =a+t =| 


7:97 
54.30 


a | 2.46 | 2.97 | 3-45 | 3-96 | 4-97 5-97 6.97 
V 67.7 65.0 63.0 61.0 58.25 56.25 55-10 


10. Speed of a vessel; H.P. is the horse power developed, v is the speed in knots. 
[y = 4 + bx’) 
12 
2300 


v 5 | is 9 UI 
EP 290 560 1144 1810 


11. Hydraulic transmission; H.P. is the horsepower supplied at one end of a line of 


pipes, u is the useful power delivered at the other end. [z =a+t bx:| 


Jee 100 150 200 | 250 | 300 
uU 96.5 138 172 196 206 


12. Magnetic characteristic of iron; H is the number of gilberts per cm., a measure 
of the field intensity, B is the number of kilolines per sq. cm., a measure of the flux 


densi ~s | 

ensity. [> =a 
Jgf | 8 | 10 | 15d ager O | 30 | 40 | 60 80 
B 13.0 14.0 TS Autos 17.2 17.8 Tin || munaoee 


13. Focal distance of a lens; p is the distance of the object, p’ is the distance of its 
; x 
image. [» = | 


p 
p 


80 
26.20 


140 
23.20 


120 
23.80 


100 
24.60 


60 
29.00 


320 240 | 180 
21.35 21.80 22.50 


14. Pressure-volume in a gas engine; p is the pressure in pounds per sq. in., v is the 
volume in cu. ft. per pound. [y = ax? + ¢] 


p | 44.7 
v 7.03 


135.8 
1.90 


53.8 | els | 85.8 | 113.2 
5.85 4.30 3.50 2.50 


15. Law of cooling; 6 is the temperature of a vessel of cooling water, ¢ is the time in 

minutes since the beginning of observation. [y = ae>” + c| 
een ean oe ae 
92.0 74:5 67.0 60.5 53-5 45.0 


I 
85.3 


2 
79-5 


20 
39-5 


0 
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16. Straw-fibre friction at 150 pounds pressure according to Goss’s experiments; y 
is the coefficient of friction for a straw-fibre driver and an iron driven wheel, « is the slip, 


x 
per cent. [y = ax + c| 


x |0.153 |0.179 | 0.213 | 0.271 | 0.313 | 0.359 | 0.368 | 0.381 | 0.386 | 0.405 
rye 0.56 | 0.58 0.61 0.78 0.99 1.10 1.04 1.22 1.40 1.75 
_% | 0.411 0.432 | 0.458 | 0.463 0.465 | 0.473 
“y 1.94 2.00 2.25 2.33 3.15 2.79 


17. Expansion of mercury according to Regnault’s experiments; y is the coefficient 
of expansion between 0° C. and#?C. [y=a-+ bx + cx?| 
t 250 


Ee 0.00018403 


200 
0.00018323 


100 
0.00018216 


fo) 
0.00018179 


150 
0.00018261 


300 | 360 
0.00018500 |0.00018641 


18. Velocity of water in Mississippi River; v is the velocity, D is the depth. 
[y =a-+ bx + cx?] 
D 


v 


0.6 
3.1807 


Co) 0.1 0.2 
3.1950 |3.2299 | 3.2532 
19. Solution of potassium chromate; s is the weight of potassium chromate which will 
dissolve in 100 parts by weight of water at a temperature of °C. [log y=a+bhx-+cx?] 


27.4 
66.3 


0.4 0.5 
3.2516 | 3.2282 


0.3 
3.2611 


On7, | Osh, S| Ose 
3-1266 | 3.0594 | 2.9759 


10 
62.1 


t 
s 


) 
61.5 


42.1 
70-3 3. 


20. Load-elongation of annealed high carbon steel wire of diameter 0.0693 and gage 
length 30 in.; W is the load in pounds, E is the elongation ininches. [y=a+bx+ce@] 


Wo 50 100 150 200 225 250 260 280 290 | 300 310 
E | 0 |0.0130 |0.0251 |0.0387 |0.0520 [0.0589 |0.0659 |0.0689 |0.0746 |0.0778 |0.0807 | 0.0842 
Wi. £320 330 340 350 360 
E | 0.0877 | 0.0916 | 0.0980 | O.11II | 0.1420 


21. Load-elongation of wire of Ex. 20 in hard-drawn condition; W is the load in 
pounds, £ is the elongation in inches. [y = a+ bx + cx4] 


Ww (o) 100 200 | 300 400 500 600 700 800 850 900 

E 0 | 0.0280] 0.0562] 0.0849] 0.1150] 0.1471| 0.1820] 0.2191] 0.2628] 0.2879] 0.3166 
22. Empirical curve. [y = ae>© + cede] 

x ) 0.3 | 0.6 0.9 nee: 1.5 | 1.8 25 2.4 3.0 

y | 3.00 | 1.89 | 1.27 0.88 0.63 0.46 0.33 0.25 0.18 0.10 


23. Magnetic characteristic of iron; H is the number of gilberts per cm., a measure 
of the field intensity, B is the number of kilolines per sq. cm., a measure of the flux density 


(cf. Ex. 12). [» = == + cet | 


6 
E1e2 


8 
13.0 


eI 
14.0 


4 
8.4 


15 
15.4 
24. Speed of a vessel; J is the indicated horsepower, v is the speed in knots. 
[y=a+bx+ cx’ + dx] 
8 9 ae) 
1000 | 1400 | 1900 


20 30 
LOL ae ee 


40 | “60 80 


| 2 
B |. 3.0 17.8) 118.514) es 


v 
I 


I2 


3250 


II 
2500 


13 


4200 


14 15 


5400 | 6950 


16 | 17 | 18 
8950 |11,450] 15,400 
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25. Test on square steel wire for winding guns; S is the stress in pounds per sq. in., 
E is the elongation in inches per inch. 


S | 5000 | 10,000| 20,000 | 30,000 | 40,000 | 50,000 | 60,000 | 70,000 | 80,000 
E © |0.00019| 0.00057 | 0.00094 | 0.00134 |,0.00173 | 0.00216 | 0.00256 | 0.00297, 
S__| 90,000 | 100,000 |110,000 . 
E | 0.00343 | 0.00390 |0.00444 


26. Flow of water over a Thomson gauge notch; Q is the number of cu. ft. of water, 
H is the head in feet. 


H m2 | 1.4 | 1.6 | 1.8 2.0 2.4 
Q 4.2 6.1 8.5 11.5 . 14.9 23.5 


27. Friction between belt and pulley; @ is the arc of contact in radians between belt 
and pulley, P is the pull in pounds applied to one end of pulley to raise a weight W at the 
other end. 


8 7 2 51 7 47 30 51 Ilr 
us ase) 2% a lege 6 ol) ot vey pa! 
2 3 6 6 B 2 3 6 
ie 5.62 6.93 8.52 10.50 12.90 15.96 19.67 24.24 29.94 


28. Electric arc of length 2 mm.; V is the potential difference in volts, 7 is the cur- 
rent in amperes. ; 

4 

V 


1.96 
50.25 


2.46 
48.70 


6.97 
44.00 


2.97 
47.90 


7.97 9.00 . 
43-60 | 43-50 


3-45 
47.50 


5:97 
45.00 


3-96 | 4.97 
46.80 | 45.70 


29. Normal induction curve for transformer steel; H is the number of gilberts per 
cm., B is the number of lines per sq. cm. 


A eEO 2.3.) 52.5 2.9 3.4 4.1 4.5 BA) "sho | ies |. OO 11.0 
B | 425 | 800} 1750 | 2850 | 4300 | 6100 | 6725 | 7800 8600 }10,200|11,150] 12,200 
30. Pressure-volume in a gas engine; p is the pressure in pounds per sq. in., v is the 
volume in cu. ft. per pound. , 
p 39.6 44.7 53.8 | 73.5 85.8 113.2 135.8 178.2 
v 10.61 9.73 8.55 7.00 6.23 5.18 4.59 3.87 


31. Melting point of alloy of lead and zinc; 6 is the temperature in degrees Centi- 
grade, x is % of lead. 
40 50 60 70 80 90 
ere ee eae ae 


x 
0 


32. Empirical curve. 
3 5 if | 9 | II | 13 
8.50° 11.03 14.03 17.53 21.55 AS 


xe I 
y 6.42 


33. Candle-power of an incandescent lamp; H is the age of the lamp in hours, C.P. 
is the candle-power. 
H 
POP: 


O) 
24.0 


1500 
14.5 


250 | 500 | 750 | 1000 | 1250 
17.6 16.5 15.8 15.3 14.9 


34. Insulation resistance-current passes through insulator and galvanometer; D is 
the deflection of the galvanometer, ¢ is the time in minutes. 


t 6 7 8 9 


D 5.0 | 4.4 | 4.0 | 3.5 


12 
2.7 


to | Ir 


3:31) 3:0 


15 
2.4 


I 


138 


2 


II 


13 
2.5 


14 
2.5 


8.0 


6.2 


5:5 
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35. Experiments with a crane; f is the force in pounds which will just overcome a 
weight w. 
w 100 | 200 | 300 | 400 | 500 
if | 8.5 12.8 17.0 21.4 25.6 
36. Copper-nickel thermocouple; ¢ is the temperature in degrees, p is the thermo- 
electric power in microvolts. 
t (0) 50 200 
p 24 25 Dyfi 
37. Law of falling body; s is the distance in cm. fallen by body in ¢ sec. 
t | 0.2477 | 0.4175 0.5533 0.6760 | 0.7477 
§ 30.13 85.26 150.39 223.60 274.20 
38. Loads which cause the failure of long wrought-iron columns with rounded ends; 
P/a is the load in pounds per sq. in., //r is the ratio of length of column to the least 
radius of gyration of its cross-section. 
L/r 180 
P/a 7500 


39. Heat conduction of asbestos; @ is the temperature in degrees Fahrenheit, C is 
the coefficient of conductivity. 

6 32 | 212 392 | 572 | 752 

G 1.048 1.346 1.451 1.499 1.548 


40. Rubber-covered wires exposed to high external temperatures; C is the maximum 
current in amperes, A is the area of cross-section in sq. in. 


600 | 700 | 800 
29.9 34.2 38-5 


100 150 
26 26.9 


140 
12,800 


300 340 380 | 420 
2800 2100 1700 1300 


220 | 260 
5000 3800 


a 
1.644 


Cc Be 5-9 9.0 22.0 42.0 68.0 84.0 102.0 
A_ | 0.001810 | 0.004072 | 0.007052 | 0.02227 | 0.05000 | 0.09442 | 0.1250 | 0.1595 
41. Pressure-volume relation for an air compressor; # is the pressure, v is the volume. 
p 18 21 | 26.5 | 2205 | 44 62 
v 0.635 0.556 0.475 0.397 0.321 0.243 


42. Power delivered by an electric station; w is the average weight of coal consumed 
per hour per kilowatt delivered, f is the load factor. 
f | 0.25 | 0.20 
w 2.843 3.012 


0.15 | 0.10 | 0.05 
3-293 3.856 5-545 

43. Temperature at different depths in an artesian well; @ is the temperature in 
degrees C., d is the depth. 


d | 28 66 
6 iegAt 12.90 


173 
16.40 


248 | 298 | 400 505 548 
20.00 22.20 23.75 26.45 27.70 


44. Resistance of copper wire; R is the resistance in ohms per 1000 ft., D is the 
diameter of wire in mils. 
D 289 
R 0.126 


182 
0.317 


102 
I.O1O 


57 
3.234 


45. Hysteresis losses in soft sheet iron subjected to an alternating magnetic flux; 
B is the flux density in kilolines per sq. in., P is the number of watts lost per cu. in. for 
I cycle per sec. 
B | 20 
ip 0.0022 


10 
105.1 


32 18 
10.26 32.80 


60 
0.0128 


80 
0.0202 


40 
0.0067 


100 
0.0289 


120 
0.0387 
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46. Volt-ampere characteristic of a 60 watt tungsten lamp: V is the number of 
volts, J is the number of milli-amperes. 


V 


2 5 10 20 30 40 50 60 70 80 |; 90 | 100 
I 49 | 80 | 117 | 180 | 227 | 272 | 311 | 348 | 383 | 414 | 443 | 473° 
V | 110 | 120 | 130 | 140 | 150 | 160 | 170 | 180 | 190 | 200 | 210 | 229 
T | sor | 526 | 553 | 577 | 597 | 618 | 639 | 663 | 682 | 702 | 722 | 743 


47. Calibration of base metal pyrometer (40% Ni and 60% Cu); V is the number 
of millivolts, ¢ is the temperature in degrees F. 
4 4 6 14 16 
255 320 


o— | 
ce) 146 


8 | 10 | 12 | 
t 396 475 553 634 714 
48. Tests on drying of twine; ¢is the drying time in minutes (time of contact of 
twine with hot drum), W is the percentage of total water on bone dry twine at any 
time, E is the percentage of total water on bone dry twine at equilibrium, d is the di- 
ameter of the twine in ins. 


(a) d = 0.102 ins., E = 18.7%. 


t | (0) 0.44 0.88 1.31 | LSE. 
W-E 29.5: | 15.4 | 9.4 | 5.1 oa 
(bd = 0:58, E-= 62%. 
t | oO PS sere | 2.23 | 3.34 | 4.45 | 5.56 
W-E 30.3 | 17.4 12.4 8.2 4.9 23 


CHAPTER VII. 
EMPIRICAL FORMULAS — PERIODIC CURVES. 


85. Representation of periodic phenomena. — Periodic phenomena, 
such as alternating electric currents and alternating voltages, valve-gear 
motions, propagation of sound waves, heat waves, tidal observations, 
etc., may be represented graphically by curves composed of a repetition 
of congruent parts at certain intervals. Such a periodic curve may in 
turn be represented analytically by a periodic function of a variable, 
i.e., by a function such that f(« + k) = f(x), where k is the period. Thus 
the functions sin x and cos x have a period 2 z, since sin (x + 27) = sinx 
and cos (x +27) =cosx. Again, the function sin5x has a period 
27/5, since sin 5 (x + 27/5) = sin(5x +27) = sin5.x, but the func- 
tion sin x + sin 5 x has a period 2 7, since sin (x +27) + sin5 (x + 27) 
= sinx + sin 5x. 

Now, any single-valued periodic function can, in general, be expressed 
by an infinite trigonometric series or Fourier’s series of the form 


y = f(x) =a + a,cosx + a,cos2x +--+ +a, cosmx+--- 
+bsinx + dsin2x+---+6,sinnx+---, 


where the coefficients a, and 0, may be determined if the function is 
known. This series has a period 27. But usually the function is un- 
known. Thus, in the problems mentioned above, the curve may either 
be drawn by an oscillograph or by other instruments, or the values of the 
ordinates may be given by means of which the curve may be drawn. 
Our problem then is to represent this curve approximately by a series of 
the above form, containing a finite number of terms, and to find the 
approximate values of the coefficients a, and b;. The following sections 
will give some of the methods employed to determine these coefficients. 

86. The fundamental and the harmonics of a trigonometric series. — 
In Fig. 86a we have drawn the curves y = a, cosx, y = b; sinx, and 
y =a,cosx + bh sin x. 

The maximum height or amplitude of y = a,cosx is a; and the 
period is 27. The amplitude of y = b; sin x is b, and the period is 27. 
Now we may write 


: ae 5 a 
y = a,cos +b sin x = Var Foe | sins + cose, 


by 
Va? + 0; Va? + 6b? 


170 


U 
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and letting Vaz +62 = C1, aoe = COS i, ares = sin q, 
1 1 1 1 


we may write 


y=asin(x+¢), where a= Va2+b2, ¢ = tan z 
Gee BS eh 


Here c, is the amplitude and ¢; is called the phase. The wave rep- 
resented by y = asin (x + ¢1) is called the fundamental wave and 
y = a, cos x, y = b; sin x are called its components. 


Fic. 86a. 
Similarly, we may represent y = a; cos kx, y = by sin kx, b u 
and y = a, cos kx + by sin kx = cz, sin (kx + gx), | | 1 
where ce = Vaz +62 and dx = tan@! az/by. 


The wave represented by y = c& sin (kx + ¢,) is called the kth har- 
monic, its amplitude is c;,, its phase is x, its period is 2 7/k, since 


ta (x 4 2) ie 6.| ees ee orn eee 


and its frequency, or the number of complete waves in the interval 2 7, 
is k. 
The trigonometric series is often written in the form 

y = cotarsin (x+¢1) + sin (2% + ¢2)+ +++ +ensin(nx+dn) +--+, 
showing explicitly the expressions for the fundamental wave and the 
successive harmonics. The more. complex wave represented by this 
expression may be built up by a combination of the waves represented by 
the various harmonics. Fig. 86b shows how the wave for the equation 


y= 2sin(x +2) + sin(2x—27)44sin(3x +52), 


or 


= VJ 2 ; : 2D 
y = cos — ~ 9 cos 2% +— cos 34+ V3sin x— } sin 2x—— “sin 3 x 
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is built up as the combination of the fundamental and the second and 
third harmonics, and how the fundamental wave is modified by the addi- 
tion of the harmonic waves. 


Fic. 86d. 


In the case of alternating currents or voltages, the portion of the wave 
extending from x = 7 to x = 27 is merely a repetition below the x-axis 
of the portion of the wave extending from x = 0 to x = 7; this is illus- 
trated in Fig. 86c where the values of 
the ordinate at x = x, +7 is minus 
the value of the ordinate at x = x,. 
Since 
an X sin (k [x + x] + dy) 
sin (kx + ¢; + kr) 
+sin (kx + ¢;) if k is even 

Fic. 86c. = —sin (kx + @¢x) if Rk is odd, 
the series can contain only the odd harmonics and has the form 
y = co +a sin (x + $1) + cssin (3% + os) +ossin(5* +5) +--+, 
or 
y = +a.cosx +a;cos3x+a;cos5x+:-:-: 
+bisinx + bssin3x+b;sin5x+---. 


les a ae 


“Tread (Ce Re seat vs Thre , 
Vv 


DETERMINATION OF THE CONSTA 
|? lont. 


87. |Determination of the constants when the/function is known. — 
If, in the series 


y = f(x) = a +acosx +a,cos2x+ - + a,cosnx +++. 
+bsinx+bosinax+--+-+b,sinnx+.-- 


we multiply both sides by dx and integrate between the limits 0 and 27, 
we have 


2 1 2 21 2 a 
i: ydx = ay ata | cost de +... an ff cosuxdx + - 
0 0 


Qa 2 1 
+b ff sinxdy +... +n f sinnxdx + ++. 
0 0 
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Qa r 
x sin Xx 
0 


+a 


2a 
tee ek Ht 
0 


Qe 


— by 


COS X 


= 2 7p, since all the other terms vanish. 


If we multiply both sides by cos kx dx and integrate between the 
limits o and 2 7, we have 


2 1 20 2 or 
ycoskx dx = ay f coskxdx + .:--:- + anf cos? kxdx +-+-- 
0 0 
2 1. 214 
tan f cosuxcoskxdx + +>: +b. f sinnxcoskxdx +++. 
0 0 
_ a]. 2a ak sin 2 kx |? 
= sin kx | + Phe x + eli 
4 on sin(n —k)x , sin(n +k) x|? 
2 n—k n+k 0 
bn cos (m — k)x , cos(n +k) x 20 
2 n—k n+k 0 


= 7a, since all the other terms vanish. 


Similarly, if we multiply both sides by sin kx dx and integrate be- 
tween the limits o and 27, we have 


2a 21 2 3 
i, y sin bx dx = a [ sinkxdx+--- +a, f cosuxsinkxdx+.+-- 
0 0 


Qa 2a 
to thf sin? kx dx +--+ +n f sin nx sinkx dx +++. 
0 0 


Qa 

=-_ cos kx ; So a6 oe — nz eel oe 
as cos (k — m)x , cos(k +) x|?* 

wes k—n Rn 0 
b,\sin (7 —k)x sim (a + k) x|2* 

Diy n—k n+k 0 


= mb, since all the other terms vanish. 
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, (> \ | \ 
Collecting our results, we have \ .| ( iA 
I wT I 2 1 it "la F 
ad = — ydx, GQ =- ycoskxdx, bh, =- y sin kx dx, 
Lc ee te BU oe 


where k = 1, 2, 3, Each coefficient may thus be independently 
determined and thus each individual harmonic can be calculated without 
calculating the preceding harmonics. 

88. Determination of the constants when the function is unknown. — 
In our problems the function is unknown, and the periodic curve is drawn 
mechanically or a set of ordinates are given by means of which the curve 
may be approximately drawn. We shall represent the curve by a trig- 
onometric series with a finite number of terms. We divide the interval 


from x = 0 to x = 27 into equal intervals and measure the first n 


ordinates; these are represented by the table 


27 T 67 2 27 
0 =e ae ee lp ee oe ee 
x nN nN nN nN 
Xo M1 Xo 33" ermw ve Xr oietonite Xn=1 
y Yo Sat Y2 aU shia tate See onc pil 


We wish to determine the constants in the equation 


y =a-+acosx+--- +a,coskx+--- 
+obsinx +--+ +&sinkx+.---, 


where the number of terms is 1, so that the corresponding curve will pass 
through the x points given in the table. Substituting the 7 sets of values 
of x and y in this equation, we get 7 linear equations in the a’s and 0’s of 
the form 


Vr = Ayo + a,Ccosx, + - ++ +axzceoskx, + - 
+ by sin x, + ~> +bysinkx,-+---, 


where 7 takes in succession the values 0, 1, 2,...,2-—1. We may 
now solve these equations for the a’s and 0’s. 

We shall first state two theorems in Trigonometry concerning the 
sum of the cosines or sines of m angles which are in arithmetic progression, 
Viz.: 


cos (a + 78) = cosa + cos (a + B) +cos(a+28)+--- 
np 


sin — pee 
> cos (a + = 8). 


sin — 
2 


+ cos (a + [n — 1] 8) = 
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Ssin (a + 78) = sina + sin (a+ 8) +sin(a+28)+--:: 


np 


Si = 


+ sin (a + [n — 1] 8) = 


If we let a = 0 and B = ou these become 


27 
cos rl — 
n 


sin lr ae L(n—1)7 


0, since sin la = 0,° 
sin — 

n 
sin Ir 


: 27 .ln—i1)t : ; 
Ddsin rl — = ——— sin ———— = 0, since sinlz = 0, 
n ale n 
sin — 
n 


* We may prove these theorems as follows: 
By means of the well-known trigonometric identities 


2cosusinyv = sin (u + v) — sin (u —v), 2sin usin v = cos (u — v) — cos (u + ») 

we may write the identities 

2 copa gal = sin(a +8) — sin (« - 8). 2 noes ='Cos (. a) — cos ( + “. 
2 2 2 2 2 2 

2cos (a+8) sin = sin( at 32) —sin (a+ 2, 2sin (w+ 8) sin =cos(a-+-2) —cos{ a+ 32). 

2cos(a+2 p)sin-=sin{ a+ 52) — sin a +58), 2sin(a+2 B)sinf—cos{ a+ 52) —cox{ at 52), 


2 2 


2 cos (a+[n—1] 8) sin = sin (a+ 2143) 2 sin (a+[n—1] B) sin’ = cos (a+2*=39] 
— sin (« + 27=5) . —cos (a+ - _— 6) 
Adding, we get | Adding, we get 


2 sin S'cos (a+ 78) = sin{ a+ 2n— 1) 2 sin dain (+78) =cos (« F 4) 


- sin(a - 8) — cos (« + 2n= | 
= 2008 (a+ *> 1A) sin Ze = 2sin(a + "=" 6) sin %. 
np 


sin ek sin = i 
B cos( a+ m1). LE >in (a+rB) = sin (a+2> 8): 


ales sin = 
sin 5 


* >, cos (a+78) = 


176 EMPIRICAL FORMULAS — PERIODIC CURVES Cuap. VII 


for all values of / except 


1 =0, when >dcos ne —— = Dicoso =H, 
1 =n, when cos r= = De Cos 27r = n. 
° 27 
Since x, = 7 — We may finally state that < 


>)cos lx, = 0, except when ] =o orl=n 
=n, when! = oorl = n. 


sin lx, = O for all values of J. 


To determine a we merely add the equations, and get 
Si, =n t-:-- +a D/cos kx, + .---. + ax Disin kx, + cee 
= nd, since all the other terms vanish. 


To determine a, we a, we multiply each c of the 2 equations by the coefficient 
of a, in that equation, 2.e., by cos kx,, and add the n resulting equations; 
we get 


DE cos kx, = ao dos kx,» + +++ +a, d)cos? Peta 


+a, >)0os px cos kept eee + bp sin px,coskx,++++. | 


Now, 
>,cos RX, = O% 
cos bx Cos kx,” = 4 >dicos (pb +k) x, +34 >d/cos (pb — k)xp= 0; 
sin Px, COS RX,* = 4 sin (p+k) x, +4 sin (p—k) x, = 0; 


Dcos? kx, os y; (1 + cos 2 kx,) = = +: >,c0s2 ke, = : if es 


Nis ioe 


* We use the trigonometric identities ns 
2 cosucosv = cos (u + v) + cos (u — 2). 
2 sin ucosv = sin (u + v) +sin (wu — 0). 
2 sin usinv = cos (vu — v) — cos (u +). 
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n n 
Hence, > cos kx, = 3 ve except when k = ; 
n 
= nay, when k = x 


To determine 5; we multiply e each of the 7 equations by the coefficient 


of b, in that equation, i.e., by sin kx,, and add the » resulting equations; 
we get 


>»): sin kx, = ad sin Rip + - ++ dy >)cos pxy sin Rx, + +66 


a by, Dysin? kx, +--+ +0, sin px, sin kx, f+ eee, 
Now, 
Ddsin Rx, = 0; 
cos pXxr sin Rey dsin (R+ p)x%+ 3 sin (R —, p) x, = 03 
dsin px, sin kx,* = 4 Sicos (pb —k) x, — 3 >/cos (bP +k) x, = 0; 


> sin? kx, = pe — cos 2 kx;) = 2 — 5 dicos 2 ke, ==, ie 


: n 
Hence, > Sih kX, = 5 bx. 


Collecting our results, we have HgUe 


= rae = = (% ea A * tant) Pak 
0g 
I ¥ 
as a Dy. cos 7 a= 7 Yor cos re zee Wb ts Is 


a = — ret cos sate == = (0 cos kx» a0 J1 COs Rey +++ + FYn-1 C08 RXn-1), 


bh = = >: Sinki, = 2 (yo sin kxp + yi sin kx, + + + + + yn-1 Sin Rxp -.). 


* We use the trigonometric identities (se 
2 cos ucosv = cos (u +7) + cos (u — 2). 
2 sin ucosv = sin (u +v) + sin (uw — »). 
2 sin u sin v = cos (u — v) — cos (u +2). 
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If 2 is an even integer, our periodic curve is now represented by the 
equation 


n 
y=a+acosx+---> +a,coskx+:-:: hn COS 5 % 
2 


4+hsinx +--+ +bpsinkx+--- + by_sin(S 1) x. 
2 


The x coefficients are determined as above. Thus — 

a is the average value of the ” ordinates. 

a, is the average value of the 7 ordinates taken alternately plus and 

2 
minus. 

a, or b;, is twice the average value of the products formed by multiply- 
ing each ordinate by the cosine or sine of k times the corresponding value 
ON 

We note that each coefficient is determined independently of all the 
others. 

If we wished to represent the periodic curve by a Fourier’s series con- 
taining ” terms, but had measured m ordinates, where m > n, we should 
have to determine the coefficients by the method of least squares. The 
values of the ordinates as computed from this series will agree as closely 
as possible with the values of the measured ordinates. It may be shown 
that the expressions for the coefficients obtained by the method of least 
squares have the same form as those derived above.t 


* We may also derive the values of the coefficients as follows: In Art. 87, we have 
shown that 


27 Qa 
{f y cos kx dx = ay A cos? kx dx, 
0 0 


since all the other terms vanish. 


If we replace the integrals by sums, and take for dx the interval 2 7/n, this becomes 
n : n 
>» cos kx, = ay > /cos" Rx, = 3 tks if k #0 or k# = 
= naz, if k =o or k =<. 


n n 
Hence, >» = Nd, >” COS 5 Xr = Nn, >» cos kx = = ah. 
2 


Similarly we may show that >> SIlietn = - br. 


+ See A Course in Fourier’s Analysis and Periodogram Analysis by G. A. Carse and 
G. Shearer. 
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We shall illustrate the use of the above formulas for the coefficients by 
finding the fifth harmonic in the equation of the periodic curve passing 
through the 12 points given by the following data (Fig. 89). 


z y cos 5 z sin5z y cos 52 ysin dz 
OF 9.3 1.000 0.000 9.30 0.00 
30° 15.0 —0.866 0.500 —12.99 7.50 
60° 17.4 0.500 —0.866 8.70 —15.07 
90° 23.0 0.000 1.000 0.00 23.00 
20 BY MnO) —0.500 —0.866 —18.50 — 32.04 
150° 31.0 0.866 0.500 26.85 15.50 
180° 1573 —1.000 0.000 —15.30 0.00 
210° 4.0 0.866 —0.500 3.46 — 2.00 
240° — 8.0 —0.500 0.866 4.00 — 6.93 
2702 —13.2 0.000 —1.000 0.00 13.20 
300° —14.2 0.500 0.866 — 7.10 —12.30 
330° — 6.0 —0.866 —0.500 5.20 3.00 
z= 3.62 — 6.14 


ds = ys > COs 5 X; = 0.60; bs = 5 > sin 5 X,= —1.02. 


Hence the fifth harmonic is 0.60 cos 5 x — 1.02 sin 5 x. 

It is evident that the labor involved in the direct determination of the 
coefficients by the above formulas is very great. This labor may be 
reduced to a minimum by arranging the work in tabular form. These 
forms follow the methods devised by Runge * for periodic curves in- 
volving both even and odd harmonics (Art. 89), and by S. P. Thompson f 
for periodic curves involving only odd harmonics (Art. go). 

89. Numerical evaluation of the coefficients. Even and odd har- 
monics. — 

(1) Six-ordinate scheme. — Given the curve and wishing to determine 
the first three harmonics, 1.e., the 6 coefficients in the equation 


y = a +a,cosx + a,cos2x +a3;cos3x+;sin x + & sin 2 x, 
we divide the period from x = 0° to x = 360° ft into 6 equal parts and 


* Zeit. f. Math. u. Phys., xlviii. 443 (1903), lii. 117 (1905); Erlauterung des Rech- 
nungsformulars, u.s.w., Braunschweig, 1913. 
1 Proc. Phys. Soc., xix. 443, 1905; The Electrician, 5th May, 1905. 
t If the period is taken equal to 2 7/m instead of 2 z, the representing trigonometric 
series has the form 
y =a +a cosm) + a2.cos2md@+:--:> 
+bhsinmé+bsin2m@++--, 
where @ represents abscissas. By the substitution m@ = x or 0 = x/m, the series be- 
comes 
y=%+acosx+a@cos2x+::- 
+bhsinx + bsina2x+---, 
and this has a period 27. The abscissas from 6 = 0 to 6 = 27/m now become the 
abscissas from x = 0 to x = 27, and we proceed to determine the coefficients in the 
second series as outlined. Having determined the coefficients, we finally replace x by 
mé. 
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measure the ordinates at the beginning of each interval; let these be rep- 
resented by the following table: 


G | 0° | 60° 120° 180° 240° 300° 
y Yo Aa ye Ya M4 5 
Here ” = 6, and using the formulas on p. 177, we have 
6 a = Vo +1 +4 +43 +4 +45 
643 = Yo —" +2 a Tr =e 


3a, =yocoso° +y,cos 60°+72cos120°+73cos 180° -+44,c0s 240°+y5 cos 300° 
3 dz =yocoso° +y1cos 120° +y2c0s 240° +y3cos 360° +44c0S 480° +y5scos 600° 
3b, =yosino° +y:sin 60° +y2sin 120° +y3sin 180° +y4sin 240° +y;sin 300° 
3 by = yosino® +y,sin 120° +y2sin 240° +y3sin 360° +y4sin 480° +5sin 600° 


We arrange the y’s in two rows, 
Yo WM Ap 43 


5 V4 
Sum Vo Vy V2 U3 
Diff. Wy We 


where the v’s are the sums and the w’s are the differences of the quantities 
standing in the same vertical column; thus, vp = yo, 11 = ¥1 + Ys, M1 = 1 
— ys, etc. Since cos 240° = cos 120°, cos 300° = cos 60°, etc. We may 
now write 

6% =H tn + U2 + 03 

6 a3 = HW — YH + vp, — Us 

31 = % + cos 60° + v% cos 120° + v3 cos 180° 

3 de = % +1 Cos 120° + vw cos 240° + v3 cos 360° 

2:0, = W, sin 60° + we sin 120° 

3b, = W Sin 120° + we sin 240° 


We arrange the v’s and w’s in rows, 


Uo V3" Wy 
U3 V2 We 
Sum Po pi tal 
Diff. do qa Si 
and we now write ; 
6% =ht+ fi, 6 ds = qo — i, 
3UN=Qt+FH, 342 = fy — $ fr, 
V2 
gh = —2n, she oa 


Example. Determine the first three harmonics for the following data 
taken from Fig. 860. 
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x 0° | 60° 120° 180° 240° 300° 
y | 0.47 ey, | 2.20 | —2,20 | = 1,64 —0.49 
0.47 177 2.20 +=2.20 
= O74 =A 
v 0.47 1.25 0.56 =A 
w 2720 ey coy 
0.47 1228 2.26 
—2.20 0.56 3.84 
p ry Wari 1.84 a 6.10 
q 2.67 Or72 Ss —=1,59 
6a =0.II, 6 a3 = 1.95, 3a, = 3.03, 
Zit, = — 2.65, 3b; =-5,28, 3 be = —1.37. 
Hence, dy = 0.02). dy = 1.01,. =de2=-—0.88, dy = 0.33, 
6b) = 1.76, 6b: = —0.46, 


and y = 0.02 + 1.01 cos x —0.88 cos 2 x + 0.33 cos 3 x 
+ 1.705in « —0.46 sin 2 x, 

The equation from which the curve in Fig. 865 was plotted was 

y = 2 sin (x = 4 + sin (2. == +<sin (3x +37) 
=cos x —0.87 cos 2 x +0.35 cos 3 x + 1.73 sin x +0.50 sin 2 x—0.35 sin3 x. 

We observe the close agreement between the two sets of coefficients, 
the small discrepancies being due to the approximate measurements of 
the ordinates for our example. 

(Il) Twelve-ordinate scheme.—Given the curve and wishing to deter- 
mine the first six harmonics, 2.e., the 12 coefficients in the equation 
W= Wy. G4 COS X + 02 C08 2'% + 03 COS 3X. -F d4c0Ss 4% + a; cos 5 X 

+a,cos6x + b:sinx + bsin2x+63;sin3x+hsin4 x +6;sin 5x, 
we divide the interval from x = 0 to x = 360° into 12 equal parts and 


measure the ordinates at the beginning of each interval; let these be rep- 
resented by the following table: 


Ba 30° | 60° | 90° | 120° | 150° ] 180° | 210° | 240° | 270° | 300° 
SON Yi te Var Yet <4 Ys 6 Mi ys Yo Yio Yu 
Here ” = 12, and the formulas for the coefficients give 
I2 dG = Yo as + Yo ree tev it 
I2 ds = Wo a4 + ye or Sie Vistas 
6 a; = yo cos 0° + 1 cos 30° + y2cos 60° + +--+ + yu Cos 330° 
6 a2 = ” cos 0° - 4, cos 60° + yecos 120° + +++ + Hi cos 660° 
aie reac sane eo? + .-.- ce I 


6 by = yo sin 0° + y; sin 60° + yesin 120° + +++ + yu sin 660° 
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If we arrange the y’s in two rows 


Yo ah y2 ys y4 Ys M6 
Yiu 10 Of) Ys 7 
Sum Uo Vy V2 ¥ U4 U5 U6 
Diff. Wi We W3 W4 Ws 
and remember that cos 330° = cos 30°, sin 330° = —sin 30°, etc., the 
above equations may be written 
I2@ =% +n + Ue afox Penta 
12 dg = % — + ve Sr a A 
6 a, = +, cos 30° + v cos 60° + -- + + v5 c0s 180° 
6 dz = % + Cos 60° + v cos 120° + - - - + u5 cos 360° 
6b; = Ww, sin 30° + w,sin60° + --- + ws; sin 150° 
6 by = Ww, sin 60° + wesin 120° + + - - + w5 sin 300° 


If we now arrange the v’s and w’s in two rows 


Vo Vy Ve V3 Wi We W3 
U6 Us U4 Ws Ws 
Sum Po Pi po ps "1 ro r3 
Diff. qo qu Qe Si S2 
the equations may be written 
I2m=OtH + @ + g 
I2 a = po — fi + hr “Ps 


6 a1 = go +H cos 30° + g cos 60° 
6 dz = po + picos 60° + pf. cos 120° + 3 cos 180° 


6 db; = “sin 30° + r.sin 60° + 73 sin 90° 
6 by = 5; sin 60° + sq sin 120° 


Finally, if we arrange the p’s, q’s, and 7’s as follows: 


Po Pi ry qo 
peo ps r3 q2 
Sum i h Diff. ty te 
the equations become 
I2Q@ = Ith. I2d, = lo = hy. 


6 a, = go +qisin 60° + gsin 30°. 6 ds = go — gisin 60° + gq sin 30°. 
6 a2 = (pbo— bs) + (fi-— pe) sin 30°. 6 ay = (fo + ps) — (pi + po) sin 30°. 
OO = iy (Oy = he 

6 db; = r, sin 30° + rosin 60° + 7s. 6 bs = 1, sin 30° — 72 sin 60° + 75. 
6 be = (1 + 52) sin 60°. 6 b; = (s1 — se) sin 60°. 
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We may now arrange the above scheme in a computing form as fol- 
lows: 


Ordinates Yo Mi Ye ys 4 Ys ¥6 
Yu Vi0 Yo ¥8 Yr 
Sum Vv) V} V2 V3 V4 Us U6 
Diff Wj We W3 W4 Ws 
Vo Vy V2 U3 W We W3 
U6 Us U4 Ws W4 
Sum Po Pp po ps Lal Yo 3 
Diff. do a Qe S1 So 
Po Pi ial qo 
pr bs "3 ge 
Sunt @ ibs. he Diff. ty ty 


Multipliers of the quan- 
tities in the same 


ener naca letore Cosine terms Sine terms 
these are entered 
sin 30° = 0.5 qe —p2 D1 rT) 
sin 60° = 0.866 Nn To | Si Se 
sin 90° = 1.0 qo Pol sa Ps te tee hts ceed CARE 
SUIMOMISUCOLUIEIM me nee ne all a Mllaecey ateee || pratiamaen |e Sareea. 
Soran eli Dalered kouanint* || eeu 8 haat eae ees el ae | Ee] RMA eer etal Bane eet al errr 
Sum 6 a 6 de 6 a3 12 ao 6 by 6 be 6 bs 
Difference 6 as 6 as 12 as 6 bs 6 bg 
| 
Checks: Yo = Ay +a, + de + a3 + a4 + as + ae. 


(bd: + bs) + V3 (be + bs) + 2 bs. 


Vie— Vin 


Result: y = a +a,cosx +a,cos2x +--+ + agcos6x 
+b,sinx + bsin2x+--+ +;sin5x. 


150° 180° } 3 360° 


120° 


Fic. 89. 
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Example. In the periodic curve of Fig. 89, the interval from x = 0° 


to x = 360° is divided into 12 equal parts and the ordinates y) to yu are 
measured. 

x | 0° | 30°| 60°| 90° | 120° | 150° | 180° | 210° | 240° | 270° | 300° | 330° 
y 19.3 |15.0|17.4}23.0| 37-0| 31.0] 15.3 | 4.0 | —8.0|—13.2|—14.2)— 6.0 
We shall determine the first six harmonics by the above scheme. 
Ordinates 9.3 15.0 17.42 | 23:07) 376 | Git SG 

—6.0° =14.2,-—132 7 =80 4.0 
Sum (v) 9.3 9.0 3.2 9:8- 29,05 “935.0 5:4 
Diff. (w) 31.0. + 3Y,6, 0) 26.2 245,08 927.0 
9.3 9.0 ee Rano ea 21.0 2126! 36:2 
1S2 - 35.0.- © 20.0 27.0 - 45.0 
Sum (p) 2469) APG. “22772 -.0.5 (7) 480 766 4a2 
Diff. (¢q) —6.0 —26.0 —25.8 (s) -—6.0 —13.4 
24:0. § “44:0 43:1 -—6:0 
32.2 9.8 46:2 = 258 
Sum (/) 56.8 53.8 Dit) 11.9 19.8 
Multipliers Cosine terms Sine terms 
0.5 —12.9 —16.1 22.0 24.0 
0.866 - : 66.3 | —5.2 —11.6 
WD) —6.0 Pasty ts} 19.8 56.8 53.8 | 36.2 9 
Sum of Ist col. —18.9 8.5 56.8 60.2 — 5.2 
Sum of 2d col. —22.5 12.2 53.8 66.3 —11:6 
Sum —41.4=6 a, 20.7=6 ag 19.8 | 110.6=12 ag 126.5=6 b, —16.8=6 by 11.9 
Diff. + 3.6=6 as —3.7=6 a =6 a; 3.0=12 ag —6.1=6 b; 6.4=6 by | =6b; 


a= — 6.90, €2=3.45, €3=3.30, @o=9.22, b}=21.08, bo = —2.80,03;=1.98, 
@;=0.60, d= —0.62, Q5=0.25, bs = —1.02, bs =1.07. 
Check: 9.3 = 9.22 — 6.90 + 3.45 + 3.30 — 0.62 + 0.60 + 0.25 = 9.30. 
21.0 = (21.08 — 1.02) + 1.732 (—2.80 + 1.07) + 2 (1.98) =21.02. 
Result: * 


y = 9.22 — 6.90 cos x + 3.45 cos 2 * + 3.30 cos 3 x — 0.62 cos 4x 
' + 0.60 cos 5 x + 0.25 cos 6 x + 21.08 sin x — 2.80 sin 2x 
+ 1.98 sin 3 x + 1.07 sin 4 x — 1.02 sin 5 x, 
or 
y = 9.22 + 22.18 sin (x — 18.12°) — 4.44 sin-(2 x — 50.93°) 


+ 3.85 sin (3 x + 59.04°) + 1.24 sin (4 x — 30.09°) 
— 1.18 sin (5 x — 30.47°) — 0.25 sin (6 x — 90°). 


* The coefficients of the fifth harmonic agree with those found by the direct process 
in Art. 88. The time and labor spent in the computation of all six harmonics by means 
of the above computing form is much less than that spent in the determination of the 
fifth harmonic alone by the direct process in Art. 88. 
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The last result was obtained by using the relations 


a, cos kx + by sin kx = cy sin (Rx + ox), 


where ch = Vae + be and o, = tan 5. 
F k 


(III) Twenty-four-ordinate scheme.— Given the curve and wishing to 
find the first 12 harmonics, i.e., the 24 coefficients in the equation 


y= % +a,cosx +a,.cos2x +--+ + a2.cosI2x 
+bsinx + dsina2x+---+ +b), sin 11 x, 


we divide the interval from « = 0° to x = 360° into 24 equal parts and 


measure the ordinates at the beginning of each interval; let these be repre- 
sented by the following table: 


x ¥ es ae | 45° can 330° | 345° 
y cn Wt meee: ys Yo2 23 


If we use the same method as that employed in deriving the 12-ordi- 


nate scheme, we shall arrive at the following 24-ordinate computing form. 
This form is self-explanatory. 


O 
Yo 


Ordinates Yo V1 Yo V1 Viz 
Vo3 Yo2 eaestes Yi3 
Sum Up V1 V2 Vn Vi2 
Diff. Wy We Wy 
Vo V1 ona 6 U5 U6 W We Cy tx 28 Ws We 
Vi2 Vi sete U7 Wi W0 a He W7 
Sum fp Pi sg 6 Ps De ‘al 9 ee 5 16 
Diff. do nN S44 5 Sj S2 sient s 55 
Po pi pr Ps Sy S2 53 
Ps Ps ps S5 S4 
Sum 1p kh ly ls ky ko kg 
Diff. Mo My, Me NM Ne 
Io h do — 4 = bo n+ 73 — % = uy 
l ls aA-B@—-Gd=h 12 — 16 = U2 
£0 21 
Multipliers Cosine terms Sine terms 
sin 30° = 0.5 me —-lk kk ky 
sin 60° = 0.866 mm, ke | m 38M 
sin 90° =1.0 | go gi | Mo lo —l3 | mo me | ks ky ks 
CU uaa ty Gye JUCRHCOELO) Ladd oo ans Sie eee ND Si sesearrasene IES, SOUR an boat reeves n Ulin neni ROG Ame eotsieata oil [ener tig 
Sh rene Oe te bier lS locate at tad Es cae sees ey Fae estovec ores a ec ereeSee a m[eicacs aecaa lutea ers ur are al Pe eee 
Sum 24 ao 12 a, 12 a 12 be 12 by 
Difference 24 aie 12 ayo 12 ag 12 as 12 dio 12 bs 12 bs 
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Multipliers Cosine terms Sine terms 
sin 15° = 0.259 ds qi | 71 15 
sin 80° = 0.5 q4 q4 T2 T2 
sin 45° = 0.707 q3 ty = 73 eis Uy es 
sin 60° = 0.866 qQ2 —q2 Pe aad 
sin 75° = 0.966 Qn qs | 75 1 
sin 90° = 1.0 qo to qo T6 U2 re 
Se bhcsWce) al Cc co) ee | rs en ee Pe be on, allt denen |) woo Sboe +2 
Guam: of Bd ol ap eee ae, I ae eee ca! este eres re ene eee eee | eee 
Sum 12 ay 12 Q3 12 Os 12 by 12 b3 12 bs 
Difference 12 ay 12 dy 12 a4, 12 by 12 by 12 b; 
Checks: Yo = Ao + ay + a2 +- neces 4 a2. 
4 (yy, — 03) = 0.259 (01 + bu) + 4 (b2 + bio) + 0.707 (b3 + bo) 
+ 0.866 (b, + bs) + 0.966 (ds + bz) + be. 
Result: 
y = + a,cosx +a,cos2x+... +a2cosI2% 


+bsinx +bosina2x+... +dysin IIx, 
or y = @+aqsin(x+ ¢1) +esin (2x+¢2) + +--+ +cpsin (12% + dy). 


We shall now pass on to the evaluation of the coefficients when only 
the odd harmonics are present.* 


“___90. Numerical evaluation of the coefficients. Odd harmonics only. — 


Most problems in alternating currents and voltages present waves where 
the second half-period is merely a repetition below the axis of the first 
half-period; the axis or zero line is chosen midway between the highest 
and lowest points of the wave (Fig. 86c).. We have shown in Art. 86 that, 
in such cases, the trigonometric series contains only the odd harmonics. 
Furthermore, since the sum of the ordinates over the entire period is 


: I : 3 
evidently zero, then a = a > = 0, and the series does not contain the 


constant term da. Again, since 


cos k (x +71) = cos (kx + kr) = — coskx, when k is odd, 
sin k (x + 7) = sin (kx + kr) = — sin kx, when k is odd, 


and Vorr = —Var “. Varx COS B (x + ar): = y, Cos Rx, 
and >» cos kx has the same value over the second half-period as over the 


* T. R. Running, Empirical Formulas, p. 74, gives similar schemes with 8, 10, 16, 
and 20 ordinates, for waves having even and odd harmonics. H. O. Taylor, in the 
Physical Review, N. S., Vol. VI (1915), p. 303, gives a somewhat different scheme with 
24 ordinates for waves having even and odd harmonics. A very convenient computing 
form for the above scheme with 24 ordinates has been devised by E. T. Whittaker for 
use in his mathematical laboratory at the University of Edinburgh; see Carse and 
Shearer, ibid., p. 22, 
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first half. Hence in finding the coefficients we need merely carry the 
summation over the first half-period; thus, 


2 2 ; 
= =e coskx, h= ms dy sin kx, 


where & is odd, x and y are measured in the first half-period only, and n 
is the number of intervals into which the half-period is divided. 

(I) Odd harmonics up to the fifth.—Given the curve and wishing to 
determine the coefficients in the equation 


y =a, cos x +a; cos 3 x +a; cos 5 x +b; sin x +3 sin 3 x+0; sin 5 x, 
we choose the origin where the wave crosses the axis, so that when xp = 0, 
yo = 0, divide the half-period into 6 equal parts, and measure the 5 ordi- 
nates 1, Ye, 3, V4, Vs. Thus we have 
x 20. 60° 90 120 
sy yA y2 Me V4 
For the coefficients we have the following equations: 
3a, = y, cos 30° + y2cos60° + y3cosgo° + y4cos 120° + ys cos 150°. 
3 a3 = y1cosgo° + y2cos 180° + y3 cos 270° + y4cos 360° + ys Cos 450°. 
3 ds = y1 COS 150° + ye cos 300° + 3 cos 450° + 4 cos 600° + ys COs 750°. 
3b; = y1 sin 30° + yesin60° + yssingo° + y4sin 120° + ys sin 150°. 
3 bs = y:singo° + y2sin 180° + y3 sin 270° + y4sin 360° + ys sin 450°. 
3 bs = y: sin 150° + yosin 300° + y3 sin 450° + ys sin 600° + ys; sin 750°. 
Simplifying and replacing the trigonometric functions by their values 
in terms of sin 30° and sin 60°, we may write 


3 a1 = (ve — ya) sin 30° + (1 — ys) sin 60°. 
3 d3 = —(ye — 4) sin 90°. 
3.425 = (ye — a) sin 30° = Gn aa Ys) sin 60°. 
3 b1 = (1 + 95) sin 30° + (v2 + 44) sin 60° + yz sin go°. 
3 b3 = (91 — Ys + Ys) sin gO°. 
3 bs = (y1 + ys) sin 30° — (ye + ys) sin 60° + y3 sin 90°. 
We may conveniently arrange the work in the following computing 
form: 


° ° 


150° 
5 


Multipliers Cosine terms Sine terms 
1 Ae 3 
Se sin 30°=0.5 | de 
Sum St Se $3 sin 60° = 0.866 dy 
: sin 90° = 1.0 
Diff. a ds Sum of Ist col. | ...... 
Checks: 0 = a + a3 + Gs. Sum of 2d col. | ...... 
pie ae Sum 38a 
SE aes ae Diff. 3 as 


Result: 
y = aicosx +a3;cos3x +a;cos5x+bisinx + b3sin3 x + bs sin 5 x. 
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The following example will illustrate the rapidity with which the co- 
efficients may be determined. 
- Example. We wish to analyze the symmetric wave of Fig. 90a, 7.., 
to find the coefficients of the 1st, 3d, and 5th harmonics. Choose the 
x-axis midway between the highest and lowest points of the wave, and 


0° 30° 60° 90° 120° 150° 18 


FIG. 90a. 


the origin at the point where the wave crosses this axis in the positive 
direction. Then divide the half-period into 6 equal parts and measure 


the ordinates y1, ..., 5. These are given in the following table: 
“ 30° 60° 90° 120° 150° 
y 10:7 2.8 20.5 26.5 16.6 


We arrange the work in the above computing form. 


Multipliers Cosine terms Sine terms 


10.7 2.8 20.5 25.37 
See .7120.5 27.3 20.5 
ee SUNLOMES COMM ieee 34.15 PALS 
Sum (s) 27.3 29.3 20.5 |Sum of 2d col.| —5.11 25.37 20.5 
Diff. (d) —5.9 —23. Sum —16.96 23.7| 59.52 
ai Diff. —6.74 8.78 6.8 
Divide by 3 jai=—5.65] as=7.9b;=19.84 y, 
as= —2.25 bs= 2.93| bs = 2.27 


Check: a + as + O& ='—5.65 + 7.00 — 2.25 = 0) 
b, — bs + bs = 19.84 — 2.27 + 2.93 = 20.5 = Ys. 
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Result: 


y = —5.65 cos x + 7.90 cos 3. x — 2.25cos 5x 
+ 19.84 sin x + 2.27 sin3 x + 2.93 sin5 x. 
(II) Odd harmonics up to the eleventh.— Given a symmetric curve and 
wishing to determine the coefficients in the equation 


y = a,cosx + a3;cos3x+ +--+ +4a@yCOosIIx 
+bsinx + b3sin3x+ --- + dy,sin11x, 


we choose the origin at the point where the wave crosses the axis, so that 
yo =_0, divide the half-period into 12 equal parts, and measure the 11 


ordinates yi, yon, ..., Yu. Thus we have 
x | 5 | 30° | 45° | 165° 
» y1 Ye ¥3 yu 
For the coefficients we have the following equations: 
6 a = yi cos 15° + yo cos 30° + + + + + yy Cos 165°. 
6 a3 = yi Cos 45° + ye cos 90° + + + + + yy COS 495°. 
6b; = yr sin 15° + yw sin 30° + - ++ + yy sin 165°. 
6 b3 — 


yi sin 45° + ye sin go° + --- + yu sin 495°. ° 


If we arrange the ordinates in two rows, 


M1 a) 3 M4 V5 M6 
Vu Yi0 9 Ys Vi 
Sum S1 Se 53 S4 S5 S6 


Diff. dy dy ds; dy ds dg 


replace the trigonometric functions by their values in terms of the sines 
of 15°, 30°, 45°, 60°, 75°, 90°, and collect terms, we may write 


6a = dssin15°+d,sin 30°+d; sin 45° +d» sin 60°-+d, sin 75°. 
6 a, = —ds sin 15°-+d, sin 30° —ds sin 45° +d, sin 60° —d, sin 75°. 
6a; = dsin15°+d,4sin 30°—d; sin 45°—d, sin 60°+d; sin 75°. 
6 a7 = —d sin 15° +d, sin 30°+d; sin 45° — dz sin 60° —d; sin 75°. 
6b, = 5, sin 15°+5» sin 30°-+53 sin 45°+5, sin 60°+55 sin 75°+5¢ sin 90°. 


6bu= 51 SiN 15°—s2 sin 30°+s3 sin 45°— 54 sin 60°-+55 sin 75° — sgsin 90°. 
6b, = Ss sin 15°+S2 sin 30°—53 sin 45°— 5, sin 60°+5; sin 75°-+5¢ sin 90°. 
6b; = 55 Sin 15° —S2 sin 30°—S3 sin 45°-+s4 sin 60°+5) sin 75°—s¢sin 90°. 
(Qu (dy — dz — ds) sin 45° —d,sin 90°. 
6 dy = — (dy — d3 — ds) sin Hee —d sin 90°. 
6 b3 = (Spchkiss <= s5) Sill 45- + (s2 — 5) sin go°. 


6 by = (si: + Ss — Ss) sin 45° — (s2 — Se) sin go°. 
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We may conveniently arrange the work in the following computing 


form: 
Ni “Ne-Yo Ge Ving (Nomar ve Sf 58 She 
Vite 10) ey 9) 7 S20 56 = 1h 
Sut 4.555 SiSe asa) Sse ss d, — d3 — ds = & 


Difi Wd," des “dy adr) ds 


Multipliers Cosine terms Sine terms 
gin 15. = "01259 ds dy $1 S5 
sin OS = ass da d4 Se S82 
sin 45° = 0.707 ds e€y —ds S83 . T1 —=S3 
sin 60° = 0.866 dy —dz S4 — 4 
sin 75° = 0.966 ; dy ds | Ss S81 
sin 90° = 1.0 —ds S6 T2 Se 


Sum pore lstieoleme | eee ae ee eee ee 
Sumivol 2dsCo] ess hee eel eee ee eg | eee 


Sum 6 a 6 a3 6 ds 6 by 6 b3 6 bs 
Diff. 6 a1 6 ag 6 az 6 bu 6 bo 6 b7 
Checks: a, +43 + a5 + a; +a + au = 0, 
b; — bs + bs — b7 + dy — bu = Ye. 
Result: yy = a,cosx + a;cos3x+ --- +a coslix 
+b,;sinx+d3sin3x+--- + dysinIIx. 


90° 105° 120° 186° 150° 165° 


45-60 75" 


Fic. 90). 


Example. Fig. 90b represents a half-period of an e.m.f. wave whose 
frequency is 60 cycles. We wish to find the odd harmonics up to the 11th 
order. Choose the x-axis midway between the highest and lowest points 
of the complete wave and the origin at the point where the wave crosses 
the x-axis in the positive direction. Divide the half-period into 12 equal 
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parts and measure the ordinates y, ye... , Yu. These are given in 
the following table: 
Diels ik ml 45 4160" 675°] 90711105... 120.1435 150° | 165° 
ees! BIO a 272 e200 2335 40k 358") 50 “1 30 33 
We arrange the work in the above computing form. 
4 21 19 27 29 33 SRO 70 not 
33 30 50 38 46 Sle —033 =18 = 7, 
Sums (s) 37 51 69 65 75 33 —29+31 +17 =19=4 


Diff. (d) —29— 9-—31—II-—I17 


Multipliers Cosine terms Sine terms 


0.259 —4.4 —7.5| 9.6 9.4 
0.5 —5.5 —5d. 25.5 25.5 
0.707 —21.9 13.4 21.9] 48.8 lO —48.8 
0.866 —7.8 7.8 56.3 —56.3 
0.966 —28.0 —16.4! 72.5 35,7 
nO 11.0 33.0 18.0 33.0 
Sum Ist col. —13.3 11.0 2.3 130.9 21.9 6.3 
Sum 2dcol.| —54.3 13.4 —2.0 114.8 18.0 DAD: 
Sum —67.6 24.4 0.3 245.7 39.9 8.5 
Diff. 41.0 —2.4 | 43 16.1 3.9 4.1 
Divide by 6 |a1 =—11.27| az= 4.07] a5=0.05 | bi =40.95 | bs3=6.65| b;=1.42 
Qu= 6.83 | a= —0.40| a7=0.72 | bu= 2.68 | bo>=0.65] b7=0.68 
Check: 
Qitdst+ -- + +ay= —I11.27 +4.07 + 0.05 + 0.72 — 0.40 + 6.83 =0, 
bi—b3+ +--+ —bdn = 40.95 —6.65+1.42 —0.68 +0.65 — 2.68 = 33.01 = ye. 
Result: 
y = —11.27 COS x=} 4.07 Cos 3. x. + 0.05 CoS 5% + 0.72: Cos 7 X — 0.40 C089 x 


+ 6.83 cosII x + 40.95 sin x + 6.65 sin3 x + 1.42sin5x 
+ 0.68 sin 7x + 0.65 sing x + 2.68 sin II x. 
(III) Odd harmonics up to the seventeenth. — Given a symmetric curve 
and wishing to determine the coefficients in the equation 
y =a, cosx + @3cos3x% +--+ +> + a1;COS1I7X 
+bsinx + b3:singx +--+ +0sin 17x, 
we choose the origin at the point where the wave crosses the axis, so that 
yo = 0, divide the half-period into 18 equal parts, and measure the 17 


ordinates yi, ye, . . - , M177 Thus we have 
t | 10° | 20° | BO% | ae | 70% 
yy Ml y2 Ys M17 


If we use the same method as that employed in deriving the 11-ordi- 
nate scheme, we shall arrive at the following 17-ordinate computing form. 
This form is self-explanatory. * 
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V1 Ve V3 V4 V5 Veo V7 Ys Yo 

iz Vi6 15 yu 413 12 Yi 10; 
Sum S1 So $3 $4 Ss S6 S7 Sg Sg 
Diff. dy ds dz ds ds de dz ds e 


51 52 53 nr de d, —ds =e 
Ss $4. Sq —173 —-d, —ds es 
Si; = Ss : —dzg —d, 
Sum ry re 13 4 ey eo e3 4 
Multipliers Cosine terms Sine terms 
sin 10° = 0.1737} ds —ds d4 $1 87 —S5 
sin 20° = 0.3420 dy —ds dy S82 —S, — 8, 
sin 30° = 0.5000} de 1 de de 83 1 $3 —S83 
sin 40° = 0.6428 ds dy —d; S4 83 So 
sin 50° = 0.7660} d, dg —dy S85 Si 87 
sin 60° = 0.8660 ds 2 —d3 —dsz 86 f2 — S65 S6 
sin 70° = 0.9397) de —d, —d3 87 —S85 81 
sin 80° = 0.9848 dy dy ds 83 S82 — SZ 
sin 90° = 1.0000 3 €4 | So 3 So —89 4 
AS LDEO TU 0) chs Ko) ran C10) Peal ee a Ieee A mtn ome liaaar me Memtotiall Son esas ob : 
Swim Gebel, Jecasculeanses es ce td ee etc Prey oor sd Be 
Sum 9 ay 9 a3 9 a5 9 az 9 ag 9 b, 9 b3 9 bs 9 bz 9 bg 
Diff. 9 Qi7 9 15 9 a3 9 ai 9 biz 9 bis 9 bis 9 bi 
Check: aqat+atat::: +a,= O, 
iE i ei Soe ye Tv 
Result: y=a,cosx +a3;cos3x +--+ +a17cos17x% 
+obsinx + b:sin3x+--- + dy;sin IIx. 


Similar computing forms may be constructed for symmetrical waves 
‘containing odd harmonics up to the seventh, ninth, etc., orders. 

o1. Numerical evaluation of the coefficients. Averaging selected 
ordinates.* — We are to determine the coefficients in the trigonometric 
series 

y= +acosx +a,cos2x+--- +a,coskx +--+-- 
+hsinx +besinax+--- +bhsinkx +---. 


Let a, and 6, represent the coefficients of any harmonic. We divide 
the period 2 7 into 7 equal intervals of width 2 +/n and measure the ordi- 
nates at the beginning of these intervals. We have the table 

x | Xo vy | X2 Xn-1 


ee 
sy) Yo v1 2 Yr Vn-1 


* These methods have been developed by J. Fischer-Hinnen, Elekrotechnische 
Zeitschrift, May 9, 1901, and S. P. Thompson, Proc. of the Phys. Soc. of London, 
Vol. XXIII, 1911, p. 334. See, also, a description of the Fischer-Hinnen method by 
P. M. Lincoln, The Electric Journal, Vol. 5, 1908, p. 386. 
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Substituting these pairs of values in our series, we have n equations 
of the form 


Vr = Ayo +a, COS xX, + a2cos2x%,+ +++ +a,coskx, +++ 
+ b,sinx, +besin2x,+ +--+ +b,sinkx,+-+-, 
where 7 takes in succession the values 0, I, 2, 3,..., ” — 1; adding 


these ” equations, we get 


>, y, = naa, cos Xrpf--- +455 008 Rxpt sss +b. dysin ka, + seey 


where the summation is carried from r = 0 tor = n — I. 


If we let p= & Se in the expressions for >, cos (a +78) and 


> sin(a + 78) derived in the note on p. 175, these become 


Yoos(a +kr 2") ees (« + eure) t) =O, since sin kr =0, 


~ sin (kr/n) cA. n 


dsin (attr 7 ere Ee sin (0 + aoa =O, since sin kr=o, 

n sin (kr/n) n 
except when & is a multiple of 1, for then both sin kr and sin (kr/n) are 
equal to zero and the fractional expression becomes indeterminate. But 
when & is a multiple of 1, 


> cos (« + kr 2") = > cos (a + multiple of 27) = > cos a=ncosa. 


2 27 ° : i ° me c 
sin (0 + kr =) = Ddsin (a + multiple of 27) = sin a=nsina. 


Hence we may state 


cos (0 + kr 27) 


0, except when k= 7,2n,3n,... 
= 7 COsa, When k = 7,21, 30; ... « 
: 27 
Ysin(« + kr 27) = 0, except when k= 17, 2”,37,..- 
n 
= siti a; Whe. = 1,2 8; 3 1, «sas 


, 27 
(1) If we start our intervals at x) = 0, then x, = 7 rat and 


cos Rx, = > cos (0 + kr =" 
sin Rx, = sin (0 + kr 27) 


Dr = 100 + dn + don + ss = (do + Gn + din + din +> °°) 


oO, except when k = 1, 2n,3n,... 


ncos0 =n, when k = 2,21,3%, .+.«- 


O- for all values of &. 


I 
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. Tv Tv Zin 
(2) If we start our intervals at xo = = then 4,’ = : + Fone and 


cos kx, = Yoos(e= +hr2*) = 0 except when k= 7,2n,3n,'... 
n n 


= kn {=n when k=2, 4%; 6n, J . 
Bn n|= —nwhenk=n, 3, 5”,... 


dsin kx, = Lsin(k a eioky 27) o for all values of k. 
n n 


>» = NOy — NAn + Nan — NOgn + °° * 


2 (do = On + Gan — Gan + din — > - > ); 
Subtracting the second set of ordinates, y’, from the first set, y, we 
have 


>- Dye = Dye — 9") = 0 0! + — 1’ +92 — a! + Pte + n-1— 9 n-1 
= 27 (dn + Gan + sn +--+: ); 
I Sf 
OF) Gn GsarGar <° - = Os oan u ays Be ee + yn-1— y'n-1)- 


The first set of x ordinates start at x = o and are at intervals of 2 7/n, 
and the second set of » ordinates, start at x = 7/n and are at intervals of 
27/n; thus, the period from x = 0 to x = 27 is divided into 2 ” equal 
parts each of width +/n (Fig. 91a). Hence, 

If, starting at x = 0, we measure 2 n ordinates at intervals of x/n, the 
average of these ordinates taken alternately plus and minus is equal to the 
sum of the amplitudes of the nth, 3 nth, 5 nth, . . . cosine components. 


FIG. gia. 


Thus, to determine the sum of the amplitudes of the 5th, 15th, 25th, 
. cosine components, merely average the 10 ordinates, taken alter- 
nately plus and minus, at intervals of 180° + 5 = 36°, or at 0°, 36°, 72°, 
. , 324° (Fig. 91c); therefore 
Os + dis + de + > + + = Py (Yo — Ys + Ye — Yios + Ysa — Vis0 +yer6 
— yos2 + yoss — Yso4). 
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If the 15th, 25th, . . . harmonics are not present, then 
as = Yo (Yo — Y35 + V2 — Yios + Vise — Vis0 + ors — Yos2 + Voss — sea). 


ers s C = Tv Ee Tv 27 
(3) Similarly, if we start our intervals at x = Bae then x, = ee ee: 
and 
cos kt, = dos (i au + kr 2") = 0 for all values of k, 
2n n 
f = C Tv 2s 
sin Ro >dysin ¢ ee + kr 2m) = 0 except whenk = 7,21,3n,.. 
= mwhen k= 1 57;-9%,-; 
. kv 
=nsin——= owhenk=2m4n, 6%, : 
2n 
= When Rf =3 0.71, Vl Wan & 
Ve DI = Nd + 1b, — nbs +t nbsn— ie oneal =n (do +b, —b3n+0in — bin + ts )s 
(4) Again, if we start our intervals at %’ = = + - then 
Xp" = 3u a eee and 
2n n 
008 is = Yyoos(#32 +hr27)— Oo for all values of , 
2n n 
P PA io 1 b oar 5 27 = =: 
ysin Lae = Yisin(23% +kr =) o except when R=n, 2,37, ... 
=—nwhenk= n,57, 9n,... 
.) 3hr | 
=nsin——;{= -owhenk=2n,4n, 6n,... 
2n 
= “when k= 37,71, 1E 2... « 
Te De! = 2A) — Mbp + Non — Nin + we =1n(do— bn t03n— D5n+bm — aa ). 


Subtracting the second set of ordinates, y’, from the first set, y, we 
have 


Dye — Dye = DG = Fe) = Ho = Hl HH — Hl A A — ot 
= 21 (0, — Usa + Osn = Omit © +); 
OF By— Ban tBon— bint + © == Go-To +— Wt + + +Fn1—Ho-)- 
The first set of » ordinates start at x = 7/2 and are at intervals of 
22/n, and the second set of 7 ordinates start at x = oe aE and are at 
intervals of 27/n; thus the period from x = r/2n tox =24+7/2n 
is divided into 2 ” equal parts each of width o Hence, 


If, starting at x = 1/2 n, we measure 2 n ordinates at intervals of x/n, 
the average of these ordinates taken alternately plus and minus 1s equal to the 
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sum of the amplitudes, taken alternately plus and minus, of the nth, 3 nth, 


5 nth, . . . sine components. 
Thus to determine the sum of the amplitudes, taken alternately 
plus and minus, of the 5th, 15th, 25th, . . . sine components, merely 


average the 10 ordinates taken alternately plus and minus, at intervals 
of 180° + 5 = 36°, starting at“x = 180° 4.50 =318 , 7.2, at * = 16, 


54° 00-, . 7 4 342 (Big. Otc); therefore 
bs — bis +953 -— - + + = to (y18 = Ved + Yoo — 126 =f Yiez — Yi9s + 234 
— Yo70 + 306 c— 342). 
If the 15th, 25th, . . . harmonics are not present, then 


bs = to (ys — Yea + Yoo — Yi26 +. 9162 — ios + Voss — Yoro + Y30e — sae). 
We may also note that the set of 2 ordinates measured for deter- 
mining the 0’s lie midway between the set of 2 ” ordinates measured for 
determining the a’s, so that to determine any desired harmonic we 
actually measure 4 ordinates, starting at x = o and at intervals of 
a/2n. Weuse the Ist, 3d, 5th, . . . of these ordinates for determining 
a, and the 2d, 4th, 6th, . . . of these ordinates for determining 0. 
If the higher harmonics are present, these must be evaluated first. 
The absolute term dp is obtained from the relation 
Y= tata+ta+:--:. 
We shall now illustrate the methods developed by an example. 
Example. Given the periodic wave of Fig. 89 and assuming that no 
higher harmonics than the 6th are present, we are to determine the co- 
efficients in the equation 
y= +acosx +acos2x +--+ +a,cos6x 
+hsinx+sina2x+--- + dgsin6x. 
To determine as and dg measure 12 ordinates at intervals of 30° be- 
ginning at x = 0° and x = 15° respectively (Fig. 91d); then 
407Y 


80 


255° 285° 815° 


\. 225° 


150° 180° 


Fic. 91d. 
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ds = ps (Yo — Y30 + Yeo — Yoo + - + + + Ys00 — ¥330) 
= zs (9.3 — 15.0 +17.4 — 23.0 +37.0 — 31.0 + 15.3 — 4.0 — 8.0 + 13.2 
114.2 6.0) = 0.25. 
bs = py (v5 = V5 V5 — Yio + + = EF Yore — Vos) 
= ys (13.0 — 16.0 + 19.5 — 31.0 + 35.3 — 23.8 + 10.5 + 5.7 — 10.0 
+ 14.5 — I1.0 — 0.5) = 0.52. 


234° 


270° 806° 


—10 


FIG. 9Ic. 


To determine a; and 6; measure 10 ordinates at intervals of 36°, 
beginning at x = 0° and x = 18° respectively (Fig. 91c) then 


as = q's (Yo — Yyag + Yr2 — Viog + * > - + Yess — Vs04) 
= zo (9.3—15.3 $18.8 — 32.8+33.0—15.3—1.0+9.5 —15.0+8.4) 
= —0.04. 

bs = to (vis 2a. oo — Nieg  * ei“ V 08: — ‘Ys42) 
= 7, (13.8—16.8+23.0— 36.8 +25.5—9.0—7.7 +13.4—13.2+1.5) 
= —0.63. 
407 Y 


Fic. gid. 


198 EMPIRICAL FORMULAS — PERIODIC CURVES Cuap. VII 


To determine a, and bs measure 8 ordinates at intervals of 45°, be- 
ginning at x = 0° and x = 224° respectively (Fig. 91d); then 


as = $ (yo — Yas + Yoo — Ys + °° 1+ Yer — ‘Ya15) 
= 1(9.3 — 16.0 + 23.0 — 35.3 + 15.3 + 5:7 — 13-2 + 11.0) = —0.03. 
by = EB (22.5 = Voreot Vie ae - + Yog2.5 — 337.5) 


= 4(14.5 — 18.0+ 35.0 — 27.7 + 7.7 + 8.8 — 14.7 + 3.0) = 1.08. 

To determine a; and }3; measure 6 ordinates at intervals of 60°, be- 
ginning at x = o and x = 30° respectively (Fig. 910); then 
as = § (Yo — Yoo + Vi20 — Yis0 + You — y300) 
2 (9.3 — 17.4 + 37.0 — 15.3 — 8.0 + 14.2) = 3.30. 
& (¥30 — Yoo + Yis0 — Yo10 + Ye70 — ‘Y330) 

= 1(15.0 — 23.0 + 31.0 — 4.0 — 13.2 + 6.0) = 1.97. 

To determine a, and b, measure 4 ordinates at intervals of 90°, 
beginning at x = o° and x = 45° respectively (Fig. 910); then 
dz + ag = § (Yo — Yo + M180 —I270) = F (9.3 — 23.0 + 15.3 + 13.2) = 3.70, 
sre 2 — 3rd Ss 
be — be = % (45 — Ii35 + 205 — Ya1s) = F (16.0 — 35.3 — 5.7 + 11.0) = —3.50, 

bp = — 2.08. 


bs 


To determine a; and 0; measure 2 ordinates at intervals of 180°, 
beginning at x = 0° and x = 90° respectively (Fig. 910); then 


a, + a3 + ds = 4 (M0 — Vis) = $ (9.3 — 15.3) = —3.00, .. a, = —6.26. 
b, — bs + bs = % (yoo — Yor) =% (23.0 + 13.2) = 18.10, .. 0; = 20.60. 
To determine a we have 
dy + a + Ge + G3 + Gy + Gs + Og = Wo = 9.3, - Ay = 8.63. 
Result: 


y = 8.63 —6.26 cos x + 3.45 cos 2 x + 3.30 Cos 3 x — 0.03 cos 4x 
— 0.04 cos 5 x + 0.25 cos6x + 20.60 sin x — 2.98 sin2x 
+ 1.97 sin3 x + 1.08 sin4 x — 0.63 sin5 x + 0.52sin6x. 


This result agrees quite closely with that of Art. 89, p. 184; the differ- 
ences in the values of the coefficients are due to the fact that by the 
method of Art. 89 only the ordinates at 0°, 30°, 60°, . . . , 330° are used, 
whereas by the method of this Art. a large number of intermediate ordi- 
nates are used. If the curve is drawn by some mechanical instrument, 
the present method will evidently give better approximations to the 
values of the coefficients; but the labor involved in using the computing 
form on p. 183 is much less than that used in measuring the selected 
ordinates above. 

92. Numerical evaluation of the coefficients. Averaging selected 
ordinates. Odd harmonics only. — If the axis is chosen midway between 
the highest and lowest points of the wave and the second half-period is 
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merely a repetition below the axis of the first half-period, then only the 
odd harmonics are present. If the ordinates at x = x,andx =x,+7 
are designated by y, and 4,4, respectively, then y,4, = —y,. In the 
method of averaging selected ordinates, the 2 ” ordinates are spaced at 
intervals of 7/n and are taken alternately plus and minus; then y,1, is 
at a distance r = n (r/n), or intervals, from y,, and since 1 is odd, y,+, 


will occur in the summation with sign opposite to that with which y, 
occurs, so that, e.g. 


On + Mn bs - + = 
I : 
sa (en Fo Bae ote bane! ae oe | 
eels Ae hice 
7, (0 24 e4- ) 
I 
Oar aa ick geet Yeo oa.) 


Hence we need merely divide the half-period into ” equal intervals and 
average ” ordinates. We may therefore restate our rules for determining 
the coefficients if the wave contains odd harmonics only. 

If, starting at x = 0, we measure n ordinates at intervals of r/n, the 
average of these ordinates taken alternately plus and minus 1s equal to the 
sum of the amplitudes of the nth, 3 nth, 5 nth, . . . cosine components. 

If, starting atx = 7/2 n, we measure n ordinates at intervals of r/n, the 
average of these ordinates taken alternately plus and minus is equal to the 
sum of the amplitudes, taken alternately plus and minus, of the nth, 3 nth, 
5 nth, . . . sine components. 

Furthermore, a) = 0 since the sum of the ordinates over the entire period 
1s Zero. 


80 


20 


160° 


720° 


FIG. 92. 


Example. Assuming that the symmetric wave of Fig. 92 contains no 
higher harmonics than the 5th, we are to determine the Ist, 3d, and 5th 
harmonics. Applying the above rules we have 
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ds = 4 (yo — Y35 + Yr2 — Vos + Yias) 
=}(0 — 86+ 6.3 — 27.7 + 19.0) = —2.20. 
bs = 3 (vis—Vs4 +00 — Vi2e-+ Vien) = # (11.3 —2.7 +20.5 — 25.5 + 10.7) = 2.86. 
a3 = 2 (yo — Yoo + Yi20) = = 42 Boge 5) =\/7.00; 
bs = 4 (30 — Yoo + Yis0) = 3 (10.7 Reg, = 2.27. 


rhs 0G Ost 00) a we qe 5.70: 
b, — b3 + bs = + (yo) = 20.5, -. 6; = +19.91. 
Result: 

y = —5.70cos x + 7.90 COs 3 X — 2.20 COS 5 ¥ 


+ 19.91 sinx + 2.27sin3 x + 2.86sin 5 x. 

We may compare this result with that obtained for the same curve 
by the use of the computing form on p. 187. 

If only the 1st and 3d harmonics had been present in the above wave, 
we should have 

= 3 (%0 — Yoo + 120) ; bs = 3 Ga — Yoo + 4150) 5 
a, + a3 = Yo = 0; bi — b3 = Yoo. 

If all the odd harmonics up to the ninth had been present in the above 

wave, we should have 


4 (Yo — Yoo + Y40 — Yeo + Yso — Yio + Yi200 — Yr4o + Y160) 3 
i (v10 — Y30 + Ys0 — Yr0 + Yoo — Yi10 + Yis0 — Yis0 + Yiz70) 3 
= 4 (yo — Ya5.n + ¥51.43 — Yr7.14 + Vio2.86 — Yi28.57 + Yi54.29) ; 
b, = +4 (Ni2.86 — 38.57 + 64.29 — Yoo = Yi115.71 — 141.43 + V167.14) 3 
= (yo — Y36 + V2 — Yio + Visa) 3 bs =4 (v1 — 54a + Yoo — Vie6 + Y162) 
3 + Ay = § (Yo — Yoo + Yi20)3 03 — bo = 3 (¥30 — Yoo + 150) ; 
ay + a3 +45 +47 +4 =%=0; b, — b3 +0; — by + bo = Yo0- 

Similar schedules may be formed for determining the odd harmonics 
up to any order. 

93. Graphical evaluation of the coefficients. — Various graphical 
methods have been devised for finding the values of the coefficients in 
the Fourier’s series, but these are less accurate and much more laborious 
than the arithmetic ones. The graphical methods, while interesting, are 
of little practical value in rapidly analyzing a periodic curve, so that we 
shall describe here only one of these methods — the Ashworth-Harrison 
method.* 

If, for example, we divide the complete period into 12 equal intervals 
and measure the I2 ordinates, we shall have the table 


©” | 302.60" 00 h1 20-1507 80" 210° | 240° | 270° | 300° | 330° 
Soo SIM ai AVS 4 Ys 6 yz | Ys Is 104 a 


* Electrician, Ixvii, p. 288, 1911; Engineering, Ixxxi, p. 201, 1906. Other methods 
are briefly mentioned and further references are given in Modern Instruments and 
Methods of Calculation, a handbook of the Napier Tercentenary Celebration. 


° 


ee 
ay 
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We have already shown (p. 181) that 


6a, = > COS X, = Yo COS O° + yy, COS 30° + + + + + y11 COS 330°, 


63, = Do. sin x, = yosino® + y,sin30° + ... +41 sin 330°. 


It is evident that if we consider the y’s as a set of co-planar forces 
radiating from a common center at angles 0°, 30°, 60°, . . . , the sum of 
the horizontal components is equal 
to 6a, and the sum of the vertical 
components is60;. To facilitate the 
finding of these sums we may draw 
the polygon of forces, starting ata 
point O and laying off in succession 
the ordinates, each making an angle 
of 30° with the preceding, as in Fig. 
93a (proper regard must be had for 
the signs of the ordinates). The 
polygon of forces may be constructed 
rapidly by means of a protractor 
carrying an ordinary measuring scale 
along the diameter. Then, OA, the 
projection of the resultant OP on the 
horizontal, is equal to 6 a, and OB, 
the projection of the resultant OP on 
the vertical, isequal to6;. Further- 
more, if we write a cosx + 0; sin x 
=c sin (x + ¢,), then the length 
of OP is 6c and the angle POB is ¢,. In Fig. 93a we have made the 
construction for the determination of a, 0;, c1, and ¢; for the periodic 
curve drawn in Fig. 89 using the table of ordinates on p. 184. We find 


OA =64a, = —41-4, OB =6),.= 126.0, OP.= 6¢.= 134, 
Yh POB = di = —18.1°; 


hence 
a = —6.9, by = 21.0, C= 22.3, gi: = —I8.1°. 


These results agree very closely with those obtained on p. 184. 

We may find a2 and by by laying off in succession the ordinates, each 
making an angle of 60° with the preceding; we proceed similarly in finding 
the other coefficients. A separate diagram must be drawn for each pair 
of coefficients. 

More generally, if we divide the complete period into » equal intervals 
of width 2 7/n and measure the x ordinates, then (p. 177) 
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a= >) Cos kx,=No coso+yy cos &(27) tour tea aaa wip 


2(n=1) 1 


b= >: sin kx, =o sinoO-+y; sin k (27) + + + +y,-181n k 


Hence, if we construct the.polygon of co-planar forces by starting at a 
point O and laying off in succession the ordinates, each making an angle 
2 kx/n with the preceding, then OA, 
the projection of the resultant OP 
on the horizontal, is equal to na;/2, 
and OB, the projection of the result- 
ant OP on the vertical, is equal to 
nb,/2, except when k = 0 ork = n/2, 
when we get the values nd, nbo, 
NOnj2, NOnj2, respectively. Further- 
more, the length of OP is n/2 (or n) 


Fic. 930. FIG. 93¢. 


times the amplitude c, and the angle between OP and OB gives the phase 
gx of the complete harmonic c; sin (kx + ¢,). 

Example. Analyze graphically the periodic curve in Fig. 860. 

As in the example on p. 181, we shall find the first three harmonics 
from the data 


x | 0° | 60° | 120° 180° 240° 300° 
y 0.47 1.77 2.20 ; —2.20 | —1.64 | —0.49 
Here 
6% =YMtMNt:-+> +¥ = OL; d= 0.02. 
643 =M—-—Mt--: —W= 1.95; a= 0.33. 
3a, = OA (Fig. 930) m= 3.008) “Gy = aT03: 
3b; = OB (Fig. 930) = 5.35; b= 1.78. 
3a = OP (Fig. 93d) = 6.25%. 4 = .2.08, ¢1 = 30°. 
3 a = OA (Fig. 93c) = —2.67; dz = —0.89. 
3 b, = OB (Fig. 93c) = 1.35; 0: = —0.45. 


3 @ = OP (Fig. 93¢) = 3.00; & = 1.00, ¢2 = —60° 


ART. 94 MECHANICAL EVALUATION OF THE COEFFICIENTS 203 


Result: 


y = 0.02 + 1.03 cos x — 0.89 cos 2 x + 0.33 CoS 3X 
+ 1.78 sin x — 0.45 sin 2x 
= 0.02 + 2.08 sin (« + 30°) + sin (2 x — 60°) — 0.33 sin (3x — 90°). 


Note the close agreement of this result with that obtained by the 
arithmetic method on p. 181. 

94. Mechanical evaluation of the coefficients. Harmonic analyzers. 
— A very large number of machines have been constructed for finding 
the coefficients in Fourier’s series by mechanical means. These instru- 
ments are called harmonic analyzers. The machines have done useful 
work where a large number of curves are to be analyzed. Among these 
analyzers we may mention that of Lord Kelvin,* Henrici,t Sharp,f Yule,§ 
Michelson and Stratton, |! Boucherot,§/ Mader,** and Westinghouse.tt 
We shall briefly describe the principles upon which the construction 
of two of these instruments depend.{{ 

The harmonic analyzer of Henrict. This is one of a number of ma- 
chines which use an integrating wheel like that attached to a planimeter 
or integrator §§ to evaluate the integrals occurring in the general expres- 
sions for the coefficients 


I 2 ar I 23 I 21 E 
a - + [ y dx, a <1 f y cos kx dx, =i y sin kx dx 
T Jo T Jo T J0 


given on p. 174. 
If the curve in Fig. 94a represents a complete period of the curve to 
be analyzed, then evidently 


2a 
ie ydx = area OABCDBO; 
0 


so that, if the tracing point of a planimeter is allowed to follow the curve 
OABCDBO, the integrating wheel will give the reading 2 ma, from which 
ay may be computed. 


* Proc. Roy. Soc., xxvii, 1878, p. 371; Kelvin and Tait’s Natural Philosophy. 

+ Phil. Mag., xxxviii, 1894, p. IIo. 

t Phil. Mag., xxxviii, 1894, p. 121. 

§ Phil. Mag., xxxix, 1895, p. 367; The Electrician, March 22, 1895. 

|| Phil. Mag., xlv, 1898, p. 85. 

{ Morin, Les Appareils d’Intégration, 1913, p. 179. 

** Elektrotech. Zeit., xxxvi, 1909; Phys. Zeit., xi, 1910, p. 354. 

{t The Electric Journal, xi, 1914, p. 91. 

tt Brief descriptions of all but the last of these may be found in Modern Instruments 
and Methods of Calculation, a handbook of the Napier Tercentenary Celebration, 1914. 

§§ For the principle of the planimeter and integrator, see pp. 246, 250. 
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Integrating by parts, we may write 


2 1 Qa 2a 
a=? [ yooskedx= | 7 ysin kx) ~ ied, "sin kedy=— 2 [" sinkx dy, 
TJ 0 k k 0 


oo coskxdy= ae coskxdy. 


Now if the planimeter carries two 
integrating wheels whose axes make at 
each instant angles kx and 7/2 — kx 
with the y-axis, and the point of inter- 
section of these axes is capable of 
moving parallel to the y-axis, then as 
the tracer point passes around the 
boundary OABCDBO, these wheels give readings proportional to 


J sinkx dy and fos (2 — ke) dy = [ cos kx dy, 


from which the values of a, and 6; can be found. 

In one form of the instrument the curve is drawn on a horizontal 
cylinder with the y-axis as one of the elements. A mechanism is attached 
to a carriage which moves along a rail parallel to the axis, by means of 
which a tracer point follows the curve while the cylinder rotates; the 
mechanism allows the axes of the integrating wheels to be turned through 
an angle kx while the cylinder ro- 
tates through an angle x. Coradi, 
the Swiss manufacturer, has per- 
fected the instrument so that sever- 
al pairs of coefficients may be read 
with a single tracing of the curve. 


b= = Ly ysin kxdx = | —Lycoske| 
J 0 


FIG. 94a. 


ire 
referencé o 


Fic. 94). FIG. 94¢. 


The Westinghouse harmonic analyzer.—This machine, constructed 
by the Westinghouse Electric and Mfg. Co., is particularly useful in 
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analyzing the alternating voltage and current curves represented by a 
polar or circular oscillogram. 

Fig. 945 gives one period of a periodic curve drawn on rectangular 
coérdinate paper. In Fig. 94c, the same curve is represented on polar 
coérdinate paper. This is done by constructing a circle of any convenient 
radius, called the zero line or reference circle and locating any point P 
by the angle @ = x and the radial distance r = yfrom the zero line. Thus 
the points marked P, A, and B in Figs. 94b and 94c are corresponding 
points. If only the odd harmonics are present, the second half-period 
of the curve in Fig. 946 will be a repetition below the x-axis of the first 
half-period; in this case, the diameters at all angles of the curve in Fig. 94¢ 
will be equal, and equal to the diameter of the reference circle. The re- 
lation between 7 and @, 


r = f(@) =a,cos@é +a,cos26+--+- +a,cosk0+:--- 
+ bsiné+d:sin20+ +--+ +b:sinkO+---, 
is the function to be analyzed. This is done as follows. 


The circular record of the periodic curve, drawn by hand from the 
rectangular record or directly by the circular oscillograph,* is transferred 


Fic. 94d. 


to a card of bristol board and a template is prepared by cutting around 
the curve. In the initial position the template M (Fig. 94d) is secured 
on a turntable T so that the axis 9 = o lies under the transverse cross-bar 
B. The turntable is set on a carriage D which slides on the rails ZL. The 


* The Electric Journal, xi, 1914, p. 262. 
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carriage is given an oscillatory motion by the motion of a crank-pin P 
(Figs. 94e, 94f) attached to a rotating gear G and sliding in a transverse 
slot S on the bottom of the carriage. The carriage thus has a simple 
harmonic motion whose amplitude is the crank-pin radius R. By means 
of a crank and a simple arrangement of gears, the carriage makes k com- 
plete oscillations while the template makes one revolution, when deter- 
mining the kth harmonic. 


Fic. 94e. Fic. 94f. 
' The cross-bar B is attached to the oscillating carriage; this bar carries 
a pin C held in contact with the edge of the template by means of springs, 
so that the bar has a transverse motion as the template revolves. Re- 
ferred to a pair of axes xx and yy, the motion of the end of the bar, 
Q(x,y), may be said to consist of two components, viz., the transverse 
motion of the bar, x = r = f(6), the function to be analyzed, and the 
simple harmonic motion of the carriage, 


(I)hiy= Ksinke or  e)ey = Koin (x0 _ 4 = —Rcos kd, 


according as the carriage is started with the slot S in the dotted position 
of Fig. 94e or of Fig. 94f. A planimeter is attached with its tracing point 
at Q. This point then describes compound Lissajous figures whose areas 
A, and Az may be read from the integrating wheel of the planimeter. 

dy 


Mae RR cos k@ and from (2) 2 = Rk sin k@, hence 


Now from (1) 


2a 2a 24 
Ages f “ay = f 7+Rk cos k0 dé = Re [ rcos k0 dé = Rk ax, 
0 0 


| 20 Qa 
A, = f iay = i rRRk sin k8 d@ = Rk di rsin k0d0 = Rkrb,, 
\ 0 0 


using the formulas for a; and b; on p. 174. 
_ Ay 2 As, 
Rkr’ * Rk 
Gears are provided to analyze for all even and odd harmonics from 1 
to 50, and the shifting of the gears is a very simple matter. 


Therefore ax 
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EXERCISES. 


1. Sketch the periodic curves 
y=2cosx; y=cos2x%; y =2cosx+cus2%.' 
2. Sketch the periodic curves 
y=1+sin4x; = —jsin2x; y = 4sin3x; 
=1+sinx —}sin2x+ 4sin3x. 
3. Sketch the periodic curves 
y = 2sin (« —40.5°); y = sin (2x + 72.3°); y = 4sin (3 x — 90°); 
y = 2sin (x — 40.5°) + sin (2x + 72.3°) + 3 sin (3 x — 90°). 
4. Sketch the periodic curve 
y =cosx + 0.4cos3% + 0.5 sinx — 0.5 sin3 x. 
5. Sketch the periodic curve 
y = cosx + 0.4cos3 x — 0.2cos5x + 0.5sinx — 0.5 sin3 x — 0.3 sin5 x. 


6. By use of the formulas on p. 177 and the direct method illustrated on p. 179, 
determine the coefficients of the third and fourth harmonics of the periodic curve in 
Fig. 89; use the table of ordinates on p. 179. 

7. Determine the first three harmonics of the periodic curve given by the following 
data; use the computing form on p. 180. 


C 0° |-+.60° 120° | 180° 2% 
y —0.85 | 0.95 0.72 2.75 Sy 


8. Determine the first six harmonics of the periodic curve given by the following 
data; use the computing form on p. 183. 


fe 30° | 60° 240° 300° 

8) BS) BS, 14 15 
9. Determine the first twelve harmonics of the periodic curve given by the following 
data; use the computing form on p. 185. (Thecurve is a graphical representation of the 


diurnal variation of the atmospheric electric potential gradient at Edinburgh during the 
year 1912.) ; 


300° 


—2.20 


0° 


—18 


210° 
10 


180° 
II 


120° 
22 


go° 
= 


150° 
22 


270° 
12 


330° 
—TI 


Pee O15, 1 330° [h452. |. 60" (750) | 90" | Tos” [120° |.135° [50° | 165" 
By || 8) XO aS) |) || Soy el Sh fl 22 24 2220 

ealt8o02 195, | 210° | 225- | 2A07 | 255° |'270- | 2852 (200° (315 .| 330. | 345- 
y II 3 10 16 14 12 12 18 15 9 —1I | -—7 


10. Devise computing forms for determining the even and odd harmonic coefficients 
using 8 and 16 ordinates respectively. 

11. Determine the odd harmonics up to the fifth for the symmetric periodic curve 
given by the following data; use the computing form on p. 187. 


fe} 


120 
1004 


x 60° | go° 150° 
y 660 940 554 


12. Determine the odd harmonics up to the fifth for the symmetric periodic curve 
from which the following measurements were taken; use the computing form on p. 187. 


0° 3 0° 
(0) 676 


° ° 


x 


y 


oO 
oO 


150° 
2 


30m | 60° 


go° | 120 
4 9:5 


9 3 
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13. Determine the odd harmonics up to the eleventh for the symmetric periodic 
curve from which the following measurements were taken; use the computing form on 


p. 190. 
x 0° | 1s? 30°. |.45°-|[800° | 75" | Gol eO5 e203 5 3 TSO" LOS. 
y | (@) TA meas 52 60 40 27 30 15 18 6 14 


14. Determine the odd harmonics up to the seventeenth for the symmetric periodic 
curve from which the following measurements were taken; use the computing form on 


p. 192. 
% 0° 10° |* 20°) 30° {40° | 50% (s G0"|, 70° 180; )90 | Boos IO: 
y | 0) 5 9 21 230) || 2 oy asie | M2e i 2%) || Ze || ah 
120° | 130° | 140° | 150° | 160° | 170° 
38 ||.-46°-[-45 | 30 | BP] 29 


15. Determine the first three harmonics for the periodic curve from which the fol- 
lowing measurements were taken; use the method of selected ordinates in Art. 91; 
assume that all higher harmonics are absent. 


0% NGOS We45 a 60cel 90n | T20rg N35. 1500 | 180 4 2TO 22542400) eo 708 
y | 10.0 | 5.0 | 5.3 | 7.2 | 6.0 |—6.8 |—10.9 |—8.9 | 10.0 | 18.5 | 10.7 |—3.4 |—25.9 
oy 300 315 330 
OL MN WES esa aes 5-1 


16. Determine the first three harmonics for the periodic curve drawn in Fig. 86); 
use the method of selected ordinates in Art. gr. 

17. Determine the first six harmonics for the periodic curve drawn in Fig. 89; use 
the method of selected ordinates in Art. 91; assume that all higher harmonics are 
absent. 

18. Determine the first and third harmonics for the symmetric periodic curve given 
by the following data; use the method of selected ordinates in Art. 92; assume that all 
higher harmonics are absent 

x On | 30° | 60° | 90° | 120° | 150° 
y oO 62.9 66.5 22.4 14.9 33:3 

19. Assuming that the harmonics higher than the fifth are negligible, determine the 
odd harmonics of the symmetric periodic curve from which the following measurements 
were taken; use the method of selected ordinates in Art. 92. 


ae 0° 30° 60° go° 120° 150° 
“y Oo 676 | 660 940 1004 | 554 
x 0° 18° 36° Bynes |) Gee 90° |, 50820) 126° | brag 1622 
y o | 470 | 719 | 678 | 702 | 940 | 1086 | 920 | 639 | 375 


20. Use the method of selected ordinates in Art. 92 to determine the ninth harmonic 
of the curve given by the table in Ex. 14. 
21. Analyze graphically the curve in Ex, 7, 


CHAPTER VIII. 
INTERPOLATION. 


95. Graphical Interpolation. — Having found the empirical formula 
connecting two measured quantities we may use this in the process of 
interpolation, 1.e., in computing the value of one of the quantities when 
the other is given within the range of values used in the determination of 
the formula. It is the purpose of this chapter to give some methods 
whereby interpolation may be performed when the empirical formula is 
inconvenient for computation or when such a formula cannot be found. 

Let the following table represent a set of corresponding values of two 
quantities 


where y is a known or an unknown function of x. Our problem is to find 
the value of y = y, for a value of x = x; between x and xp. 

A simple graphical method consists in plotting the values of x and y 
as codrdinates, drawing a smooth curve through or very near the plotted 
points, and measuring the ordinate y, of the curve for the abscissa x,. 
The value of y, thus obtained may be sufficiently accurate for the purpose 
inhand. Thus from the curve in Fig. 72), weread t = 10, A = 77.0, and 
t = 30, A = 45.0. If we use the empirical formula derived on p. 133, 


A = 100.1 ¢°-%5*, or log A = 2.0005 = 0.0115 f, 


we compute ¢ = 10, A = 76.8 and ¢ = 30, A = 45.2. By comparison 
with the table on p. 132 we note that the measured values of A fort = 10 
and ¢ = 30 agree about as closely with the computed values as the neigh- 
boring observed values agree with their corresponding computed values. 
Here, the last significant figures in the values of A were used in construct- 
ing the plot. 

On the other hand, in Fig. 71c, we read v = 40, p = 10.00, whereas 
the empirical formula on p. 131 gives v = 40, p = 9.42. The residual 
is 0.58, much larger than the residuals in the table on p. 130 for neighbor- 
ing values of v. Here, the plot was constructed without using the last 
significant figures in the values of the quantities. It is of no advantage 
to construct a larger plot since the curve between plotted points is all 
the more indefinite. 

209 
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For most problems the arithmetic or algebraic methods to be explained 
in the following sections give much better results. 

96. Successive differences and the construction of tables. — Given a 
series of equidistant values of x and their corresponding values of 4, 


x Xo x1 Xe x3 Xn 
Xo xo th xX + 2h xo + 3h Xo + nh 
ay Yo y1 y2 ¥3 ns 
we define the various orders of differences of y as follows: 
Ist difference = A‘: ad) = Yi — Yo, @1 = Yo — Vy - + y On1= Yn — Vn—-13 
2d difference = A?: bo = a; — A, by = 09 = GQ, 3 1 5 ye y= AQn—2 — An-1; 
3d difference = A?: Cc = b; at bo, qQ = by = bi, Soo eee Sa bn—3— bn-2} 
kth difference = A*: Ry = jo — fu, Ri = je — jy - 
These may be tabulated as follows: 
x y Al A? A3 A‘ Ak 
Xo = Xo Yo 
a 
Mm =x +h VM bo 
aq co 
X. =X + 2h Y2 bi do 
a2 C1 
x3 = Xo + 3 h 3 be . ‘ 
, a3 5 5 . ko 
X= % +4 h M4 ° 
: ky 
Xn = ao ia — i)h Ves 
Qn—-1 


Xn = X0 + nh Vn 


where a quantity in any column of differences is written between two 
quantities in the preceding column and is equal to the lower one of these 
minus the upper one. 

We may apply the above definitions in the formation of the differences 
of y when y = f(x); thus, 

Ay = f(x + h) — f(x) = Af(x); Ay = Af(x + h) — Af(x) = A*f(x); 
etc. £.g., if 


y= —2x+2, Ay = [e+ hk)? — 2 (e+ h) + 2] — [x — 2474-2] 
=2hx + (hk? — 2h); 

[2h (x +h) + ( — 2h)| — [2hx + (2 —2h)] 

= a) [(P. 


We note that A?y = 2 h?, so that the second differences are constant for 
all values of x. 


A*y 
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Similarly, if y = x" where 7 is a positive integer, 


ee Bee ah 


Aye h)* = x= x nx h + 


= nxf + eo! deo iy 
aty=[n (e+ hth + ne) (Mm + | 
= | mh a oe an-2fe te | 
2 


n(n—1)x™ 72+ --.-, 
A’y= n(n —1)(n—2) x74 ---, 


Any= n(n —1)(n—2)...3+2+1k* =(nh"; 


hence the nth differences of x”, where 1 is a positive integer, are constant, 
and hence the mth differences of any polynomial of the mth degree 


Axi Bee Cx ea he Le 


where » is a positive integer, are constant. If in forming the differences 
of a function some order of differences, say the mth, becomes approxi- 
mately constant, then we may say that the function can be represented 
approximately by a polynomial of the mth degree, where n is a positive 
integer. , 

The formation of the differences for various functions is illustrated 
in the following tables: 


GQjy=x = (2) y= % 
x y Al A? A3 | x y Al A? 
I I 5-16 137-39 
7 4.03 
2 8 12 5.21 141.42 0.08 
19 6 4.11 
a 27 18 5.26 145.53 0.08 
37 6 4.19 
4 64 24 5.31 149.72 0.08 
61 6 4.27 
5 125 30 5.36 153-99 0.08 
91 6 4.35 
6 216 36 5.41 158.34 0.08 
127 6 4-43 
7 343 42 ; 5.46 162.77 0.08 
169 6 4.51 
8 512 48 5.51 167.28 0.09 
217 4.60 
9 729 5.56 171.88 
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(3) y= Vx (4) y = Vx 
x y Al a2 x y Al 
20 2.7144 611 8.4856 
445 46 
21 2.7589 —14 612 8.4902 
431 46 
22 2.8020 —12 613 8.4948 
419 46 
23 2.8439 =13 614 8.4994 
406 46 
24 2.8845 ily 615 8.5040 
395 46 
25 2.9240 616 8.5086 
(5) Train-resistance (6) Speed of a vessel 
j j AS 2 Lip | t A2 3 
i 1 in 
ors ne a ; oe er hr. | 208se-power i 3 
20 5.5 8 1,000 
3.6 400 
40 9.1 2.2 9 1,400 100 
5.8 500 fe) 
60 14.9 Deli 10 1,900 100 
Fc) 600 50 
80 22.8 2.4 II 2,500 150 
10.5 750 50 
100 238) DP 12 3,250 200 
127 950 50 
120 46.0 13 4,200 250 
1200 100 
14 5,400 350 
1550 100 
15 6,950 450 
2000 50 
16 8,950 500 
2500 
17 11,450 
Cay —slogs (8) y = log sin x 
x y Al x | y Al a2 A8 
500 2.6990 Lear Os 8.2419-10 
8 669 
501 2.6998 i ikoy. 8.3088-10 —89 
9 580 20 
502 2.7007 T° 20' 8.3668-10 —69 
9 511 ‘16 
503 2.7016 30! 8.4179-10 —53 
8 : 458 8 
504 2.7024 1° 40 8.4637-10 —45 
; 9 413 10 
505 2.7033 1° 50’ 8.5050-10 —35 
9 378 
506 2.7042 Boao! 8.5428-10 


In the above tables we note the following: 

In (1); -ye= + and: A® is constant: 

In (2), y = x* and A? is constant since we have carried the work to two 
decimal places and A* does not sensibly affect the second decimal place. 
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If the computation had been carried to six decimal places, A? would not 
be constant but A’ would be. 

In (3), A? is approximately constant, so that if we desire to work to 
four decimal places, Wx could be represented by a polynomial of the 
second degree within the given range of values of x. 


In (4), A‘ is approximately constant so that Wx could be represented 
by an equivalent polynomial of the first degree. 

In (5) and (6), A? and A® are approximately constant, so that R may 
be approximately represented by a polynomial of the second degree in V, 
and I by a polynomial of the third degree in V. 

In (7), log « may be approximately represented by a polynomial of the 
first degree, and in (8), logsinx by a polynomial of the third degree 
within the given range of values of x. 

In general, it is evident that we may stop the process of finding suc- 
cessive differences much sooner the smaller the number of digits required 
and the smaller the constant interval h. We should stop immediately if 
the differences become irregular. 

The formation of differences is often valuable where a function is to 
be tabulated for a set of values of the variable. Thus, suppose we wish 
to form a table for y = rx?/4, expressing the area of a circle in terms of 
the diameter, for equidistant values of x. Since we have a polynomial 
of the second degree, A’y is constant, and if h = 1 and the work is to be 
carried to 4 decimal places, we need merely compute the values of y for 
x = 1, 2, 3 and form the corresponding differences; proceeding back- 
wards, we repeat the value of A?’y = 1.5708, add this to Ay = 3.9270 and 
get 5.4978, add this to 7.0686 and get 12.5664, which is the value of y for 
x = 4. We proceed in the same manner to get the values of y for suc- 
cessive values of x. 


x y =1x?/4 Al a2 x y =1x?/4 Al A? 

I 0.7854 69 3739.28 
2.3562 109.17 

2 3-1416 1.5708 70 3848.45 1.57 
3.9270 110.74 

a 7.0086 1.5708 71 3959.19 1.57 
5.4978 112.31 

4 12.5664 1.5708 72 4071.50 1.57 
7.0086 113.88 

2) 19.6350 73 4185.38 


For larger values of « where we wish to work to two decimal places 
only, we take A’y = 1.57 and proceed as above. 

Suppose we wish to tabulate the function y = x3. Here A® is con- 
stant so that we merely compute the part of the accompanying table in 
heavy type. Then we extend the column for A’ by inserting 6’s, extend 
the columns for A? and A! by simple additions and subtractions, and thus 
determine the values of «x? for all integral values of x. 
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x y= x8 Al A2 A3 
—I -I 
I 
O O (0) 
I 6 
I I 6 
7 6 
2 8 12 
19 6 
3 27 18 
37 6 
4 64 24 
61 6 
5 125 30 
gI 
6 216 


The same procedure may be followed in the construction of a table for 
a function where a certain order of differences is only approximately con- 
stant. Thus, in forming table (4) of cube roots, we note that for that 
portion of the table Ay is approximately 0.0046 so that we can find the 


values of “x by simple additions; we must check the work by direct 
computation every few values in order to find when A’y changes its value. 

97. Newton’s interpolation formula. — We shall now express the 
value of y for any value of x. From the definitions of successive differ- 
ences we have 


V1 = Yo + a; Yo = V1 + G1 = (Yo + Go) + (Go + bo) = Yo + 2 yo + Dy; 
Ys = Yo + de =(yo + 2 a0 + bo) + (do + 2 bo + Go) = Yo + 3G + 300 + &; 
Ys = ¥3 + Gs = (Yo + 3A + 300 + co) + (Go + 350 + 3 Co + do) 

= yo + 44) + 6b) + 469 + dy; 


We note that the coefficients are those of the binomial expansion, and 
this suggests that 


n(n — 1) 


Ya = Wot ray + "= Dp, 4 WO = A) 2) eee) 


Cot:::, (I) 


where 1 is a positive integer. If this equation is true, then, replacing y 
by a, the first difference, we may also write 


Hea yeas! Dog p BO 2 ay eo. 
Ynti = Yn + An = Yo + (um + 1) a + [P= + nlite 


n(n —1)(m—2) , n(n—1) 
(| (See eae 


= 0+ (nt 1) a +E", ee ee 
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where the coefficients are again those of the binomial expansion with n 
replaced by » + 1. Thus we have shown that if equation (I) is true for 
any positive integral value of 1, it is true for the next larger integral value. 
But we have shown (I) to be true when = 4, therefore it is true when 
nm = 5;since itis true for = 5, therefore it is true form = 6;etc. Hence 
(I) is true for all positive integral values of n. 

Now if some order of differences, say the kth order, is constant, ‘.e., 
A*'y = ko, then y is a polynomial of the kth degree in m, and equation (I) 
may be written 


DeoBy Cr» => Kut = 30+ mag + =) 4 Yee 
n(n—1I)...(nm—k+1) 
[k 
The right member of this equation is also a polynomial of the kth degree 
in m, and since these polynomials are equal for all positive integral values 
of n (t.e., for more than k values of 1), they must be equal for all values of 

n, integral, fractional, positive, and negative. 

Hence if the &th order of differences is constant, we have 
Yn = 90+ 1p +25, 4 2 oil olen 
n(n—1)...(n—k+1) 

[Rk 
for all values of nm. This fundamental formula of interpolation is known 
as Newton's interpolation formula. In this formula, yo is any one of the 
tabulated values of y and the differences are those which occur in a line 
through y) and parallel to the upper side of the triangle in the tabular 
scheme on p. 210. 

Newton’s formula is approximately true for the more frequent case 
where the differences of some order are approximately constant; all the 
more so if ~<_1.. We can always arrange to have  < 1; for if we wish 
to find the value of y = Y for x = X, where X lies between the tabular 
values x; and x;41, we use Newton’s formula with y; and the correspond- 
ing differences a;, bi,¢c;, . . . ,sothat X = x; + nhandn = aos <a ne 
The values of the binomial coefficients occurring in the formula have 
been tabulated for values of 2 between 0 and 1 at intervals of 0.01.7 

Let us now apply Newton’s formula to the illustrative difference- 
tables (1) to (8). é 


* The ordinary interpolation formula of proportional parts disregards all differences 
higher than the first, so that y = yo + do, where n = (X — xo)/h. This simple formula 
will often give the desired degree of accuracy if the interval h can be made small enough. 

{+ See H. L. Rice, Theory and Method of Interpolation. 


+ 


Ro: 


a Ro (N) 
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(1) To compute (3.4)*; yo=27, h=1, n=(3.4—3)/1=0.4; 


.4) (—0.6 0.4) (—0.6) (—1.6) 
(3.4)8=27+ (0.4) (37) + oe (24) + (0-4) ( 6 (6) 
= 39.304. 
(3) To compute V23.5; yo = 2.8439, h = 1, n =(23.5 — 23)/1 = 0.5; 
V'23.5 = 2.8439 + 4 (0.0406) + } (0.0011) = 2.8643. 
If we use the ordinary interpolation formula of proportional parts, 


VW23.5 = 2.8439 + 3 (0.0406) = 2.8642, which would be correct to 
three decimals only. 


(4) To compute V612.25; yo=8.4902, h=1, n=(612.25—612)/1=}; 
V612.25 = 8.4902 + } (0.0046) = 8.4914. 


(5) To compute R when V=65; Ro=14.9, h=20, n= (65 —60)/20=}; 
R= 14.9 + 3 (7.9) — #5 (2.4) = 16.7. 

(7) To compute log 501.3; yo= 2.6998, h=1, n=(501.3—501)/I =0.3; 
log 501.3 = 2.6998 + 0.3 (0.0009) = 2.7001. 

(8) To compute log sin 1° 16’; y = 8.3088 — 10, h = 10’, n= 
(1°16) = 17.10) /105 = 0.6; 

log sin 1° 16’ = (8.3088 — 10) + 0.6 (0.0580) — 0.12 (— 0.0069) 
+ 0.056 (0.0016) = 8.3445 — I0, correct to 4 decimals. 

If we use the ordinary formula of proportional parts, we have 
log sin 1° 16’ = 8.3088 — 10 + 0.6 (0.0580) = 8.3436 — 10, correct to 2 
decimals only. 

If the value of x for which we wish to determine the value of y is near 
the end of the table we may not have all the required differences. To 
take care of this case Newton’s formula is slightly modified. If we 


invert the series of values of x in the tabular scheme on p. 210, and form 
the differences, we have 


x y Al a2 A3 As 
Xn Yn 
An 
Xn—-1 Yn—-1 ; 
x4 M4 
—d3 
X3 Ms be 
—a2 Cr 
Xo Ye by dy 
7a —60 
x1 al bo 
—d 
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Starting at y, and applying Newton’s formula, we get 


a = een (ay) + Os 4 BE) (+ rar 
Se et Pes 


Comparing the result with the scheme on p. 210, we note that the differ- 
ences are those which occur along a line parallel to the lower side of the 
triangle in that scheme. Here ys is any value of y, and if X lies between 
x4 and x3, then X = x4 — nh, and n = (x, — X)/h. 

Example. To compute V24.8. In table (3), ys = 2.9240, h = 1, 
j= (25 24.8)/ 1 = 0.2; 

“. W24.8 = 2.9240 — 0.2 (0.0395) sp 22098) (0.0011) = 2.0162. 

If a series of corresponding numerical values of two quantities are 
given, we may use Newton’s formula for finding the polynomial which 
will represent this series of values exactly or approximately. For this 
purpose we replace m by (x — x) /h. 

Mus sitietapie (1), b=" 9, X= 1, 2 = 4 — 1: 


= = Cai?) (et) (ven?) 3) 
y= 1 + (x lana Hes Mi aa I2+ B O= 


Intable (5), h=20, Vo=20, n = ———=—-I; 


V & a :) & i 2) 
20 20 
R=55+(5—1)3.6+ E 2 
= 4.1 + 0.015 V + 0.00275 V?. 
The values of R computed by this formula agree quite closely with those 
in the table. 
In table (6),4 = 1, Vo = 10, » = V — to; 


I = 1900 + (V — 10) 600 + Wess ee a) 150 


fe (V — 10) —- 11) (V — 12) =o 
= —6850 + 2042 V — 200 V? + 83 V?. 
The values of J computed by this formula agree quite closely with those 
in the table; thus, V = 12 gives I = 3254. 
Various formulas of interpolation similar to Newton’s have been de- 
rived which are very convenient in certain problems. Among these 
may be mentioned the formulas of Stirling, Gauss, and Bessel.* 


* For an account of these formulas, see H. L. Rice, Theory and Practice of Interpola- 
tion, and D. Gibb, Interpolation and Numerical Integration. 
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98. Lagrange’s formula of interpolation. — Newton’s formula is 
applicable only when the values of x are equidistant. When this is not 
the case, we may use a formula known as Lagrange’s formula. Given 
the following table of values of x and y, 

x | _o _|_t a3 yt | |e 
y | on Yo Ys 4 Ys Ms 
we are to find an expression for y corresponding to a value of x lying be- 
tween a; anda,. We take for y an expression of the (n — 1)st degree in 
x containing m constants, and determine these m constants by requiring 
the 2 sets of values of x and y to satisfy the equation. But instead of 


assuming the form y = A + Bx + Cx? +--+ + Nx", we may assume 
the equivalent form 
y = A (x — d) (x — a3) («© — as) . . . (% — Gn) 
+ B(x — a) (« — a3) (wv — ay)... (4% — An) 
+ C(x — a) (x — ae) (x — a4)... (% — Gn) 
+ N (x — a) (x — de) (% — a3)... (& — Gn), 
where the 7 terms in the right member of the equation lack the factors 
(% — a), (x — de), . . . (« — a,) respectively. 
Since (a1, yi) is to satisfy this equation, 
yi = A (a, — de) (Qi — G3) (@; — ay) . . . (4) — Gn), 


since all the other terms contain the factor (a; — a;) and therefore vanish. 
Similarly, 


yo = B (dy — ay) (dz — a3) (€2 — a4) . . . (G2 — Gn), 
ys = C (a3 — a1) (@3 — ae):(@3 — a4) . . . (€3 — Gy), 
Na = N Cae Be (dn i a2) fon a a3) . Shas ie es Ge). 
Hence, , 
ee 
(a1 — dy) (a1 — a3) (1 — ay) . . . (Q) — ay)’ 
Te er elc., 
(dz a) (de a3) (de = ds) ose (a2 = Qn) 
and, finally, 
fete (x — a2) (x—4s) .. « (xan) (x—a1) (x—a3) . . . (x—an) 
(@1—*t2) (@1— as) . . . (@1—Gn) |?” (@2—a) (G2—as) . . . (Ga—Ga) 
ee (=p) Gere Gy) 2, (Stee 
7 Ga) Gea ane 


We note that in the term containing 4,, the numerator of the fraction 
lacks the factor (x — ax) and the denominator lacks the corresponding 
factor (a, — ax). Lagrange’s formula is in convenient form for logarith- 
mic computation. 
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Example. In the table on p. 132 we have 


A 68.7 64.0 44.0 39.1 
and we are to find the value of A when t=27. Using Lagrange’s formula, 
= (27-17) (27-31) (27—35) (27—14) (27-31) (27-35) 
A= 87 477) (gst) 435) 1 7 —14) (731) 735) 
(27 —14) (27-17) (27—35) , ,, , (@7—14) (27-17) 7-31) 
(3114) (31-17) (1-35) | >?” (35-14) (35-17) (35-31) 
= —20.5 + 35.2 + 48.0 — 13.4 = 40.3, 
which agrees exactly with the observed value. 
Example. In the table on p. 157 we have 
t 0.1 0.2 0.4 0.8 
1 2.48 | 2.66 ~ 2.58 | 2.00 
and we are to find the value of 7 when t = 0.3. Using only the values 
= 0.2 and t = 0.4, 
: (0.3 — 0.4) (0.3 — 0.2) 
ae (0.2 — 0.4) tee (0.4 — 0.2) 
Using all four values of t, 7 = 2.68. Using the empirical equation 
t = 4.94 e71-7t — 2.85 e-3-76¢ (on p. 159), we get 7 = 2.66. 

Gauss’s interpolation formula for periodic functions. — When the data 
are periodic we may find the empirical equation as a trigonometric series 
by the method of Chapter VII and use this equation for purposes of in- 
terpolation, or we may use an equivalent equation given by Gauss: 


+ 44.0 


= 1.33 + 1.29 = 2.62. 


‘et sin 4 (x — a2) sin} (x — a3) . . . sin } (x — ay) 
ame? sin 4 (a, — dz) sin} (a; — a3) . . . sin} (1 — Gn) 
sin 4 (x — a) sin} (x — a3) . . . sin} (x — an) 
se ain as Sues he 
ae 


It is evident that y = y,; when x = a, y = y2 when x = a, etc., so that 
the equation is satisfied by the corresponding values of x and y. 
99. Inverse interpolation. — Given the table 
a Xo | x1 Xe | x3 | Gn 
ty Yo By y2 ¥3 Yn 
we may wish to find the value of x corresponding to a given value of y. 
If the values of x are equidistant we may use Newton’s interpolation 
formula. Here we know yn, Yo, Go, 0, Co, . - - ,» and substituting these 
values in the formula we have an equation which is to be solved for n. 
If only the first order of differences are taken into account, then 
| the ordinary formula for inverse inter- 
0 
polation by proportional parts. 


Yn = Yo + Nd, and n = 
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Example. In table (7), given log « = 2.7003, to find x. 
n= 270083 oe = 0.56,andx = % + nh = 501 + 0.56 (1) = 501.56. 
If only the first and second differences are taken into account, then 
(He 
2 


y bo, a quadratic equation which can easily be 
/ 


Yn = Yo + Nd + = 

solved for n. 
Example. In table (5), given R = 27.3, to find V. 

Here 27.3 = 22.8 +n (10.5) + ese (2.2), 

or Linv+o94n — 4.5 = 0; 

hence » = 7 = 0.455 and x = Vo + nh = 80 + (0.455) 20 = 89.1. 

The empirical formula R = 4.62 — 0.004 V + 0.0029 V? on p. 149 
gives Vo ="860.1, ck =" 27-3: 

But if the third and higher orders of differences have to be taken into 
account, the method would require the solution of equations of the third 
and higher degrees. In such cases as well as in the case where the values 
of x are not equidistant, we may use Lagrange’s formula and merely in- 
terchange x and y; 1.¢., 


(Vie=aGs) (Gea Cs) eat we (y—ai)Ww-—as)... 
(a) — dz) (a, — ds) . . Prete ne eae nae 7 ek ae 


x= X11 


Example. In table (8), given log sin x = 8.3850 — Io, to find oe 
Using only the following values, 


log sin x 8.3088 — Io | 8.3668 — 10 | 8.4179 — 10 


x 70’ 80’ go’ 
we have 
RS (0.0182) (— 0.0329) , (0.0762) (—0.0329) 
(—0.0580) (—0.1091) (0.0580) (—0.0511) 
eae! (0.0762) (0.0182) 


(0.1091) (0.0511) 
= 70’ (—0.0946) + 80’ (0.846) + 90’ (0.249) 
=) 83:47" == 1° 23.47", 
We may also use a method of successive approximations as follows: 
From Newton’s formula we write 
y — Yo ' 
a +4(n— Tho titn— 1) (1 — 2)e+ Sees 
Applying this to the above example, and taking only the first differences 
into account, we get asa first approximation, 
onal A eal (8.3850 — 10) — (8.3668 — 10) _ 182 
do 0.0511 ast E 


n= 


= 0.356. 
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Taking also second differences into account and introducing the value of 
m, for m in the denominator, we get as a second approximation, 

~ ye sf 0.0182 18 
do +3(m—1)b) 0.0511 + 0.0017. 52 


n Bee 
2 3 = 0-345: 


We may continue in this way approximating more and more closely to the 
value of x. In this example it will be unnecessary to carry the work to 
third differences since A® is negligible. Hence 


nm = 0.345, and x = % + nh = 1° 20’ + (0.345) (10’) = 1° 23.45’. 
We may check this by direct interpolation. Here 


yo = 8.3668 — 10, h = 10’, and n = 0.345; 
hence, 


y = 8.3668 — 10 + 0.345 (0.0511) — 0.113 (—0.0053) = 8.3850 —10. 


Example. Find the real root of the equation x3+5x—1=0. 
We form a table of differences of the function y = #8 + 5x -— 1, 


c y Al A? As 
—2 —I19 
I2 
—I 7, —6 
6 6 
0) —I (0) 
6 6 
I 5 6 
12 6 
2 17 12 
24 
3 4! 


The root lies between x = 0 and x = I, and we are to find the value of x 
when y = 0. Using the method of successive approximations we have 


evi Orit iE 


n= a @ = B = 0.1667, 

= eet = : a2 0.1428 
tim —tb 6F4G-6 7 °°” 

= a = : = 49 ~ 0.1968 
ue ~ do+ 4 (ne — 1) bo +3 (m2 — 1) (m — 2) 6—45+48 249 Sesh 

= BS ea it 
4 ay + 3 (ms — 1) bo + 3 (ms — 1) (Ms — 2) 6 — 2.4096 + 1.4483 

I 
= —— = 0.1985. 
5.0387 is 

Hence, x = xX) + nh = 0.1985. 


From the table 
0.1984 
— 0.00019 


a9 0.1985 | 0.19845 
y | 0.00032 | 0.00006 


we note that x = 0.1984 is the root correct to 4 decimals. 
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EXERCISES 


1. Tabulate the values and differences of the following functions; / is the common 
interval. 
(a) x’, from x = 5tox = 12 whenk = 1;and fromx = 3tox = 3.1 when / = 0.01. 
(b) Vx, from x + 1 to x = 10, when h = 1, and from x = 563 to x = 570 when 
h=1. 


(c) = from x = 60 to x = 70 when # = 1, and from x = 260 to x = 262 when 


POD, 


(d) = (volume of a sphere), from D = 1 to D = 1.8 when h = 0.1. 


(e) log x, to 4 decimals, from x = 356 to x = 362 when h = 1. 

(f) tan x, to 4 decimals, from x = 32° to x = 33° when h = I0’. 

(g) log cos x, to 4 decimals, from x. = 88° 10’ to x = 89° 20’ when hk = 10’. 

(h) e*, to 4 decimals, from x = 0.8 to x = 0.9 when h = 0.01. 

(i) 3 (e — sin a), to 4 decimals (area of a segment of a circle subtending a central 
angle a, in radians) from a = 25° to a = 32° when h = 1°. 

2. Tabulate the differences for the following experimental results and indicate for 
each case the degree of the polynomial that would best express the relation between the 
variables. 

(a) S = stress in lbs. per sq. in. in steel wire used for winding guns, E = elongation 
in inches per inch. 


S 
E 


80,000 
0.00297 


60,000 
0.00216 


70,006 
0.00256 


50,000 
0.00173 


10,000 
0.00019 


20,000 
0.00057 


30,000 | 40,000 
0.00094 | 0.00134 


(b) Q = cu. ft. of water per sec. flowing over a Thomson gauge notch; H =ft. of head. 
Jak | 1 
Q 4.2 
(c) P/a = load in lbs. per sq. in. which causes the failure of long wrought-iron 


columns with round ends, //r = ratio of length of column to least radius of gyration of 
its cross-section. 


2.0 
14.9 


1.4 | 1.6 __ 1.8 
6.1 8.5 11.5 


l/r 140 18} = 220 260 300 340 380 420 
P/a| 12,800 7500 5000 3800 2800 2100 1700 1300 
(d) e = volts, p = kilowatts in a core-loss curve for an electric motor. 
e 40 | 60 | 80 100 120 140 160 
p O03 eal GO 2.18 3.00 3.93 6.22 8.59 


(e) A = amplitude of vibration in inches of a long pendulum, ¢ = time in min, 

since it was set swinging. 
t 0 I | 2 | 3 | 4 5 6 
4.97 2.47 Te22 0.61 | 0.30 | 0.14 


A 10 


(f) V = potential difference in volts, A = current in amperes in an electric circuit. 


A | 2.97 | 3-97 4.97 | 5.97 6.97 7:97 
V 65.0 61.0 58.25 56.25 55.1 54.3 
(g) 

=| - 1 s =| 1 13 
y 6.42 8.50 11.03 14.03 17.53 21.55 26.12 
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(hk) 
0.6 
1.27 


2.1 
0.25 


12 
0.63 


x 


y 


to) 
3.00 


0.3 
1.89 


0.9 
0.88 


24 | 2.7 
0.18 | 0.05 


1.5 | 1.8 
0.46 0.33 


3. By the method of differences explained in Art. 96, extend the tabulation of the 
functions in Exs. 1 a, b, d, e, h, i, for several values of the variables beyond the range 
of values for which the tables were constructed. 

4. Apply Newton's interpolation formula to the tables in Ex. 1. 

(a) In Ex. 1 a, find x? when x = 7.3 and x = 3.056. 

(b) In Ex. 1 b, find Wx when x = 566.2. 

(c) In Ex. 1 d, find rD?/6 when D = 1.452. 

(d) In Ex. 1 e, find log x when x = 361.4. 

(e) In Ex. 1 g, find log cos x when x = 88° 43’. 

5. Apply Newton’s interpolation formula to the tables in Ex. 2. 

(a) In Ex. 2 a, find EF when S = 42,000. 

(b) In Ex. 2 b, find Q when H = 1.7, and compare with the value given by the em- 
pirical formula Q = 2.672 H?-*. 

(c) In Ex. 2 c, find P/a when //r = 327, and compare with the value given by the 
empirical formula P/a = 417,000,000 (J/r)?-! 

(d) In Ex. 2f, find V when A = 4.07. 

(e) In Ex. 2 g, find y when x = 6. 

(f) In Ex. 2h, find y when x = 1.3 and x = 2.46. 

6. In the following table (taken from p. 129) 


0 288 293 313 | 333 
S 35.2 BE: 45.8 55-2 
S is the number of grams of anhydrous ammonium chloride which dissolved in 100 grams 
of water makes a saturated solution of 6° absolute temperature. Use Lagrange’s formula 
of interpolation to find S when @ = 300°, using (1) only two values of 6, (2) three values 
of 6, (3) all four values of 6. Compare the results with the value given by the empirical 
formula S = 0,000000882 63-. 
7. In the following table (taken from p. 141) 
4 2 | 4 | 8 
Ue 94 75 62 
4 is the current and V is the voltage consumed by a magnetite arc. Use Lagrange’s 
formula to find V when 7 = 3, and compare the result with the value given by the 
empirical formula V = 30.4 + 90.4 7-57, 
8. Use the methods of inverse interpolation (Art. 99) in the following: 
(a) In Ex. 1 a, find x when x? = 39 and when x? = 9.34. 
(b) In Ex. 1 e, find x when log x = 2.5542. 
(c) In Ex. 1 g, find x when log cos x = 8.3946 — I0. 
(d) In Ex. 2 a, find S when E = 0.00192. 
(e) In Ex. 2 c, find //r when P/a = 4000. 
(f) In Ex. 2 g, find « when y = 15.25. 
g. Approximate to the real roots of the equations: 
(a) #8 —2x+3=0. 
(b) ##—4x+2=0. 
(c) ee +22 —4=0. 
(d) rologx —x —2=0. 
(e) sinx + «2 —1.5 =0. 
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CHAPTER-IX; 
APPROXIMATE INTEGRATION AND DIFFERENTIATION. 


100. The necessity for approximate methods. —In a large number 
of engineering problems it is necessary to determine the value of the 


b 
definite integral, Af f(x) dx. Geometrically, this integral represents the 


area bounded by the curve y = f(x), the x-axis, and the ordinates x = a 
and x = b. Physically, it may represent the work done by an engine, 
the velocity acquired by a moving body, the pressure on an immersed 
surface, etc. If f(x) is analytically known, the above integral may be 
evaluated by the methods of the Integral Calculus. But if we merely 
know a set of values of f(x) for various values of x, or if the curve is drawn 
mechanically, e.g., an indicator diagram or oscillograph, or even where the 
function is analytically known but the integration cannot be performed 
by the elementary methods of the Integral Calculus — in all these cases, 
the integral must be evaluated by approximate methods — numerical, 
graphical, or mechanical. The planimeter is ordinarily used in measuring 
the area enclosed by an indicator diagram and in certain problems in 
Naval Architecture; such approximations often have the desired degree 
of accuracy. Where a higher degree of accuracy is required or where a 
planimeter is not available numerical methods must be used. 

In certain problems it becomes necessary to determine the value of 
the derivative, a Geometrically, this represents the slope of the curve 
y = f(x) at any point. Physically, it arises in problems in which the 
velocity and acceleration are to be found when the distance is given as a 
function of the time, in problems involving maximum and minimum 
values and rates of change of various physical quantities, etc. To 
evaluate the derivative we may use the methods of the Differential Cal- 
culus if the function is analytically known. Otherwise we are forced to 
use approximate methods — numerical, graphical, or mechanical. 

It is our purpose, in the following sections, to develop some of the 
numerical, graphical, and mechanical methods used in approximate in- 
tegration and differentiation. 

tor. Rectangular, Trapezoidal, Simpson’s, and Durand’s rules. — 
Suppose we wish to find the approximate area bounded by the curve 
y = f(x), the x-axis, and the ordinates x = x) and x = x, (Fig. 101). 

224 
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We divide the interval from x = x to x = x, into ” equal intervals of 
width h, and measure the (n + 1) ordinates yo, V1, Yo, . ~~ 5 Yn—1y Vn: 

(1) Rectangular rule. — If, starting at Po, we draw segments parallel 
to the x-axis through the points Po, Pi, Po, . . . , Pai, the area enclosed 
by the rectangles thus formed is given by 


Arp =h(yotmutyet::+ + Yn). 


Fic. Tor. 


If, starting at P,, we draw segments parallel to the x-axis through the 
points P,, Pr, ... , Pe, Pi, the area enclosed by the rectangles thus 
formed is given by 


Ap =hmtytos+ st oid + Yn) 


It is evident that the smaller the interval /, the better the approxima- 
tion to the required area. 

(2) Trapezoidal rule. — If the chords PoPi, PiP2, ... 5; PnP, are 
drawn, then the area enclosed by the trapezoids thus formed is 


ee n(% : 2) ns p(t ale p(eee um %) 4 tne p(t 2) 
=e so te Ve) eit ee Eval. 

This expression for the area is the average of the two expressions given 
by the rectangular rules. It is evident that the smaller the interval h 
and the flatter the curve, the better the approximation to the required 
area. If the curve is steep at either end or anywhere within the interval, 
the rule may be modified by subdividing the smaller interval into 2 
or 4 parts; thus, subdividing the steep interval between x,-; and x, in 
Fig. 101 


1 h n— ~Jm h m n 
into 2 parts: Ap = b(™5*) + te + 5 (tee) 4 (ee), 
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h (Yn + h (Ye + Ym 
into 4 parts: Ap = h(™2%) + te + 2 (ttm) 4 Ae 


h (9m + »*) (# + %). 
iter yet 

(3) Simpson’s rule. — Let us pass arcs of parabolas through the points 
P.PiP2, PoP3Ps, .« ., Pn-2Pn-iPn. Let the equation of the parabola 
through Py)P\P2 be y = ax* + bx +c. Then the area bounded by the 
parabola, the x-axis, and the ordinates x = x) and x = x, is 


= AG = Xo°) Se : (29? a Xo") 


Z2 2 Le 
A=] (ax? + bx + c) dx = e+e tee 


+ ¢ (x, — x) = ww) a oe [2 @ (x2? + 2X0 + X07) + 30 (x2+%0) + 6c]. 


Xe — Xo 


Now, Yo = ax? + bx +¢, n= am? + bate, h=— Tn 


ar? ton to=a(2t*) 49(2#ES)4¢ 


and we may easily verify that 
A= §h(w +4" + 42). 


If we have an even number of intervals and apply this formula to the 
successive areas under the parabolic arcs, we get 


Ag=3 h (yo+4 N+42) +3 h (Yo+4 Vs+4) + +s +3 An 24 YnttIn) 
=jh(ot4nt+2rm+4y3t+2y+ oe ae + 2 Yn-2 + 4 Yn—-1 + Va) 
=$§ Al (vot yn) +4 (M+ Yat s+ ++ + FVn—1) +2 (2+ Yet Yet - -- +Yn-2)]. 


To apply Simpson’s rule we must divide the interval into an even 
number of parts, and the required area is approximately equal to the 
sum of the extreme ordinates, plus four times the sum of the ordinates 
with odd subscripts, plus twice thé sum of the ordinates with even sub- 
scripts, all multiplied by one-third the common distance between the 
ordinates. 

(4) Durand’s rule.* — If we have an even number of parts and apply 
Simpson’s rule to the interval from x; to x,-1 and the Trapezoidal rule to 
the end intervals, 


A=hlSmtent+Gmtiv test - 
+2 3 Yn-3 + ¢ 3 Vn—2 a ry + (4 Yn +4 BUN IE 


Applying Simpson’s rule to the entire interval from xo to xn, 


A=h[3y +t emt+exwtiost - - + Yaz + F Yn—2 + FVn-1 + $F Ine 

Adding, 

2A =h[Fyt¥yn +2 +293 +... +2 yas +2 Yn-2 +42 ya + Syn. 
* Given by Prof. Durand in Engineering News, Jan., 1894. 
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Hence, 


Ap = h las (yo + Yn) +43 (1 + Yn) + 2 Fs t+ + + Yn-2] 
= =e ee cee (v1 + Yn—1) + 2+ ys t+ + + ne]. 


Collecting our rules, we have 


(1) Ap =AQot+s + e+ +++ + 4n-1), 
or Ar’ =h(ytyetyst e+ + Hy). - 
(2)0 Ae = hls oF yn) +t 2 + es Yn). 


(3) Ags = 310+ ye) £401 + 9s bos Es Vn) 
+2(2+yw+yte-: + Yn—2) J. 
(4) Ap = h [0.4 (yo + yn) + 1.1 (1 + Yn) + 2+ Ya e+ + ne). 


102. Applications of approximate rules. — We shall give some ex- 
amples illustrating the application of these rules. 


AG Cas ‘ ody 
1. Area. — Evaluate Hi) This is equivalent to finding ‘the area 
2 


between the curve y = I/x, the x-axis, and the ordinates x = 2 and 
x = 10. If we divide the interval into 8 parts, then # = 1; we have the 
table 


ej] $j 4] $j 9) 7) 4 | 9 |e 
Af 5 3 z 3 é 7 3 3 to 
ee aa tae te ty) = 1382905 
Ag =i ot tt a 4 Po) = 124200; 
Ap= 1G +r) +34 +++ +4) = 1.6290; 
Asa HG W) TAG TET ETD A 2G+E+D) = 161001 
Ap=1fo4¢+%)+11G64+4)+2+4+4+- + 3] = 1.6134. 
10 10 
By actual integration, [ oe Ine«| =Inio—Ing2= las = 1.6094. 
2 2 


We note that Simpson’s rule gives the best approximation (within 
0.1 % of the true value), with Durand’s next. 
If we take h = 3, 


Be Gt a)t sa; /2 


fe aie 43) 
+a(z area ares + £)| = 1.6096. 


Thus the Trapezoidal rule with 16 ordinates does not give the accuracy 
given by Simpson’s rule with 8 ordinates. 
2. Area. — The half-ordinates in feet of the mid-ship section of a 
vessel are 
Teh 2.01 2-0,4.2 4.0) 014-0; 12.0, 12:4, 110; 10.4, 6.8, O15; 
and the ordinates are 2 feet apart; find the area of the whole section. 


+34 $e++ +E] = 16144; 
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2 Ap = 2[3 (12.5 +0.5) + 12.8+ -->- + 6.8] = 224.8; 
4 Ag = 2[(12.5 + 0.5) ds Oso 1s Seo 
4+ 2 (12.9 + 13.0 + 12.4 + 10.4)] = 226.1. 
Hence, Ap = 449.6 sq. ft., Ag = 452.2 sq. ft. 
3. Work. — Given the following data for steam 
v 2 4 6 8 | 10 
pb | 68.7 31.3 19.8 14.3 LES 


where v is the volume in cu. ft. per pound and is the pressure in pounds 
per sq. in.; find the work done by the piston. 


10 
Work = a p dv; this is equivalent to finding the area under the curve 
2 


obtained by plotting (v, p). 
Wr = 2 [4 (68.7 + 11.3) + 31.3 + 19.8 + 14.3] = 210.80; 
Wg = 3 [(68.7 + 11.3) + 4 (31.3 + 14.3) + 2 (19.8)] = 201.33. 
By the methods of Chapter VI we find the empirical formula con- 
necting v and p to be pu! = 148, and hence, 


W= iP pdv = 148 f° y 12 dy = 148 


This last value differs from the value given by Simpson’s rule by about 1%. 
4. Mean effective pressure. Indicator diagram. Fig. 102a is a re- 
production of an indicator diagram; to find the mean effective pressure. 


yo: 12 


se vs = 199.31. 


0.89 
0.78 
0.64 


ee) 
be 
s 


0:977-~-~ \ 0-88 


i 
8 ase 
a = _ Ss 
i H ! 

| H | 

| H ! | 

i ! H 1 

1 | ! I 

' | 1 1 

! | ! 1 

' | ' | 

; | | | 

i ! | | 

I H I ! 

{ | | H 

| | 1 

| | 


Fic. 102a. 


The mean effective pressure P is the area of the diagram divided by 
the length of the diagram, since the area represents the effective area of 
the piston in sq. in. and the length represents the length of the stroke in ft. 
Since the total area enclosed by the curve is the difference between the 
area bounded by a horizontal axis, the end ordinates, and the upper part 
of the curve, and the area bounded by the same straight lines and the 
lower part of the curve, we need merely measure the lengths of the ordi- 
nates within the curve. The diagram is 3.5 ins. long. We divide the 
interval into 8 parts; then h = 7%, and we measure the ordinates 


0,  0.40,° 0.63; 0.91, 0.98; 1,00; > =-0.025" 30:74. 10; 
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Ar = ye [0.40 + 0.63 +-- - + + 0.74] = 2.44; 
Ag = 7s [4 (0.40 + 0.91 + 1.00 + 0.74) + 2 (0.63 + 0.98 + 0.92)] = 2.52. 
Hence, = As eect 0.72. 

ay |) he) 


We divide the interval into 14 parts; then # = }, and we measure the 
ordinates 


0, 0.30, 0.42, 0.54, 0.68, 0.88, 0.96, 0.98, 1.00, 1.02, 0.97, 0.89, 0.78, 0.64, O. 


Ar= {[0.30 + 0.42 + --- + 0.64] = 2.52. 
Ag=ys [4(0.30+0.54+ --- +0.64)+2](0.42+0.68+ +--+ +0.78)]=2.55. 
Hence, P= as eons 0:73: 

ot) oro 


We note that Ag with 9 ordinates has the same value as Ar with 15 
ordinates. 

5. Velocity. — Given a weight of 1000 
tons sliding down a 1% grade (Fig. 1020) 
with a frictional resistance of Io lbs. per ton 
at all speeds. The total resistance is 30,000 
Ibs. (a frictional resistance of 10,000 Ibs. and 
a grade resistance of 20,000 lbs.). Let the following table express the 
accelerated force F as a function of the time ¢ in seconds: 


Fic. 1020. 


t Co) 100 200 | 300 |400| 500 | 600 | 700 - 800 900 1000 
F |20,000| 19,000 Terao 5000] — 1000| — 5000] — 8500] — 11,000] — 13,000] — 15,000 
Find the velocity acquired by the body in 1000 seconds. 
: 2,000,000 ‘1,000,000 
Since F=mXa, and m= ———= ——; 
g 16.1 
F 16.1 dv 1000 
therefore, a@ = — = ——; and —=a, hence, v= a dt. 
m 1,000,000 dt 0 
We form a table for the acceleration a. 
gecall- 0 100 | 200 | 300 | 400 500 600 700 800 900 1000 
@ |0.322| 0.306! 0.258|0.177 | 0.081 | —0.016 | —0.081 | —0.137 | —0.177 | —0.209 |—0.242 
Here, i = 100, so that 
vp = 100 [4 (0.322 —0.242) + (0.306+0.258+ - - - —0.209)] = 24.2 ft. per sec. 


vg = 192 [(0.322 — 0.242) + 4 (0.306 + 0.177 — 0.016 — 0.137 — 0.209) 
+ 2 (0.258 + 0.081 — 0.081 — 0.177)] = 24.2 ft. per sec. 
6. Volume. —-- If S, is the area of a cross-section of a solid made by a 
plane perpendicular to the x-axis, then the volume of the solid included 


between the planes x» and x, is V = i "S,dx. In order to integrate, 
Xo 


we must know the analytical expression for S, as a function of x. 
Otherwise we employ the approximate formulas; the values of S, are the 
ordinates and h is the common distance between the cutting planes. 
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A buoy is in the form of a solid of revolution with its axis vertical, 
and D is the diameter in ft. at a depth p ft. below the surface of the water. 


p (0) 0.3 0.6 0.9 1.2, 1.5 1.8 
D 6.00 5.90 5.80 5,55 eS 4.70 4.20 
Dh 336.00 34.81 33.64 30.80 27.56 22.09 17.64 


Find the weight of water displaced by the buoy (1 cu. ft. of sea water 
weighs 64.11 lbs.). 


oof fil} 
Here, V= i yD ae: and h = 0.3, 
ae 


hence, Vg= ao [(36.00 + 17.64) + 4 (34.81 + 30.80 + 22.09) 


+ 2 (33.64 + 27.56)] = 41.38 cu. ft., 
and the weight of water displaced = 2652.87 Ibs. 
The areas in sq. ft. of the sections of a ship below the load-water plane 
and 3 ft. apart are’ 


7500, 7150, 6640, 5680, 4225, 2430, 260, 
where the load-water plane has an area of 7500 sq. ft. Find the dis- 
placement in tons (35 cu. ft. of sea water weigh I ton). 

Vr =3 [3 (7500+260) + (7150+6640+ 5680+ 4225 +2430) ]= 90,015 cu. ft. © 
‘Vs=3|(7500 +260) +4(7150+5680+2430) +2(6640+4225)] =90,530cu. ft. 


Hence, the displacement is 2572 tons by the Trapezoidal rule and 2587 
tons by Simpson’s rule. 
7. Moment of inertia. — The moments of inertia of an area about the 


axes are 
Sn e Ln 
Jz = 3 y° dx, Jy = xy dx. 
X Xo 


The evaluation of these integrals is equivalent to finding the areas under 
the curves with 4 y® or xy as ordinates and x as abscissas. 
The half-ordinates in ft. of the mid-ship section of a vessel are 


12.5). 12.8; (12.0, 13.0;° 13.0,/ 552.8, (12:45 910.5. 1040) 16.05 Oo 
and the ordinates are 2 ft. apart. Find the moment of inertia of the 
entire section about the axis. 

20 
Here, J, = 2 i 3 y°dx, h=2, and the values of 4% are 
1953-1, 2097 2, 2146.7, 2197.0, 2197.0, 2097.2, 1906.6, 1643.0, 1124.9, 314.4, 0.I, 
_ and applying Simpson’s rule, 


Jz = 3 (%) [(1953.1 + 0.1) + 4 (2097.2 + + + - + 314.4) 
+ 2 (2146.7 + +--+ + 1124.9)] = 22,266.1. 
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8. Pressure and center of pressure. — The pressure on a plane area 
perpendicular to the surface of the liquid, between depths xo and xn, is 


p=w fe xy dx, where w is the weight of the liquid per unit volume, y is 
the width of the area at a depth x beneath the surface. The depth of 


ff : xy aX 


| > xy dx 
Tq J 
these integrals can be evaluated approximately. 


9. Center of gravity. — The coérdinates of the center of gravity of an 
area are 


the center of pressure of such an area is given by x = All 


1 v2 
ae Je a _ Moment about OY ae fp ae _ Moment about OX | 


IF ne Area i WE Area 


The half-ordinates in ft. of the mid-ship section of a vessel are 


Deyo O et 20,7 13.0, 13:0," 12.5; 124, V 11.8," T0457 6:5,)_ 0.5, 


and the ordinates are 2 ft. apart. Find the center of gravity of the 
section. 


20 
ae sf ae _ Moment about OY 
4 *20 “i Area i 
y dx 


0 


and applying Simpson’s rule to the table, 


“| 0 2 4 6 8 10 12 14 | 16 18 20 
Nara S14 228 1 1250 |1 3:0) 13:0 | 12.8 012.4 | Li-S., 10,4) 6.8, | 0:5 
xy} O | 25.6] 51.6 |78.0/104.0] 128.0] 148.8 | 165.2 | 166.4 | 122.4 | 10.0 
Ms = 3[(o+ 10.0) +4 (25.6+ +--+ +122.4)+2(51.6+ --- +166.4)] 
= 2018.9. 
Ag = %[((72.5+0.5) +4(128+ --- +6.8)+2(12.9+ --- +10.4)] 
= 226.1. 
awe 2086,9 ie 
Hence, ara 8.93 ft. 


103. General formula for approximate integration. — We may derive 
a general formula for approximate integration by integrating any of the 
formulas of interpolation. Thus, Newton’s formula (p. 215), 
n(n —I1 n(n—1I)...(n—k+ 
Yn = Jot rag + 2 = Ny + Fag ( ) S = 


* 


1) i 
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where x = x) + mh, is true for all values of if some order of differences 
is constant or approximately constant. Multiplying by dn and inte- 
grating term by term between the limits o and n, we have 


n TA n n *n 
[yan = 0 f dn +5 {"ndn+% f n(n — 1) dr 
0 0 Jo [2 Jo 


+f [7 a@— 0 —2)an+ .s 


: x— x I 
Since x = x» + ah, therefore, 2 = mee and dz= ; dx. Hence, 


3 2 b 4 
{3 yds = b| ny + Say + ( - yea Cotamad le 


We 3p TL yee (F- 54.35% - 2" ‘\F vl 
ere nyse ee) 


Thus, if the differences after some order, as the kth, are negligible, we 
may use this formula to get the approximate area between the curve, the 
x-axis, and the ordinates x = x» and x = x,. The process is equivalent to 
approximating the equation of the curve by a polynomial of the kth 
degree. The differences ao, do, Go, . . . are those which occur in a line 
through yp parallel to the upper side of the triangle in the scheme on p. 210. 
Similar integration formulas can be derived from the other interpolation 
formulas. 

If the interval from xp to x, is large, it is well to divide this into smaller 
intervals, apply the formula to each of the smaller intervals, and add the 
results. In this way we may derive the formulas of Art. 1o1 and similar 
formulas as special cases of the above general formula. 

Let us first note that by means of the rule for the formation of the 
successive differences of a function (p. 210) we may express the differences 
G@o,-0o, Co; . ... in terms of yo, 41, yx; . . . ~- Thus, 


Oi = SAL SDDS 
by = a1 ~ A = (92 — M1) — (11 — Ho) = 32 —-2N +H, 

Co = by — bo = [(¥3 — 22 + 91) — (92 — 291 + 0) = ¥3 — 392 +3 V1 —Yo, 
dy = M4 — 493 +62 — 441 + Yo, 

&) = 5 — 534 + 10.43 — pacer 


k(k — 1) 


Ro = Yu — Rye + ae 


Sense 


where the coefficients in the right members of these equations are the 
binomial coefficients, taken alternately plus and minus. 
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(1) Let m =1 and bo, @, . . . all zero, .e., approximate the curve 
(Fig. 101a) from x» to x, by a straight line, y = A + Bx. Then 


fora = h [yo + 4 a] = h{[yo ts (yn — yo)] = | eee. 


Applying this result to each interval and adding, we get the Trapezoidal 
rule: 


Ar= ["yde=hUb Q0+ 9») ++ 92+ Pere Vert |: 


(2) Let 2 = 2 and ©, do, . . . all zero, 7.e., approximate the curve 
(Fig. 101a) from xo to x2, by a parabola, y = A + Bx + Cx®. Then 


Jo yd = lz 90 + 205 + 40] =h 1299+ 201 — 90) +3 02-2 1 +90] 
h 
ee: + ye]. 


Applying this result to an even number of intervals, two at a time, and 
adding, we get Simpson’s rule: 


As= ["y de ="[nb yn) +4019 oe py 1) H2(yot yt ++ + yn—2)]. 


(3) Let n = 3 and d,a, .. . all zero, 7.e., approximate the curve 
(Fig. 1o1a) from x to x3 by a parabola of the 3d degree, y = A+ Bx + 
Cx? + Dx’. Then 


ip ‘y dx = h[3 yo Got? bot 3 Co] = h[3 vot3-(yi—Y0) + $ 2-2 +40) 


+ 2 (93 — 392 +39 — ¥o)] = Fh [0 +39 +392 + a). 


Applying this result to 7 intervals, where 1 is a multiple of 3, and adding, 
we get Simpson's three-eighths rule: 


As'= [ “yde=Bhlu+2 (yet ver Verte 3 Gia vet vet rele 


(4) Let ~ = 6 and the differences beyond the 6th order negligible, 
1.€., approximate the curve (Fig. 1o1a) from xo to xs by a parabola of the 
6th degree, y = A+ Bx + Cx? + --- + Hx’. Then 


” y dee = h [6 yo + 18 dp + 27 by + 24 Co + 4%? do + $3 Co + Hes fol. 


TS xo 


Substituting the values of a, bo, . . . , fo in terms of the y’s and re- 
placing 745 fo by 3%5 fo, thus neglecting +49 fo which will be fairly small, we 
get Weddle’s rule: 


Ay = [“yde = tybbut5nt yet 69 ++ 5 6 + Hl 


We may apply this rule to 7 intervals where m 1s a multiple of 6. 
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2.6.x 
Example. Apply the approximate rules (1) to (4) to evaluate iy ae 


We divide the interval into 6 equal parts, so that h = 0.1. From the 
table 


x 2 2.1 2:2 2:3 2.4 2.5 2.6 
elena! i as ai as ee a3 
y 2 2.1 2.2 2.3 2.4 2.5 2.6 
1A me) I I I i sean 
dn=er[s Gag) tau baat a5 tea sl 0120244057 
o.1f /1 I I I PSV : 
4s= (E42) +4(5 +3 +5) +2(s5+ 5g) |= 02623644: 
Tg 2 zs ee ee oe) (e3s 
As eee aaa) 


3 
I 
+3 (=) +3, = 0.2623645; 


Ay =Son|3+5(7 : +3 : 5 +6(3 : i) + : poole =) +1, | =0.2623643. 


ee rest 


By integration, A= [4 Pes 6—In2=In 1.3 =0.2623637. 


2 


Ar agrees with A to 4 decimals, while As, As’, and Aw agree about 
equally well with A to 6 decimals. 

104. Numerical differentiation. — We are to find the slope of the 
curve y = f(x) at any point when the curve is drawn or a table of values 


eas : . ; d : 
of equidistant ordinates are given, 7.e., we are to find when the analyti- 


cal form of the function is unknown. Graphically, we must construct 
the tangent line to the curve at the given point. The exact or even 
approximate construction of the tangent line to a curve (except for the 
parabola) is difficult and inaccurate.* 


: , d : Ble ; 
We may derive an expression for = by differentiating Newton’s in- 


terpolation formula. Newton’s formula 


na wt nay tPF Dag... - RHR DOOR ED, 


is true for all values of m if some order of differences, as the kth, is con- 
stant or approximately constant. 

: “ es dy _idy dy _ 1 dy 
Since x = x) + nh, therefore, dx = hdn, and de hoa an 


* See Art. 106 on graphical differentiation. 
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Hence, ; 
oO - jot onde + Gn Ona) + (4m —18 wt +220— oat | 
ay 


dn ja| 2+ 2 — ot (On? — 18" +11) 2 + +} 


The values of these coefficients are tabulated for values of 2 between 
O and I at intervals of 0.01.* 

For the tabulated values x, 1, . . - , Xn, we have n = 0, so that for 
these values of x we have the simpler formulas 


I I I 
Tea Go ge tee — st foe ee i 
TEAS | II 
a = J | — + Toot i 


If the value of x for which = is required is near the end of the table, 


we may use similar formulas derived from the modified Newton’s formula 
for end-interpolation - 217). 

Example. Find 2 = and - > forx = 3andx = 3.3 from table (1) on 
p. 211 and check the ene by differentiating y = x’. 

Since x = 3 isa tabulated value we apply the second set of formulas: 


oN = t C7 Pap a ee 
= |37 5 (24) +26)| = 27; 7 = [24 — 6] = 18. 
From y=, = =3xH7 =27,, =5 = 6x = 18. 


For x = 3.3 we apply the first set of formulas, where a) = 37, bo = 24, 
(= 6,7 = 0.3; Then 


d 
a [37+ (- 0.4) t+ (one |= 32.67; $3 = [24+ (- 0.7)6|=198. 
dy d’y 
yh eae es ES ae a 
From A St os eas 3x7 = 32.67, FI 6x = 19.8. 


Example. Rate of change.—The following table gives the results of ob- 
servation; § is the observed temperature in degrees Centigrade of a vessel 
of cooling water, /is the time in minutes from the beginning of observation. 


t 0) | I | 2 | a | ee 5 
0 92.0 85.3 79.5 74.5 70.2 67.0 
To find the approximate rate of cooling when t = 1 and ¢ = 2.5. 


* See Rice, Theory and Practice of Interpolation. 
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From the table of differences 


t ) Al A? As 
(0) 92.0 
—6.7 
I 85.3 0.9 
— 5.5 Onl 
2 79-5 0.8 
—5-.0 —O.1 
3 74-5 0.7 
mAs 0.4 / 
4 70.2 1st 
—=3.2 
5 67.0 


when ¢=1, #=o0 and ae | -58 —< (08) +4 (-o.) |= — 6.23; 


when ¢ = 2.5, 7 = 0.5 and 7, = | -50 +o-+ (— 0.25) (+24) = 5.02. 


Example. Maximum and minimum.— The following table gives the 
results of measurements made on a magnetization curve of iron; B is the 
number of kilolines per sq. cm., uw is the permeability (Fig. 104). 


B 3 Sino LAS |S 13 


oO 4 Bescon 
yb 865 | 985 |1090|1175|1245|1295|1330 1120 


370 


2 
730 


I 
570 


10 
1340 


12 
1250 


II 
1320 


14 
930 


15 
725 


aS Bie Bie a Fe abe a YETI IC 


FIG. 104. 


To find the maximum permeability. In Fig. 104 the maximum perme- 
ability appears to be in the neighborhood of B = 10. We therefore tabu- 
late the differences of u in the neighborhood of B = to. 


B KL Al mx} A8 
9 | 1330 Ea 
10 I = 
340 — 20 3° —20 
II 1320 a eee ee 
12 1250 7 —60 e 


13 1120 soe 
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For values of B between B = 9 and B = 10, we have 


For a maximum, oe =0, hence 6”?+6n—11 =0, and n= 0.94. 
Therefore, B = By + nh = 9.94. 

We find the corresponding value of uw by the interpolation formula, 

= 1330 + (0.94) (10) + (0.0282) (—30) + (0.0100) (—20) = 1340. 

If we take account of A! and A? only, we get 
a = 10+ (2 — )(- 0) — 0, or n= : = 0.83, .and B =9.83. 
Then —_# = 1330 + (0.83) (10) + (0.0275) (—30) = 1337.5. 

105. Graphical integration. — Let us find the value of the definite 
integral a: ‘ f(x) dx or the area under the curve y = f(x) by graphical 


methods. Wedraw thecurve y = f(x) (Fig. 105a) and along the ordinate 
at P (x, y) erect the ordinate y’ whose value is a measure of the area under 


the curve y = f(x) from the initial point A (x = a) to the point P, 7.e., 
ay ez ic dx. Thus for every point P (x, y) we have a corresponding 


point P’ (x, y’). The curve traced by the point P’ (marked J in the 
figure) is called the integral curve and the curve traced by the point P 
(marked A in the figure) is called the derivative curve. Evidently, if P 
and Q are two points on the A-curve and P’ and Q’ are their correspond- 
ing points on the J-curve, the difference of the ordinates of P’ and Q’, 
y’’ — y’, is a measure of the area under the arc PQ. 

The practical construction of the integral curve consists of the follow- 
ing steps (Fig. 1050). , 

(1) Divide the interval from xp to x, into ” equal or unequal intervals 
and erect the ordinates yo, yi, .. - , Vn: 
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(2) Measure the areas %AoAim1 = yy’, xXoAoAoxe = yo’, ..., 
xoAvAnkn = Yn’. These areas may be found by means of a planimeter 
or by the construction of the mean ordinates. Thus, the area xAoAim 
is equal to the area of a rectangle whose base is xx; and whose altitude 
is the mean ordinate m, within that area. Similarly, the area 1,A;Aox1 
is equal to the area of a rectangle whose base is 1x2 and whose altitude is 
the mean ordinate m,. within that area. Estimate the mean ordinates 
M\, M2, M3, . . . , Mm, Within the successive sections. Then 
yy! Sint (xo%1), yo" = yy! + me (x1%2), ys" 7 yo" = Wile (eis) oes 

Yn = Yn —1 + Mn (Xn 1X0). 
If the intervals are all equal, 4.¢., xox, = “it, = ..- = Xn-1%a = AX, 
then y’ = XmAx. (We shall later give a more exact construction for the 
mean ordinate.) 


(3). :At [hy 03) te, eet 2 ee erecimon- 
dinates... 11Bi, -%2Boy 2, Anbae reama! 
respectively to y1’, yo’, ..-, Yn, and 
draw a smooth curve through the points 
Bo, Bi, Bs se 4 Ba Fis testectayve 


will approximate the required integral 
curve. 

Example. Construct the integral 
curve of the straight line y= 1— «x be- 
tween x = Ooandx = 2. (Fig. 105c.) 

Divide the interval from x = 0 tox =2 
into Io equal parts and erect the ordinates 
given in the table; here, Ax = 0.2. 


FIG. 105¢. 


: 2 5 8 3 ones : 18 
er a6 C7 0.14 ae a 
0.6 ols 28 ee o 2 
0.8 one 2:3 oo 0.48 
ax js Ont 0.02 0.56 
1.2| —o0.2 aw Dee 0.48 
io a Ons —o.06 oe 
1.6 —o.6 Se ens 0.32 
1.8 —o.8 7 ; 0.18 
2.0| —r.0 5 O58 =o ° 


It is evident that the mean ordinate in each section is merely one-half 
the sum of the end ordinates, so that the values of m are easily found. 
Erect the ordinates y’ and draw a smooth curve through the ends of the 


ordinates. The curve will approximate the parabola y’ = i ; (1 — x) dx 
0 - 


eT eae ei 
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Example. The following table gives the accelerations a of a body 
sliding down an inclined plane at various times ¢, in seconds. To find 
the velocity and distance traversed at any time, if the initial velocity 
and initial distance are zero. 


t 
a 


oO 
0.320 


100 
0.304 


200 
0.256 


300 
0.176 


400 
0.080 


500 600 
—0.016 | —0.080 


700 
——Ouns6) 


800 900 1000 
—0.176| —0.208} —0.240 


Since v = fe dtand s = fe dt, the time-velocity curve is the integral 


curve of the time-acceleration curve, and the time-distance curve is in 
turn the integral curve of the time-velocity curve. 
In Fig. 105d, we have plotted ¢ as abscissas and _a as ordinates. The 


units chosen are I in. = 100 sec., and I in. = 0.16 ft. per sec. per sec. 
70,000 
60,000 
sill ean ie 
Xs 
0.32 > == 100150,000 
0.16 ee 40,000 
(e (2) 
0 — “: 6.0|30,000 
@) 
@ 
~0.16 4 lay 40\20,000 
| fame irs ies =a =I 
20,32 20110,000 
bil zal | 
ae al 0} g 
0 100 200 300 400 iy Sade 700 800 900 1000 
Fic. 105d. 
t a AVQ. ACC. An A, Al Vv = Da,At | avg. vel. Up, Vm At S = Dp, At 
° ©.320 2 ° 6 6 ° 
Baa BT ne Tis 1560 
toe pees 0.280 28.0 ae 45.2 4520 1,560 
shes 0.256 0.216 216 59-2 70.0 7000 6,080 
300 0.176 : ; 80.8 ; 13,080 
0.128 12.8 Ouse 8720 
400 0.080 93.6 21,800 
= 0.032 Bee 68 Q5.2 9520 
500 0.016 ee: Be) 96. 94.4 9440 31,320 
eC 26. toh e408)" 2 36.6 8660 | 40,760 
700 | —O.136 7 81.2 49,420 
800 | —0.176 0-156 15.6 65.6 at as 56,760 
goo | —o0.208 Dis ee 46.4 ee eas 62,360 
1000 | —0.240 : ; 24.0 k 65,880 
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In each interval of 100 sec. we have estimated the mean acceleration 
as the average of the accelerations at the beginning and end of the in- 
terval; thus, in the first interval, am = } (0.320 + 0.304) = 0.312. This 
is equivalent to replacing the arcs of the curve by their chords or to find- 
ing the area by the trapezoidal rule. Since the initial velocity is zero, 
the (i, v) curve joins ¢ = 0, »v = 0 with ¢ = 100, v = 31.2, etc. We have 
drawn the (¢, v) curve with a unit of I in. = 20 ft. sec. 

In each interval of 100 sec. we have estimated the mean velocity as the 
average of the velocities at the beginning and end of the interval; thus in 
the first interval, vu = 4 (0 + 31.2) = 15.6. Since the initial distance 
is zero, the (¢, s) curve is drawn through the points ¢ = 0, s = 0,¢ = 100, 
Ss = 1560, etc. The unit chosen is I in. = 10,000 ft. 

The tables for v and s give the velocity and distance at the end of each 
100 seconds, and we may interpolate graphically or numerically for the 
velocity and distance at any time between ¢ = 0 and ¢ = 1000. 

In the foregoing discussion the accuracy of the construction of the 
integral curve depends largely upon the construction of the mean ordi- 
nates in the successive intervals. If the intervals are very small, we may 
get the required degree of accuracy by replacing 
the arcs by their chords and taking for the mean 
ordinate the average of the end ordinates. 

The approximation of the mean ordinate 
for the arc AjAi (Fig. 105e) is equivalent to 
finding a point M on the arc such that the area 
under the horizontal CoC; through M is equal 
to the area under the arc AoA; or such that the 
shaded areas ApC)M and A1CiM are equal. By means of a strip of cellu- 
loid and with a little practice the eye will find the position of M quite 
accurately, for the eye is very sensitive to differences in small areas. 

We may draw the integral curve by a purely graphical process. Let 
us first consider the case when the derivative curve is the straight line 
AB parallel to the x-axis (Fig. 105f). Choose a fixed point S at any con- 
venient distance a to the left of O. Extend AB to the point K on the y- 
axis anddraw SK. Through A’ (the projection of A on the x-axis) draw 
a line parallel to SK cutting the vertical through B in B’. Then, the 
oblique line A’B’ is the integral curve of the horizontal line AB. For, if 
P and P’ are two corresponding points, then 


FIG. 105e. 


YC AL OVE NG AO ROr ye == (y x A’Q) = ; X (area under AP). 


Similarly, for another horizontal CD, with C and B in the same verti- 
cal line, extend CD to the point L on the y-axis and draw SL; through B” 
draw a line parallel to SL cutting the vertical through D in C’’; then, the 
oblique line B’C’”’ is the integral curve of the horizontal CD. Finally, 


ART. 105 GRAPHICAL INTEGRATION 241 


draw B’C’ parallel to B’’C” or to SL; then the broken oblique line A’B’C’ 
is the integral curve of the broken horizontal line ABCD. 

Consider, now, any curve. Divide the interval from x to x, into n 
parts and erect the ordinates (Fig. 105g). Through Ao, Ai, Ay... , 
draw short horizontal lines. Cut the arc AoA, by a vertical line making 


M4 

SSF eas eee a ae te C D 

cee 10” 
“ : 

a Ae re eel 12 B nn 

ae eee, Yj ; / aii 
ce pee y, fed ae 
cogs [he Re Se co 
Biger jj as 
é Q a Q BY x 


Fic, 105f. 


the small areas bounded by this vertical, the arc, and the horizontals 
through A» and A;, equal. Proceed similarly for the succeeding arcs. 
Then construct the integral curve of the stepped line by the method 
explained above. Choose a point S at a convenient distance a to the left 
of O and join S with the points Co, Ci, C2, . . . , in which the extended 


FIG. 105g. 


horizontals cut the y-axis. Then, starting at Bo, draw a line through Bo 
parallel to SCp until it cuts the first vertical; through this point draw a 
line parallel to SC, until it cuts the second vertical, etc. The points 
where the resulting broken line cuts the ordinates at Ao, A1, Ao, ..., 
t.e., the points Bo, Bi, Bz, . . . , are points on the required integral curve; 
for at each of the points Ao, A1, A2, . . . , the area under the curve from 
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Ay to that point is equal to the area under the stepped line; so that a 
smooth curve through the points Bo, Bi, Bz, . . . will be the required 
integral curve. 

dy’ 


Pegs vk. 
Since y’ = E fy dx, therefore, ps 


integral curve at any point is proportional to the ordinate of the deriva- 
tive curve at the corresponding point. Furthermore, by the construc- 
tion, the slopes of the oblique lines through Bo, Bi, Bs, . . . are propor- 
tional to the ordinates yo, V1, y2, . . . , so that these oblique lines are 
tangent lines to the required integral curve at these points. We can 
thus get a more accurate construction of the integral curve by drawing 
the curve through the points Bo, Bi, Bs, . . . , tangent to the oblique 
lines through these points. 

The polar distance SO = a is constructed with the same scale unit as 
the abscissa x, and the ordinate y’ is measured with the same scale unit 
as the ordinate y. 

Example. Determination of the mean spherical candle-power of a mazda 
lamp. —In testing a lamp for the m. s.c. p., the intensity of illumination 
is measured every 15° by means of a rotating lamp and a photometer. 
The following table gives such measurements for a particular case: 

Angle 6° | o 15 1-301 45. | 000 275 90 | 105 | 120 | 135 | 150 | 165 | 180 
c-p 11.55 |13-0|15.4 | 22.4 |31.0|] 38.8 | 42.7 | 43.9 | 45.2 | 32.0] 21.8 | 9.1 | O 
According to the well-known Rousseau diagram, a semicircle is drawn 

(Fig. 105%) and divided into 15° sections, and perpendiculars are dropped 
from the points of division to the diameter, xo, x1, %2, ..., %2. Upon 
these perpendiculars the values of c-p are laid off as ordinates. The area 
under the curve AoAiA2 . . . Aig determined by these ordinates divided 
by the length of the base is the m. s.c. p. of the lamp, and this value 
multiplied by 4 7 will give the flux in lumens. 

To measure the required area we have constructed the integral curve 
(Fig. 1052) by the method described above. We chose 7 in. for the 
length of the diameter of the circleand I in. = 10c-f in laying off the ordi- 
nates. The y-axis or axis to which the horizontals are extended is drawn 
5 in. to the right of the point Ao, so that the polar distance is 4,0 = a 
= 5 in. 

The area under the curve ApAiA2 . . . Ai is measured by the ordi- 


= “yy, so that the slope of the 


: I 
nate X»Biz2 = 4.66. Since y’ = - Xx area, therefore area = a X y’ = 5 


X 4.66 = 23.3 sq. in. Since 1 in. on the scale of ordinates represents 

OF 2 16 
. . 7 

= 33.3 ¢-p. The straight line AoBy will cut the y-axis in a point D such 

that OD read on the c-p scale will also give the m. s. c. p., for 


10 c-p and the base of the diagram is 7 in., the m.s.c. p. = 
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mM. S.C. p. 


area _ 
base 


OD = X2Bie xa ps 
base 


33-0 c-p. 


We measure OD 


GG ee ee ee ee ee a 


pS 


N 


aA 


E. L. Clark 


LIS. 105h. 


Having drawn the integral curve we may immediately find the 


Thus, for 15° 


m.s.c. p. of any portion of the lamp between two sections. 
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on each side of the vertical, the m. s. c. p. isfound by drawing AoE parallel 
to B;B; and reading OE = 42.0 c-p on the candle-power scale, since 


OE a x7;B7 — sBs One area under A;A7 
a X5X7 base 


=mM.S.C. p. 


Similarly the m.s.c.p. of the section above a horizontal plane 
through the lamp is measured by OF = 37.0 c-p, and the m.s.c. p. of 
the section below a horizontal plane through the lamp is measured by 
OG = 29.5 c-p. 

106. Graphical differentiation. —If the integral curve y’ = f(x) is 


; eh dy’ : : 
given we may construct the derivative curve y = ae by using the prin- 


ciple that the ordinate of the derivative curve 
at any point P (x, y) (Fig. 106a) is equal to 
the slope of the integral curve or of the 
tangent line P’T at the corresponding point 
PNG 

The practical construction of the deriva- 
tive curve consists of the following steps: 
(1) Divide the interval from x to x, (Fig. 
1060) into parts and erect the ordinates 


FE . 6 . , 
Teo! Noy Va 9 Nel 9 << oy Da (2) GConstaietrine 
tangents at Bo, Bi, Bo, . . . , B, and measure their slopes. (3) At xo, 


1,2...) Xm erect ordinates xAy = Yo, 41A1 = uy 2. ) HaAn = Ynys Where 


Fic. 1060. 


the y’s are proportional to the corresponding slopes, and draw a smooth 
curve through the points Ao, Ai, A2,..., An. This curve will ap- 
proximate the required derivative curve. | 
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Example. The following table gives the pressure » in pounds per 
sq. in. of saturated steam at temperature 6° F. Construct the curve 
showing the rate of change of pressure with respect to the temperature, 


dp/do. 


> 
-) 


e 


CNH MO00 MAY 
ANAnNnANnnnn 
WWWWWAA LA A 
wWHhaNOOnN AN 
Le A oe Oe oe 


In the above table we have approximated dp/dé@ by Ap/AQ, 1.e., we 
have replaced the (0, p) curve by a series of chords, and the slopes of the 
tangents by the slopes of these chords. We then plotted (6, Ap/A@) and 
joined the points by a smooth curve (Fig. 106c). 


120 


170 


100 
() 
390 


Fic. 106c. 


It is evident that the difficulty in the construction of the derivative 
curve lies in the construction of the tangent line to the integral curve. 
The direction of the tangent line at any point is not very well defined by 
the curve. Asa rule it is better to draw a tangent of a given direction 
and then mark its point of contact than to mark a point of contact and 
then try to draw the tangent at this point. A strip of celluloid on the 
under side of which are 2 black dots about 2 m.m. apart may be moved 
over the paper so that the two dots coincide with points on the integral 
curve and so that the secant line which they determine is practically 
identical with the tangent line. If the arc AB (Fig. 106d) is approxi- 
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mately the arc of a parabola, we have a more accurate construction of the 
tangent; the line joining the middle points M@ and M’ of two parallel 
chords AB and A’B’ intersects the curve in P, the point of contact, 
and the tangent PT is parallel to the chord AB. 

We may also draw the derivative 
curve by purely graphical methods. 
The process is the reverse of the 
process described for constructing 
the integral curve (Art. 105). Let 
Bo, Bi, Bo, ... be the points of 
contact of tangent lines to the in- 
tegral curve (Fig. 105g). Choose a 
fixed point S at a convenient dis- 
tance a to the left of the y-axis and 
drawthe lines iS@6, Ci; Ge, ~ 2 os 
parallel respectively to the tangent lines at Bo, Bi, Bo, . ... Project 
the points Co, Ci, Cx, ..., horizontally on the ordinates at Bo, Bi, 
Bo, ..., cutting these ordinates in Ao, A1, Ao, .... The points 
Ao, Ai, Ao, ..., are then points on the required derivative curve, 
since ByAy + a = slope of SCy = slope of tangent at Bo, etc. We may 
now join the points Ao, A;, Az, . . . by a smooth curve, or we may get 
greater accuracy by using the stepped line of horizontals and verticals. 
Thus, we draw the horizontals through the points Ao, Ai, Ao, .. . , and 
the verticals through the points of intersection of consecutive tangents 
to the integral curve. The arcs AgA;, AiA2, ... , are now drawn so 
that the areas bounded by each arc, the horizontals, and the vertical, 
are equal. 

107. Mechanical integration.* The planimeter.— This is an in- 
strument for measuring areas. Consider a line PQ of fixed length / 
moving in any manner whatever in the plane of the paper. The motion 
of the line at any instant may be thought of as a motion of translation 
combined with a motion of rotation. Suppose the line PQ sweeps out 
the elementary area PQQ’P’ = dS (Fig. 1072). This may be broken up 
into a motion of translation of PQ to P’’Q’ and a motion of rotation 
from P’Q’ to P’Q’. If dn is the perpendicular distance between the 
parallel positions PQ and P’Q’ and d@¢ is the angle between P’’Q’ and 
P'Q’, then 


x 


Fic. 106d. 


dS =1dn +432 dg. 


* For descriptions and discussions of various mechanical integrators see: Abdanke 
Abakanowicz, Les Intégraphes (Paris, Gauthier-Villars); Henrici, Report on Planimeters 
(Brit. Assoc. Ann. Rep., 1894, p. 496); Shaw, Mechanical Integrators (Proc. Inst. Civ. 
Engs., 1885, p. 75); Instruments and Methods of Calculation (London, G. Bell & Sons); 
Dyck’s Catalogue; Morin’s Les Appareils d’Intégration. 
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Now if PQ carries a rolling wheel W, called the integrating wheel, whose 
axis is parallel to PQ (Fig. 1076), then, while PQ moves to the parallel 
position P’Q’, any point on the circumference of this wheel receives a 
displacement dn, and while P’’Q’ rotates to the position P’Q’, this point 
receives a displacement a dd, where a is the distance from Q to the plane 


FIG. 1074. Fic. 1070. 


of the wheel. So that, as PQ sweeps out the elementary area dS, any 
point on the circumference of the wheel receives a displacement 

ds = dn +adg. 
Therefore, dS =Ilds —aldo¢+4P do. 


Hence the total area swept out by PQ is 


S=1fas—al fas+4r fas. 


FIG. 107¢. 


Now, if PO comes back to its original position without turning com- 
pletely around, then the total angle of rotation [es = 0, so that 
Si hs: 


where s is the total displacement of any point on the circumference of the 
integrating wheel. 
But if PQ comes back to its original position after turning completely 


around, then 
S = Ils — 27ral + wi. 
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The most common type of planimeter is the Amsler polar planimeter * 
(Fig. 107c). Here, Fig. 107d, by means of a guiding arm OQ, called the 
polar arm, one end Q of the tracer arm PQ is constrained to move in a 
circle while the other end P is guided around a closed curve c-c-c- . 
which bounds the area to be measured. Then the area Q’/P’PP”Q” “ey 
is swept out twice but in opposite directions and the corresponding dis- 
placements of the integrating wheel cancel, so that the final displacement 
gives only the required area c-c-c- .... The circumference of the 
wheel is graduated so that one revolution corresponds to’a certain definite 
number of square units of area. 

p! 


Fic. 107d. 


The ordinary planimeter used for measuring indicator diagrams has 
l = 4 in. and the circumference of the wheel is 2.5 in.; hence one revolu- 
tion corresponds to 4 X 2.5 = 10sq.in. The wheel is graduated into Io 
parts, each of these parts again into Io parts, and a vernier scale allows 
us to divide each of the smaller divisions into Io parts, so that the area 
can be read to the nearest hundredth of a sq. in. The indicator diagram 
‘on p. 228 gives a planimeter reading of 2.55 sq. in., which agrees with the 
result found by Simpson’s rule with 15 ordinates. 

The polar planimeters used in the work in Naval Architecture usually 
have a tracer arm of length 8 in., and a wheel of circumference 2.5 in., so 
that one revolution corresponds “6 20 sq. in., thus giving a larger range 
for the tracing point. If the area to be racasereds is quite large, it may be 
split up into parts and the area of each part measured; or the area may 
be re-drawn on a smaller scale and the reading of the wheel multiplied 
by the area-scale of the drawing.t 


* This instrument was first put on the market by Amsler in 1854. 
+ If PQ (Fig. 107 d) turns completely around, the required area is S + 7 (OQ)?. 
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If very accurate results are required, account must be taken of several 
errors. (1) The axis of the integrating wheel may not be parallel to the 
tracer arm PQ. This error can be partly eliminated by taking the mean 
of two readings, one with the pole O to the left of the tracer arm, the other 
with the pole to the right * (Fig. 107e). ° 


This cannot be done with the ordinary 
Amsler planimeter because the tracer 
arm is mounted above the polar arm, 
but can be done with any of the Coradi 
or Ott compensation planimeters; one 
of these instruments is illustrated in 
Fig. 107f. (2) The integrating wheel 
may slip; some of this slipping may be 
due to the irregularities of the paper 


and has been obviated by the use of | 


disc planimeters, in which the recording 
wheel works on a revolving disc instead 
of on the surface of the paper. 


FIG. 107e. 


Various types of linear planimeters have been constructed. These 
differ from the polar planimeters in that one end of the tracer arm is 


= 
om 
] 
—— 


G. Coradi, Ziirich 


Fic. 107f. 


constrained to move in a straight line instead of in a circle. Planimeters 
of the linear type form part of the integrators described in Art. 108. 


Fic. 107g. 


Various other types of planimeters 
have been constructed, which do not 
have an integrating wheel. One of the 
best known of these is that of Prytz, 
also known as the hatchet planim- 
eter.t| In this form of the instrument 
(Fig. 107g) the end Q forms a knife- 
edge so that Q can only move freely 
along the line PQ. When P traces the 


given curve, Q will describe a curve such that PQ is always tangent to it. 


* For a proof of this statement, see Instruments and Methods of Calculation, p. 196. 
+ For the theory of this instrument, see F. W. Hill, Phil. Mag., xxxvili, 1894, p. 265. 
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Prytz starts the instrument with the point P approximately at the 
center of gravity G of the area to be measured, moves P along the radius 
vector to the curve, completely around the curve, and back along the 
same radius vector toG. The required area is then given approximately 
by 22¢, where / is the length PQ and ¢ is the angle between the initial 
and final positions of the line PQ. 

108. Integrators. — The Amsler integrator is practically an extension 
of the linear planimeter. In the latter instrument, the end Q of the tracer 
arm PO of constant length J, is constrained to move in a straight line X’X, 
while the tracing point P describes a circuit of the curve. If the axis of 
the integrating wheel attached to PQ makes a variable angle ma with 
X'X (Fig. 108a) at each instant, the point P will have for ordinate 


Vm = lsin ma, and the area described by P will be | lsin madx. On 


the other hand, the area described by P is equal to / times the displace- 
ment of any point on the circumference of the integrating wheel; hence 


He sin ma dx is equal to the displacement of a point on the circumference of 


an integrating wheel whose axis makes an angle ma with X'X. 


P 
2 ons 7 : 


Fic. 108a. Fic. 1080. 


Now, given a curve c-c-c- . . . (Fig. 1080), 


Area = fydx= [Isinade =1 fin ade. 
Moment of ar I - 2 
Sia ae tf tde=t [rsintady =< fo — cos 2 a) dx 
= {a = fs - d 
5 4 x fi sin (90° — 2 a) dx 


2 


l ; 
= —5 f sin (90° — 2a) dx, since dx = 0; 


the arm PQ returning to its original position when P makes a complete 
circuit of the curve. 


Moment of inertia 
of area about X/X _ £ [» de= fe sin’ a dx = = (3 sin a — zn 3 a)de 


=F fsinedr — 5 sing ade. 
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Now, f sin a dx, fein (90° — 2 a) dx, and fosin 3.adx, and hence 


the area, moment, and moment of inertia can be measured by three in- 
tegrating wheels whose axes at any instant make angles a, 90° — 2a, 
and 3 a, respectively, with X’X. 

The Amsler 3-wheel integrator (Fig. 108c) consists of an arm PQ and 
3 integrating wheels A, M,and J. The instrument is guided by a carriage 
which rolls in a straight groove in a steel bar; this bar may be set at a 
proper distance from the hinge of the tracer arm by the aid of trams. The 


axis of moments 


FIG. 108c. 


line X’X, which passes through the points of the trams and under the 
hinge, is the axis about which the moment and moment of inertia are 
measured. The radius of the disk containing the M-wheel is one-half 
the radius, and the radius of the disk containing the J-wheel is one-third 
the radius of the circular disk D to which they are geared. Therefore, 
the axis of the M-wheel turns through twice, and the axis of the J-wheel 
turns through three times the angle through which the tracer arm PQ or 
the axis of the A-wheel swings from the axis X’X. 

The integrating wheels are set so that in the initial position, 7.e., when 
‘P lies on X’X, the axes of the A- and J-wheels are parallel to X’X while 
the axis of the M-wheel is perpendicular to X’X. Then, when the tracer 
arm PQ makes an angle a with X’X, the axes of the A-, M-, and I-wheels 
make angles a, 90° — 2a, and 3a, respectively, with X’X. Further- 
more, the graduations of the M-wheel are marked so that these gradua- 
tions move backward while the graduations on the other wheels move 
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forward. Hence, when P has completed the circuit, and if a, m, and vare 
the displacements of points on the circumferences of the A-, M-, and I- 
wheels, respectively, we have 

P . Ip ee 

Area =la; Moment = rae Moment of Inertia = ae eee 
The wheels are graduated from 1 to Io so that a reading of 5, for 
example, means 5/10 of a revolution. The constants by which these 
readings are multiplied depend upon the length of the tracing arm and 
the circumferences of the integrating wheels. In the ordinary instru- 
ment, / = 8 in. and the circumferences of the A-, M-, and J-wheels are 


Ci = 2.5-in., Cy = 2.5 in., Cy = 2734375 10. 
Thus, to find the 


| 


area, a must be multiplied by 8 X 2.5 = 20; 

moment, Mi ‘ s st X 2.5 = 40; 

momentofinertia,a “ “ SS ce 2.5 = 320, 
2 


and 4 X 2.34375 = 100. 


12 
Finally, if a1, d2, my, m2, and 11, 72 are the initial and final readings of 
the A-, M-, and J-wheels, we have 


Area = 20 (dz — a); Moment = 40 (mz. — ™); 
Moment of Inertia= 320 (a2 — a1) — 100 (t2 — %). 


109. The integraph. — This is a machine which draws the integral 
curve, y’ = fre dx, of the curve y = f(x). The most familiar type of 


such machines is the one invented by Abdank-Abakanowicz in 1878. 
The theory of its construction is very simple. A diagram of the machine 
is given in Fig. 109a. The machine is set to travel along the base line of 
the curve to be integrated, and two non-slipping wheels, W, ensure that 
the motion continues along this axis. The scale-bar slides along the main 
frame as the tracing point P, at the end of the bar, describes the curve 
y = f(x) to be integrated. The radial-bar turns about the point Q which 
is at a constant distance a from the main frame. The motion of the re- 
cording pen at P; is always parallel to the plane of a small, sharp-edged, 
non-slipping wheel w, and by means of the parallel frame-work ABCD, 
the plane of the wheel w is maintained parallel to the radial bar [since w 
is set perpendicular to AB which is parallel and equal to CD throughout 
the motion, and the radial bar is set perpendicular to CD]. As the point 
P describes the curve y = f(x), the angle @ between the radial-bar and the 
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axis, and consequently the angle 0 between the plane of the wheel and the 
axis, are constantly changing, and the recording pen at P; draws a curve 
with ordinate y’ such that its slope 

dy’ 


> re a a 


CAND) 
te om 
and therefore, 
f I 
= x fs dx = q X area ORP, 


so that the curve drawn by P, is the integral curve of the curve traced 
by. 


Ww 


Main Frame 


os 


Ys R 


"mae 
il Keg, 
h 


1 Framework for 
. Parallel Motion D 


O 


: 
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If we now set the machine so that the point P traces the integral curve, 
then the recording pen P, will draw its integral curve 


r= forax = ff vax) ax ioe 


We may thus draw the successive integral curves y’, y’’, y’", .... Fig. 
1090 gives the integral curves connected with the curve of loads of the 
shaft of a Westinghouse-Rateau Turbine. The curve of loads is repre- 
sented by the broken line in the figure. By successive integration we get 
the shear curve, the bending moment curve, the slope curve, and the 
deflection curve. The distance marked “‘offset”’ is the distance OO, in 


Fig. 109a. 
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110. Mechanical differentiation. The Differentiator.— This is a 
machine which draws the derivative curve y’ = 2 of the curve y = f(x). 
Since the ordinate of the derivative curve is equal to the slope of the in- 
tegral curve, it is necessary to construct the tangent lines at a series of 
points of the integral curve. -We have already mentioned (Art. 106) the 
use of a strip of celluloid with two black dots on its under side to deter- 


U 


HIG TL; 


mine the direction of the tangent. This scheme is used in a differentiating 
machine constructed by J. E. Murray.* In a differentiating machine 
recently constructed by A. Elmendorf,{ a silver mirror is used for de- 
termining the tangent. The mirror is placed across the curve so that 
the curve and its image form a continuous unbroken line, for then the 
surface of the mirror will be exactly normal to the curve, and a perpen- 
dicular to this at the point of intersection of the mirror and the curve will 
give the direction of the tangent line. If the surface of the mirror de- 


* Proc. Roy. Soc. of Edinburgh, May, 1904. 
+ Scientific American Supplement, Feb. 12, 1916. ' 
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viates even slightly from the normal, a break will occur at the point where 
the image and curve join. It is claimed that with a little practice a re- 
markable degree of accuracy can be obtained in setting the mirror. 

Fig. 110 gives a diagram illustrating the working of this machine. 
The tracing point P follows the curve y = f(x) so that the curve and its 
image in the mirror MP form a continuous unbroken line; then the arm 
PT, which is set perpendicular to the mirror, will take the direction of the 
tangent line to the curve. The link PR, of fixed length a, is free to move 
horizontally in the slot X’X’ of the carriage C. The vertical bar SU 
passes through R and is constrained to move horizontally by heavy rollers. 
The point Q slides out along the tangent bar PT and also vertically in the 
bar SU, carrying with it the bar QP’. If we choose for the x-axis a 
line XX whose distance from X’X’ is equal to QP’, then the point P’ 
will draw a curve whose ordinate is equal to y’ = RQ. But RQ/a is 


the slope of the tangent PT, hence, y’ = a X a and the curve drawn 


by P’ is the derivative curve of the curve traced by P. 

The machine is especially useful for differentiating deflection-time 
curves to obtain velocity-time curves, and by a second differentiation, 
acceleration-time curves. It is also helpful in solving many other 
problems. 


EXERCISES. 
Apply the approximate rules of integration (Art. 101) to the following examples: 


1.0 dx , 
1. Evaluate ( = when h = 0.1, and when h = 0.05, and compare the results with 


the values obtained by integration. 
2. Evaluate ie sin x dx, when h = g and when h= = and compare the results 


with the values obtained by integration. 


3. The arc of a quadrant of an ellipse whose eccentricity is 0.5 is given by 
us 


Fee ae ara 
if Vi — 0.25 sin? ¢ dg. Evaluate the integral when h = 9°. 


3 
4. Evaluate if : a , when hk = 0.5. 


0 Ves —x +1 
5. The semi-ordinates in ft. of the deck plan of a ship are 
3; 10:6,5, 92525, 28:0,08 (20.8 9mn sonmezo.e, 29°5,, 928.50) 24.2 eel Se 
these measurements are 28 ft. apart. Find the area of the deck. 
6. Given the following data for superheated steam 
v | 2 | 4 | 6 
p 105 42.7 25.3 
Find the work done. 


8 | 10"%. 
16.7 13 
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7. The length of an indicator diagram is 3.6 in. The widths of the diagram, 0.3 in. 
apart, are 


OMO:4 O-==- O15 20,0380! 210056 0.005 1200.8 1.00; 1.00) 11.00}10.97; 0; 
Find the mean effective pressure. 


8. The length of an indicator diagram is 3.2 in. The widths of the diagram, 0.2 in. 
apart, are 


1.00, 1.68, 1.62; 1.00, 0.64, 0.48, 0.36, 0.26, 0. 
Find the mean effective pressure. 


9. The speed of a car is v miles per hour at a time ¢ seconds from rest; 
t fo) 


5 10 15 20 | 25 30 
v Co) 207, | Piet 10.9 LA TOMm nS) 14 
Find the distance traversed in 30 seconds. 
10. s is the specific heat of a body at temperature 6° C. 
6 oO 2 4 6 8 | 10 12 
5 1.00664 1.00543 1.00435 1.00331 1.06233 1.00149 1.00078 


Find the total heat required to raise the temperature of a gram of water from 0° C. to 


85 
12-1; (total heat-= f s d6). 


L 
11. The areas in sq. ft. of the sections of a ship above the keel and two feet apart are 
2690, 3635; 4320, 4900, 5400. 
Find the total displacement in tons. 


12. A reservoir is in the form of a volume of revolution and D is the diameter in ft. 
at a depth of p feet beneath the surface of the water. 


p 16 32 | 48 | 64 80 2 
D 105 100 86 | 66 48 27 


Find the number of gallons of water the reservoir holds when full. 


Oo 
IIo 


13. A plane board is immersed vertically in water. The widths of the board in ft. 
parallel to the surface of the water and at depths 3 ft. apart are 


4.0, 3.6, 3.0, L.7% 13; Oy, 0.8, 0.6, 0.1. 


Find the pressure on the board and the depth of the center of pressure when the surface 
of the water is level with the top of the board. 


14. The half-ordinates in ft. of the mid-ship section of a vessel at intervals 2 ft. 
apart are 


1D WRG OS, TAR GR TOR, ao ios SKOrn ey les CO 
Find the position of the center of gravity of the section. 


15. The shape of a quarter-section of a hollow pillar is given by the following table. 
The axes of x and y are the shortest and longest diameters. 


xin. , O |0.25 |0.50 |0.75 |1.00 |1.25 |I 50|1.75 | 2.00 | 2.25 | 2.50 | 2.75 | 3.00 | 3.25 |3.50 


side | 6 [5.95 |5-90 |5-83 |5.76 |5-64 |5.48 |5.22 | 4.99 4.68 | 4.35 | 3.88 | 3.25 | 2.34] 0 


side | 5 |4.90 /4.78 |4.65 |4.45 |4.22 |3-80 |3.40 | 2.77 | 2.08 | 0 


Find the moments of inertia of the section about the x- and y- axes. 
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16. Apply the formulas for numerical differentiation (p. 235) to table (2) y = x* on 
p. 211, and find a and © when x = 5.31 and x = 5.33. Check the results by actual 


differentiation. 


17. Apply the formulas for numerical differentiation (p. 235) to table (8) y = log sin x 
dy 
dx 
actual differentiation. 


2. 
on p. 212, and find and > when x = 1°20’ and x = 1°24’. Check the results by 


18. Inthe following table, s is the distance in ft. which the projectile of a gun travels 
along the bore in ¢ sec. 


s| o I 2 3 4 5 Gin if | 8 9 10 

t| 0 | 0.0360 | 0.0490] 0.0598] 0.0695 | 0.0785 | 0.0871 | 0.0953 | 0.1032 | 0.1109 | 0.1184 
; ; _ds__ /dt : _ Oss) tty. =) e 

Find the velocity v ar Pa 5) ds? and the acceleration a = fe ia) ee when s = 


iB ili 


19. Use the data given in Ex. 6 to find the rate of change of the pressure with re- 
spect to the volume, dp/dv, when v = 4 and v = 5. 
dv 


20. Use the data given in Ex. 9 to find the acceleration, a = a’ when ¢ = 10 and 


t = 12. 
21. Find the minimum value of the polynomial which has the values 
€ (o) 2 4 6 
ay 3 | 3 II | 27 
22. The following table gives the results of measurements made on a normal in- 
duction curve for transformer steel; B is the number of kilolines per sq. cm.; u is the per- 
meability. 
B 


u 


6 
1465 


8 
1480 


I 
625 


2 
870 


10 
1370 


ai 
1280 


12 
1130 


9 
1430 


i 
1520 


4 
1210 


3 5 
1035 1350 
Find the maximum permeability. 


23. Construct the integral curve of the parabola y = x — 4 x2 asx varies from 0 to 2. 
24. Construct the integral curve of the sine wave y = 2 sin 2 x as x varies fromo tom. 


25. The following table gives the accelerations a of a body sliding down an inclined 
plane for various distances s in ft. 


100 
0.304 


Use the method employed in the illustrative example on p. 239 for drawing the integral 


200 
0.256 


600 
—0.080 


300 
0.176 


400 
0.080 


500 
—0.016 


1000 
—0.240 


y (e) 
a '0.320 


700 | 800 goo 
—0.136| —0.176 | —0.208 


curves and determining the velocity, v = \/2 fa ds, and the time, ¢ = f- ds, for any 
distance, if v = o and t = o whens = 0. 


26. The following table gives the accelerations a of a body at various velocities v in 
ft. per sec. 
I 
0.360 


2 
0.283 


3 
0.179 


4 
0.069 


5 
0.013 


v | (0) 
a 0.405 


: . . I : 
Draw the integral curves to determine the time, tf = f- dv, and the distance, s =f- dv, 
a 


for any velocity, if tf = o and s = 0 whenv = 0. 
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27. In the following table 


oO 4 
38,000 38,500 


Ss 
P 


I 
38,500 


8 
27,500 


6 15 20 
35,500 15,700 3850 


P is the resultant pressure in pounds on the piston of a steam engine at distances s inches 
from the beginning of the stroke. Draw the integral curve to find the work done as the 


piston moves forward. (Work = fe as.) 


ie5 
19,000 


19 
11,000 


28. Acar weighs 10 tons. It is drawn by a pull of P lbs.; ¢ is the time in seconds 
since starting. 
t (0) 
di 1020 


2 5 
980 882 


8 10 
720 702 


13 | 16 19 22 
650 713 YEP | 805 
If the retarding friction is constant and equal to 410 lbs., draw the integral curve to find 


the speed of the car at any time. (Momentum ={( — 410) at.) 


29. In the following table 


t 
v 


0.00490 
869 


0.00785 
1142 


0.00953 
1242 


0.00598 | 0.00695 0.00871 0.01032 err 0.01184 
987 1074 1195 1277 1309 1335 
v is the velocity of projection in ft. per sec. in the bore of a gun at time # sec. from the 


beginning of the explosion. If s = 2 ft. when ¢ = 0.00490 sec., draw.the integral curve 
to show the relation between the distance and the time. 


30. A beam 10 ft. long is loaded as in the following table, where w is the weight per 
unit length at distances x ft. from the free end. 


2 | 3 | 4 5 
Sp a 55 eT ony 


ie) 
2 


ye 
Des 


6 
II.2 


x 
WwW 


8 
rng} 


ei 
12.2 


9 | 10 
10.2 7.2 


Draw integral curves to show (1) the shearing force, s = fe dx and (2) the bending 


moment, M =fs dx. 


31. The following table gives the measurements for every 15° of the intensity of 
illumination of a lamp. 
Angle é@°| o | 15 (30 

c-p 60.5 ;88.0| 99.5 


Apply the method of the illustrative example on p. 242 to find the m.s.c.p. for various 
sections of the lamp. 


120 
28.0 


180 
12.5 


45 | 60 
86.5 | 50.0 


75 
31.0 


90 
29.0 


105 
29.0 


165 
13.0 


135 
20.0 


150 
15.0 


32. In the following table 
Ce) 10 20 30 
O 156 608 | 1308 


t 
§ 


40 | 50 60 70 | 80 | 90 | 100 
2180 | 3132 | 4076 | 4942 | 5676 | 6236 | 6588 
s is the distance in ft. traversed by a body weighing 2000 Ibs. in¢ sec. Draw the deriva- 
tive curves to show the velocity and acceleration at any time. Also draw the curve 
showing the relation between the kinetic energy and the force. 


33. The observed temperature @ in degrees Centigrade of a vessel of cooling water 
at time ¢ in min. from the beginning of observation are given in the following table: 


t | (o) I 2 3 | 5 7 10 
9 92.0 85.3 79:5 74-5 67.0 60.5 53-5 


Draw the derivative curve to show the rate of cooling at any time. 


15 | 20 
45.0 39-5 


Adiabatic expansion formula, 48 
chart for, 33, 49 
Alignment or nomographic charts (see 
also Charts, alignment or nomo- 
graphic) 
fundamental principle of, 44 
with curved scales, 106 
with four or more parallel scales, 55 
with parallel or perpendicular index 
lines, 87, 91, 97 
with three or more concurrent scales, 104 
with three parallel scales, 45 
with two intersecting index lines, 68 
with two or more intersecting index 
lines, 76 
with two parallel scales and one inter- 
secting scale, 65 
Approximate differentiation, 224 
Approximate integration, 224 
Area, 
by approximate integration rules, 227 
by planimeter, 246 
Armature or field winding formula, 90 
chart for, 90 


Bazin formula, 101 
chart for, 102, 116 


Center of gravity, by approximate inte- 
gration rules, 231 
Chart, alignment or nomographic, for 
adiabatic expansion, 49 
armature or field winding, 90 
Bazin formula, 102, 116 
Chezy formula for flow of water, 58 
D’Arcy’s formula for flow of steam, 81 
deflection of beams, 72, 73, 86 
discharge of gas through an orifice, 89 
distributed load on a wooden beam, 83 
focal length of a lens, 106 
Francis formula for a contracted weir, 
109 
friction loss in pipes, 94 


Chart, Grasshoff’s formula, 51 
Hazen-Williams formula, 60 
horsepower of belting, 54 
indicated horsepower of a steam en- 
gine, 63 

Lamé formula for thick hollow cylin- 
ders, 92 

McMath “run-off,’’ formula, 49 

moment of inertia of cylinder, ree 

multiplication and division, 47 

prony brake, 70 

resistance of riveted steel plate, 103 

solution of quadratic and cubic equa- 
tions, I12 

specific speed of turbine and water 
wheel, 75 

storm water run-off formula, 108 

tension in belts, 54 

tension on bolts, 67 

twisting moment in a cylindrical shaft, 
78 

volume of circular cylinder, 49 

volume of sphere, 49 

Charts, hexagonal, 40 

Chart with network of scales, for 
adiabatic expansion, 33 
chimney draft, 38 
elastic limit of rivet steel, 34 
equations in three variables, 28 
equations in two variables, 20 
multiplication and division, 30, 31 
solution of cubic equation, 36 
temperature difference, 39 

Chezy formula for flow of water, 56 
chart for, 58 

Chimney draft formula, 37 
chart for, 38 


Coefficients in trigonometric series evalu- 


ated, 
by six-ordinate scheme, 179 
by twelve-ordinate scheme, 181 
by twenty-four-ordinate scheme, 185 
for even and odd harmonics, 179 


xii 
Coefficients in trigonometric series evalu- 
ated, 
for odd harmonics only, 186 
for odd harmonics up to the fifth, 187 
for odd harmonics up to the eleventh, 
189 
for odd harmonics up to the seventeenth, 
IgI 
graphically, 200 
mechanically, 203 
numerically, 179, 186, 192, 198 
Constants in empirical formulas deter- 
mined by 
method of averages, 124, 126 
method of least squares, 124, 127 
method of selected points, 124, 125 
Coordinate paper, 
logarithmic, 22 
rectangular, 21 
semilogarithmic, 24 


D’Arcy’s formula for flow of steam, 79 
chart for, 81 

Deflection of beams, 70, 71, 84 
chart for, 72, 73, 86 

Differences, 210 

Differentiation, approximate, 224 
graphical, 244 
mechanical, 255 
numerical, 234 

Differentiator, 255 

Discharge of gas through an orifice, 89 
chart for, 89 

Distributed load on a wooden beam, 80 
chart for, 83 

Durand’s rule, 226 


Elastic limit of rivet steel, 32 
chart for, 34 
Empirical formulas, 
determination of constants in, 124, 125, 
173, 174 
for non-periodic curves, 120 
for periodic curves, 170 
involving 2 constants, 128 
involving 3 constants, 140 
involving 4 or more constants, 152 
Equations, solutions of (see Solutions of 
algebraic equations) 
Experimental data, 120, 170 
Exponential curves, 131, 142, 151, 153, 
156 


INDEX 


Focal length of a lens, 
chart for, 35, 40, 106 
slide rule for, 15 
Fourier’s series, 170 
Francis formula for a contracted weir, I10 
chart for, 109 
Friction loss in pipes, 94 
chart for, 94 
Fundamental of trigonometric series, 170 


Gauss’s interpolation formula, 219 
Graphical differentiation, 244 
Graphical evaluation of coefficients, 200 
Graphical integration, 237 
Graphical interpolation, 209 
Grasshoff’s formula, 50 

chart for, 51 


Harmonic analyzers, 203 
Harmonics of trigonometric series, 170 
Hazen-Williams formula, 57 
chart for, 60 
Hexagonal charts, 40 
Horsepower of belting, 52 
chart for, 54 
Hyperbola, 149 
Hyperbolic curves, 128, 135, 137, 140 


Index line, 44 : 
Indicated horsepower of steam engine, 61 
chart for, 63 
Integraph, 252 
Integration, approximate, 224 
applications of, 227 
by Durand’s rule, 226 
by rectangular rule, 223 
by Simpson’s rule, 226, 233 
by trapezoidal rule, 225 
by Weddle’s rule, 233 
general formula for, 231 
graphical, 237 
mechanical, 246 
Integrators, 250 
Interpolation, 209 
Gauss’s formula for, 219 
graphical, 209 
inverse, 219 
Lagrange’s formula for, 218 
Newton’s formula for, 214, 217 : 
Isopleth, 44 


Lagrange’s interpolation formula, 218 
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Lamé formula for thick hollow cylinders, 
gI 
chart for, 92 
Least Squares, method of, 124, 127 
Logarithmic codrdinate paper, 22 
Logarithmic curve, I51 
Logarithmic scale, 2 


Maxima and minima by approximate 
differentiation formulas, 236 
McMath ‘‘run-off”’ formula, 48 
chart for, 49 
Mean effective pressure by approximate 
integration rules, 228 
Mechanical differentiation, 255 
Mechanical integration, 246 
Moment, by integrator, 250 
Moment of inertia, 
by approximate integration rules, 230 
by integrator, 250 
Moment of inertia of cylinder, 99 
chart for, 100 
Multiplication and division, charts for, 
30, 31, 41, 47 


Newton’s interpolation formula, 214, 217 
Nomographic or alignment charts (see 
Alignment or nomographic charts) 
Numerical evaluation of coefficients, 179, 

186, 192, 198 
Numerical differentiation, 234 
Numerical integration, 224 
Numerical interpolation, 215 


Parabola, 145 
Parabolic curves, 128, 135, 140 
Periodic phenomena, representation of, 
170 
Planimeter, 
Amsler polar, 248 
compensation, 249 
linear, 249 
principle of, 246 
Polynomial, 159 
Pressure and center of pressure, by 
approximate integration rules, 231 
Prony brake, 69 
chart for, 70 


Rates of change, by approximate differ- 
entiation formulas, 235 


- 


Rectangular coérdinate paper, 21 
Rectangular rule, 225 

Resistance of riveted steel plate, 101 
’ chart for, 103 


Scale, 
definition of, 1 
equation of, 2 
logarithmic, 2 
representation of, I 
Scale modulus, 2 
Scales, 
network of, 20 
perpendicular, 20 
sliding, 7 
stationary, 5 
Semilogarithmic codrdinate paper, 24 
Simpson’s rule, 226, 233 
Slide rule, 
circular, 16 
for electrical resistances, I5 
for focal length of lens, 15 
Lilly’s spiral, 18 
logarithmic, 9 
log-log, 13 
Sexton’s omnimetre, 17 
Thacher’s cylindrical, 18 
Solutions of algebraic equations, 
by means of parabola and circle, 26 
by means of rectangular chart, 35 
by means of alignment chart, I10 
by method of inverse interpolation, 221 
on the logarithmic slide rule, 11 
Specific speed of turbine and water wheel, 
73 
chart for, 75 
Storm water run-off formula, 107 
chart for, 108 i 
Straight line, 122, 125 


Tables, construction of, 213 
Temperature difference, 37 
chart for, 39 
Tension in belts, 52 
chart for, 54 
Tension on bolts, 66 
chart for, 67 
Trapezoidal rule, 225 
Trigonometric series, 170 
determination of constants in, 173, 174 
Twisting moment in a cylindrical shaft, 77 
chart for, 78 
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Volume of sphere, 50 


Velocity, by approximate integration rules, 
chart for, 49 


229 
Volume, by approximate integration rules, 
229 
Volume of circular cylinder, 48 
chart for, 49 


Weddle’s rule, 233 
Work, by approximate integration rules, 
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